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Beitrag zur Ermittlung der hydrodynamischen Dimpfungs- 
. und Federeigenschaften von Gleitlagern. 


Von E. Pestel. 


1, Einleitung. Die Kenntnis der hydrodynamischen Dimpfungs- und Federeigenschaften 
von Gleitlagern ist wiinschenswert fiir die genaue Ermittlung der kritischen Drehzahlen einer 
in Gleitlagern laufenden Maschinenwelle und fiir die Voraussage ihres Verhaltens in der Um- 
gebung einer solchen Drehzah]. Dariiber hinaus ist sie aber eine notwendige Voraussetzung fiir 
die Beurteilung der Stabilitaét der rotierenden Welle gegen selbsterregte Schwingungen. Das 
Auftreten kritischer Wellenstérungen infolge der Nachgiebigkeit des Schmiermittels in Gleit- 
lagern wurde zuerst von Newkirk} und von Stodola? erkannt. Letzterer wies auch einen Weg 
zur Ermittlung der Federeigenschaften des sogenannten Halblagers. Unter Vernachlassigung 
der Diampfungseigenschaften dieses Lagers stellte Stodola Stabilitatskriterien fiir die ,,starre“ 
und spater auch fiir die elastische rotierende Welle auf. Die Arbeiten von Newkirk und 
Stodola wurden vor einigen Jahren von Hagg ® aufgegriffen mit dem Ziel, der Praxis zuver- 
lassigere Unterlagen fiir die rechnerische Vorausbestimmung der Laufstabilitat zu liefern. 

Da die zahlenmaBige Kenntnis der Dampfungseigenschaften allein fiir die Aufstellung des 
Routhschen Stabilitatskriteriums * unzureichend ist, wird am SchluB dieser Mitteilung auch 
das Verfahren zur Bestimmung der Federzahlen des sogenannten 120°-Lagers erlautert, 
das im wesentlichen auf einer experimentellen Arbeit von Needs und deren Auswertung 
durch Hagg unter Benutzung des Grundgedankens der Methode von Stodola beruht. 

Die Betrachtung des in einem Gleitlager rotierenden Zapfens 
einer als starr angenommenen, ideal ausgewuchteten, zweifach- — Ghiohpeminhislyeas|__raumfeskes 
gelagerten Welle mit waagerechter Achse fihrt wohl am ein- “*#ilefpunktes Leila 
fachsten in die dieser Arbeit zugrundeliegende Fragestellung ein. 
Infolge der Rotation des Zapfens in einem viskosen Schmiermittel 
kann auch die ,,starre‘‘ Welle gedampfte Higenschwingungen 
ausfiihren. Bei translatorischer Bewegung der Wellenachse 
kénnen wir die Masse der Welle auf die beiden Zapfen ent- 
sprechend den Auflagerdriicken verteilen und fiir die Bewegung 
des Zapfenmittelpunktes, in dem wir uns den entsprechenden 
Massenanteil m konzentriert vorstellen, folgende Differential- 
gleichungen aufstellen (Abb. 1). 


M% Be % Dee ZH Cee XH + Oe, Z=— 0, (la) 
mz+b,*+b,z+og,x*+¢,2=0. (1b) 
Darin sind die DampfungsgréBen b und die Federzahlen c abhaingig von der dimensionslosen 
Lagerkennzahl G = VoE , welche die verschiedenen Einfliisse des radialen Lagerspiels 6, 


des mittleren Lagerdruckes p, der Zahigkeit ~ des Schmiermittels und der Drehzahl n/min 
zusammenfakt. Die Gleichgewichtslage des rotierenden Zapfens, gekennzeichnet durch die 


1 B. L. Newkirk u. H. D. Taylor, General Electric Review 28 (1925), S. 559. 

2 A. Stodola, Schweiz. Bauztg. 85 (1925), S. 265. 

8 A. C. Hagg, Trans. A.S.M.E. (1946), A — 211—220. Der Verfasser verdankt A, C. Hagg die An- 
regung zu der im folgenden abgeleiteten Methode fiir die Berechnung der Dampfungseigenschaften von 
Gleitlagern. 

4 Weer ctaie fiir die Anwendung dieses Kriteriums ist natiirlich die Konstanz der Groen b und c 
in den Differentialgleichungen (1). Im Maschinenbau wird nun bereits das Einsetzen instabilen Verhaltens 
als unerwinscht betrachtet. Fiir die Bestimmung der Stabilitatsgrenze ist dann aber die Annahme kon- 
stanter b und cin den meisten praktischen Fallen zulassig. 


iL 
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Exzentrizitat ¢ des Zapfenmittel punktes und den Winkel y (Abb. 2), ist namlich eine Funktion 
von G; wir werden sehen, daB die Gleichgewichtslage die GréBe von b und c entscheidend 


beeinfluBt. 


2, Rechnerische Ermittlung der Diampfungseigenschaften. a) Grun dlagen des 
Verfahrens. Die Grundlagen des Verfahrens mégen am Fall des Plattenstreifens der 
Breite 2a erlautert werden, der sich senkrecht zur Zeichenebene nach beiden Seiten ins Unend- 
liche erstreckt und auf einer Schmiermittelschicht gelagert ist, deren Trager eine unendlich 
groBe Ebene ist (Abb. 3). Plattenstreifen und Ebene sind einander parallel und werden relativ 
zueinander lings ihrer Normalen mit endlicher Geschwindigkeit bewegt. Wird dabei die 
Spalthéhe zwischen Streifen und Ebene verringert, so wird Fliissigkeit seitlich herausgedriickt. 
Dadurch entsteht in dem Spalt ein lings der Plattenbreite verdnderlicher Uberdruck gegeniiber 
dem Druck p, in der Schmiermittelschicht auSerhalb des Plattenstreifens. Bei umgekehrter 
Bewegungsrichtung, stroémt Fliissigkeit in den Spalt ein, und es entsteht ein Sog, der der Ver- 
gréBerung der Spalthéhe entgegenwirkt. Somit wird der 
Energiequelle, welche die Bewegung des Plattenstreifens 
und der Ebene relativ zueinander hervorruft, wahrend 
der Bewegung dauernd Energie entzogen. Wir wollen fiir 
diese zweidimensionale Strémung den der Geschwindigkeit 
der Relativbewegung ent- 
gegenwirkenden Uber- bzw. 
Unterdruck im Spalt er- 
mitteln. 

Die Navier-Stokesschen 
Gleichungen der zahen 
Flissigkeit, mit denen wir 
den zuvor erérterten Vorgang mathematisch beschreiben kénnen, reduzieren sich auf folgende 
Form, da die Tragheitsglieder und die itbrigen Reibungsglieder gegeniiber dem jeweils an- 
gegebenen Reibungsglied zuriicktreten: 


Ort aller Glefah- 
gewlehtslagen 


abhinngig von 6 Mittepumkt 


des Lagers 


Abb, 2. Abb. 3. 


Op O2u 
ao aes (3a) 
dp Ow 
a Oe (3b) 


In jeder der beiden Gleichungen verhalt sich das gréBte der vernachlassigten Glieder zu dem 
mitgenommenen Reibungsglied gréBenordnungsmaBig wie 


Darin bedeuten h die zeitlich verinderliche Spalthéhe und y die kinematische Zahigkeit des 
Schmiermittels. Abschatzungen zeigen, da® Re in praktischen Fallen etwa 0,1 betragen kann. 

Zu den Gleichungen (3) fiir die schleichende Strémung tritt noch die Kontinuitatsgleichung 
fiir die inkompressible Fliissigkeit 


Ou Ow 
PERI 


sowie die Randbedingung, daB au®erhalb des Spaltes itberall der Druck Po Vorherrscht }. 
Der Ansatz 


erfiillt die Kontinuitatsgleichung. Einsetzen in (3) zeigt, da® der Uberdruck P — Py von folgen- 
der Form sein muB: ; ze 


P — Po =p! [x* + F(z)] . 


* Wegen des unbedeutenden Einflusses der Rander x = -- a auf das Spaltinnere werden wir uns damit 
begniigen, die Randbedingung nur fiir die vier Kanten « = + a,z = 0 undz = h zuerfiillen 
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Nach kurzer Rechnung erhalten wir 
F(z) = — 22+ C,2+ C,, 
fae kf + Kye K;; 
Zwischen C, und K, besteht wegen (3b) die Beziehung 
(t 


Die Haftbedingung uv = 0 und w = 0 fiir z = 0 ergibt sofort K, = 0 und K, = 0, wahrend die 
Haftbedingung u = 0 fiir z = h (t) liefert 


K,=—* v(t) und somit C, = A(t). 
Mit Hilfe der Haftbedingung w = dh/dt fiir z = h(t) erhalten wir schlieBlich 
p(t) = 12,8 
Das Geschwindigkeitsfeld wird daher durch folgende Geschwindigkeitskoordinaten beschrieben: 
SER we -folah aay. 


Der Druckunterschied im Spalt gegeniiber den Verhaltnissen auBerhalb des Spaltes ergibt 
sich dann wie folgt: 


2a? +hz—2s dh 
Pe) Ob ee (4) 


Wir erkennen, daB er der Geschwindigkeit dh/dt entgegengesetzt gerichtet ist und daher der 
Energie quelle, welche die Bewegung erzeugt, Energie entzieht. Der Druckunterschied ist ferner 
der Geschwindigkeit proportional. 
Fiir unsere weiteren Untersuchungen ist nun wesentlich, daB bis auf eine unbedeutende 
Randzone in der Umgebung von x = + a mit guter Naherung ! 
x — a? dh 


Die Po OU ae 


: (4a) 


gesetzt werden darf. Nun ist aber (4a) gerade der Ausdruck fiir den Druckunterschied, der sich 
bei eindimensionaler Behandlung des Problems = 0 ergeben hatte. Wir kénnen uns daher 
im Folgenden bei der Behandlung ahnlicher ebener Probleme auf die eindimensionale Be- 
trachtungsweise beschranken, die die Geschwindigkeitskomponente w ganzlich vernachlassigt 
und unter Beriicksichtigung der Haftbedingung an Stelle der Geschwindigkeitskomponente u 
durch Mittelbildung in z-Richtung die an jeder Stelle x in x-Richtung strémende Fliissigkeits- 
menge q, einfithrt. Damit erhalt man bekanntlich aus (3a) 


d 12 
= ae (9) 
Aus Kontinuitatsgriinden mu8 dann gq, = — x dh/dt sein. Damit ergibt sich als Lésung fiir den 


Druckunterschied die Forme! (4a), so da®B der Nachweis fiir die Berechtigung eindimensionaler 
Betrachtungsweise erbracht ist. 


b)Dampfungseigenschaften des ebenen Gleitschuhes in ein- 
dimensionaler Betrachtungsweise. Der ebene Plattenstreifen sei nun gegen 
die Ebene um den Winkel f geneigt und fiihre eine translatorische Bewegung mit endlicher Ge- 
schwindigkeit lings der Normalen zur Ebene, also in 2-Richtung, aus. Dabei bewege sich die 
Ebene in x-Richtung mit der konstanten Geschwindigkeit U (Abb. 4). Unter der Voraussetzung 
schleichender Strémung erhalten wir bei eindimensionaler Betrachtungsweise folgende Glei- 


1 Die Formel (4a) gilt sogar (vgl. FuBnote 1 von S.148) genau fiir die Ebenen z = 0 und z = h(t), auf 
die es uns ja bei der Ermittlung der Dampfungseigenschaften des Schmiermittelfilms ankommt. 


ike 
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chung fiir den Druck p im schragen Spalt: 


0 U Ee 


Ox 


Darin ist hfolgende Funktion von x und t, wenn wir im Hinblick auf die Ermittlung der Damp- 
fungskonstanten gleich annehmen, daB der Plattenstreifen gegeniiber der Ebene x eine periodi- 
sche Bewegung ausfithrt: 


h(x, t) = (a) — x) B + esinwt = [a(t) — x|)B mit a(t) =a +5 sinwt. (7) 


Fir die Bestimmung von q, auf Grund der Kontinuitatsbedingung ist die Kenntnis der Stelle & 


der ,, Wasserscheide* erforderlich, fiir die dP _Qwird. Dort ist dann entsprechend (6) und (7) 


dx 
U-h(&, t) U [a(t) — 6) 6 
Yl &, t) = 5} —— 5} >; 
Damit 1aBt sich q, wie folgt schreiben: 
h(é, : 
get) SAE) 4 he — 2). (8) 
Nach kurzer Rechnung erhalten wir als Bestimmungsgleichungen fiir h(é, t) baw. fiir € 
6 KO; hil, 2) 
het) = AO, 1) (lt) (9} 
bzw 
_ _a(t) 1 
areas 


Durch Einsetzen dieser Werte in (8) 14Bt sich (6) mit 
folgendem Ergebnis integrieren: 


Me U/2 —a [h(0, t)—h(x,t) h(E, t) h?(0, t) — h?(«x, t) 
P—Po= 12h B h(0,t)h(x,t) 2 h?(0, t) h?(x, t) 


Dieser Ausdruck nimmt nach einigen weiteren Umformungen unter Benutzung von (7) und (9) 
schlieBlich folgende Form an: 
Cel G2 =e) alee) 
ie ran ara Varna Sy © (10) 
Bei festgehaltenem A(0, t) geht fiir a—> oo diese Forme] in die fiir den Druckunterschied bei 
parallelem Spalt iiber. 


Durch Integration von (10) itber die Spaltbreite erhalten wir folgende Gleichung fiir die 
Resultierende/Langeneinheit des Druckunterschiedes im Spalt: 


__ 12 w (U/2 —a) a(t) 
D) = Fire Gan| 2 ee) — Mage 2ah. (1) 


Das Integral 
2 n/w 


E=—|[Dhd 
0 


liefert nun den Verlust an mechanischer Energie je Langeneinheit des Plattenstreifens, der bei 
einer vollen Schwingung entsteht. Unter der Voraussetzung kleiner Schwingungen ergibt sich 


angenahert 
= 12 u a 21 
i De e BB ln(;, —— ) Bs q — | wae. (12) 


Der bei dieser Integration gemachte Fehler ]aBt sich durch Mitnahme des nachst hdheren Glie- 
des der bei der Rechnung durchgefithrten Reihenentwicklung abschatzen. Dieses betragt 


12ufl 1 1 21 \/e 2 5 
Paleo: a wa | (8) OTe 


| 


| 
| 


| 
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und kann unter der Voraussetzung kleiner Schwingungen gegeniiber (12) als vernachlassighbar 
klein angesehen werden, weil dann e/B <hy/Pp = apist. 

Durch Reihenentwicklung von (12) kénnen wir noch folgende Form fiir die durch Dampfung 
verlorene Energie erhalten: 


pk 1\2 1\3 1\4 : 
E =m |1+ 150— + 1,65(-) ae 1,63 (--) a 1,42(7) +o Jame, (12a) 


aus der wir durch den Grenziibergang a,—> oo, B > 0 und a, f = hy, = konst. den entsprechen- 
den Ausdruck fiir den parallelen Spalt gewinnen kénnen: 


3 
Beat om et. (12b) 
0 


Unter der Voraussetzung kleiner Schwingungen kann also die Dampfung entsprechend (12) 
bzw. (12b) als geschwindigkeitsproportional betrachtet werden. Indem wir die Dampfungs- 
konstante mit b bezeichnen, ergibt sich der Energieverl ust iiber eine volle Periode zu E = bwze?. 

Durch Vergleich mit (12) erhalten wir also als Dampfungskonstante/Langeneinheit fiir den 
ebenen Gleitschuh 


_12p ay 21 
eer R tn) Fa Ce 
bzw. bei parallelem Spalt 
b als 
= Ee (13b) 


Wir erkennen, daB die Geschwindigkeit U nicht in die Dampfungskonstante (13) eingeht. Ihr 
Kinflu8 ist nur insofern wesentlich, als bei gegebener Belastung des Gleitschuhes die die Spalt- 
abmessungen bestimmenden Gréfen a, und f von U abhangen. Bei der nun folgenden Berech- 
nung der Dampfungskonstanten fiir den im Gleitlager schwingenden Zapfen kénnen wir also 
von vornherein seine Rotation unberiicksichtigt lassen; sie spielt dann lediglich bei der 
Bestimmung der Gleichgewichtslage des Zapfenmittelpunktes im Lager eine Rolle. 

c) Dampfungseigenschaften des im Gleit- 
Pacer rotierenden Zapfens (ebenes Problem). 
Wir untersuchen zunachst den Fall, daB8 der Zapfen entsprechend 
Abb. 5 in der Lagerschale gelagert ist. Wenn wir mit 6 das radiale 
Lagerspiel, mit ¢ die Exzentrizitat der Zapfenachse und mit e wieder- 
um die Amplitude der Zapfenschwingung bezeichnen, ergibt sich die 
Spalthéhe als Funktion des Zentriwinkels ~ und der Zeit t zu 


h = h(g, t) = 6 — (e — esin wt) cosy 
und die Geschwindigkeit der Spalthéhenanderung zu 


h=ewcosa@tcos@. 


Analog (5) gilt hier 


dp 12 uw 
gee Serie 
Mit q, =— rewcoswtsin y (Kontinuitatsgleichung) folgt 
dp sin p 
= 0S Wt. 
De 12 u te ee cag yeep 


Indem der Druck auBerhalb des Spaltes der Einfachheit halber gleich Null gesetzt wird, ergibt 
sich folgender Druck im Spalt: 


” 6 wre ; ] I 
es é€—esinat | [6 — (e — e sin o t) cosy? [6 —(e —esinw t) cos 9}? | ; 


Die Resultierende dieses Druckes wird daraus wie folgt gefunden: 


” 
D = 2 f{ p(y) cos pr dy. 
0 
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Unter der Voraussetzung kleiner Schwingungen (e < 6 — ¢) 1aBt sich dann auf dem gleichen 
Wege wie in b) folgende Dampfungskonstante ermitteln: 


Yeoh sii (cosy—y)siny | (14) 
b = 12 u(~ = arc tg (ptt tg (1 —x cos wp)? (1 — x) | 


e E 
mit We es 


Beim rotierenden Zapfen ist nun die Verbindungslinie zwischen Lager- und Zapfenmittel- 
punkt nicht mehr Symmetrielinie der Lagerschale ; somit ist die engste Spaltstelle keine ,, Was- 
serscheide“. Die Folge ist, daB der Versuch, fiir eine Lagerschale, z. B. fiir das 120°-Lager, bei 
unsymmetrischer Zapfenlage die Dampfungskonstante zu berechnen, einen erheblichen Rechen- 
aufwand mit sich bringen wiirde. Auf Grund der folgenden einfachen Uberlegung diirfte je- 
doch dieser Aufwand, zumindest fiir die praktische Anwendung der vorliegenden Theorie, un- 
noétig sein. 

Bekanntlich ist an der Tragwirkung des Lagers nur ein kleiner Bereich um die engste Spalt- 
stelle maBgeblich beteiligt ; jedem Praktiker ist jedoch der groBe HinfluB gelaufig, den auch das 
Spiel der oberen Lagerschale auf die Laufruhe einer rotierenden Maschinenwelle hat. Es er- 
scheint daher sinnvoll fiir die Berechnung der Dampfungseigenschaften des Lagers, dieses als 
vollumschlossenes Halslager (Abb. 6) zu betrachten, weil dadurch angendhert auch die Wir- 
kung der oberen Lagerschale mitberiicksichtigt wird, die besonders fiir groBe Werte der Lager- 
kennzahl] G — also gerade fiir die Falle gefaihrdeter Stabilitat — an der Dampfung erheblichen 
Anteil hat. Unter dieser Voraussetzung erhalt man dann aus (14) fiir y = a die Dampfungs- 
konstante fiir Zapfenschwingung in der 2’-Richtung wie folgt: 


r \3 I 
= 1245) ae (15) 


Fir die Dampfungskonstante in der «’- Richtung 
wird angenommen, das diese praktisch unab- 
hangig von der Exzentrizitat ¢ des Lagerzapfens 
ist und daher durch Nullsetzen von y aus (15) 
gewonnen werden kann: 


b= in as, : (15a) 


Abb.6. Abb.7. 


d) Diampfung der ebenen Platte und des Zapfens endlicher 
Lange. Bei den bisher behandelten Aufgaben war stets vorausgesetzt, daB in y-Richtung 
(Abb. 7) keine Strémung stattfindet. Beim wirklichen Lager tritt aber ein AusfluB des Schmier- 
mittels nach beiden Lagerenden auf, der die dimpfende Wirkung herabsetzt. Da das Problem 
fiir den Zapfen endlicher Lange bisher nicht gelést wurde, soll hier nur eine Abschatzung ge- 
wonnen werden, indem der EinfluB des Endausflusses beim Zapfen der entsprechenden Lé- 
sung dieser Aufgabe bei der ebenen Platte entnommen wird. 

Das Problem findet seine mathematische Formulierung unter den zulassigen Voraussetzungen 
der schleichenden Strémung in folgenden drei Gleichungen: 


coe = a Op PA 0qx Oqy : 
pe fe ay eas aga ce oe 
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Wegen 0q./0y = 0q,/0x lassen sich diese Gleichungen zu folgender inhomogenen Potential- 
gleichung vereinen: 


Op aC pee 12h 
Ox? = ayie BB (mo, 


Das Einsetzen des Liésungsansatzes 


Gof (2i +1) 2% 


Gases) 2 
P=G,> ara ps (28+ YE = + C, (a? — a7) + C, 
as eae Cof(2i+ 1) >— 
ergibt 
_— 6 ee 


Die Erfillung der Randbedingung p = p, fiir x = -+ aerfordert C,; = py. Im Hinblick auf die 


Summenforme] 
1 (=D 9: wx ff a\2 mu x\2 
peers Lg ae (a) tal 


i=0 


folgt wegen der Randbedingung p = p, fiir y = +b 


oo : Gof (27 + 1) = 
me | L92) a2 (— 1): ; 2 2\| 7 
P—Po=|—aae pages QE NE ae Hof (eel, (17) 
at Gof QE +I 


Gleichung (17) geht fiir b > oo in die entsprechende Formel fiir den unendlich langen Platten- 
streifen (4a) iiber. 
Durch Integration von (17) iiber die Plattenflache 


ba 
D=4] | (p— py) dx dy 
00 


ergibt sich mit 1 = 2a und L = 2b folgender Ausdruck fir den geschwindigkeitsproportionalen 
Widerstand der Platte gegen Auslenkung in z-Richtung !: 


fb, 
: Wa (Oa 


D 
Bola Lay - (21+ 1) 


Fir den Widerstand je Langeneinheit des unendlich langen Plattenstreifens wiirde sich ergeben 


haben 
ee (19) 


eg ae 
Vergleichen wir hiermit die entsprechende Forme! fiir den unendlich langen Zapfen 


OTS ge see lZura; 
SEO Ha IA NO 


so lassen sich diese in (19) iiberfihren, wenn wir 
jl2ar=1 und h=s~l—y bew. Ylar=l und h=d 


setzen. Indem wir von diesen Beziehungen zwischen unendlich langem Zapfen und unendlich 
langer Platte fiir die naiherungsweise Ermittlung der Dampfungskonstanten: des endlich 


1 Wenn in (18) an Stelle von — 3 wu h/h® das Produkt G @ (elastischer Schubmodul G- Drillwinkel #) 
gesetzt wird, entspricht die rechte Seite von (18) dem in einem Drillstab mit Rechteckquerschnitt | + L 
durch Schubspannungen iibertragenen Drillmoment. Es wird also dem bekannten Strémungsgleichnis 
fiir die Stabtorsion ein weiteres hinzugefiigt. 
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langen Zapfen Gebrauch machen, erhalten wir mit Hilfe von (18) folgende Ausdriicke fiir diese 
Konstanten in den 2’- und x’-Richtungen: 


Py Toei ty 7 | 
12uaLr' 1923 ro 2V12 a7. 
arg = RI pa l—— laa es (2i+1)5 ? (20) 


bee = (1— 22) be 
Aus diesen ergeben sich die entsprechenden Dampfungskonstanten fir die x- und z-Richtungen 
wie folgt: 
bez = 5,,, sin? y + b, _ cos? y, | 
b, — b, 2! cos? y-+ by x sin? VY > ( (21) 
b,, = bz-= cos y sin y by, [1 — (1 — x*)??]. | 
Die Auswertung der Ausdriicke (21) ergibt, daB die DampfungsgréBen b,, und b,, praktisch 
identische Funktionen der Lagerkennzahl G sind, und daB ihnen gegeniiber b,, vernachlassig- 
bar klein ist. Die Dampfungseigenschaften werden daher mit ausreichender Genauigkeit durch 


nur eine DampfungsgréBe gekennzeichnet, die wir nun ohne Indices mit b bezeichnen wollen. Thr 
Verlauf ist in dimensionsloser Darstellung fiir einige Verhaltnisse d/L der Abb. 8 zu entnehmen. 


wachsendes 
ae 


x \ ak 
cas, 
0 50 100 150 200 


Abb. 8. Abb. 9. 


3. Ermittlung der Federzahlen fiir das 120°-Lager. Dieses besonders in der amerikanischen 
Praxis haufig verwendete Lager hat S. J. Needs! in einer ausfithrlichen experimentellen Arbeit 
untersucht. Seine Untersuchungsergebnisse hat A. C. Hagg fiir die Ermittlung der Federzahlen 
mit Hilfe des Verfahrens von Stodola ausgewertet. Der Grundgedanke dieses Verfahrens sei im 
folgenden kurz erlautert. Die Punkte A und B der Abb. 9 stellen zwei benachbarte Gleich- 
gewichtslagen der Zapfenmitte dar. Die entsprechenden vertikalen Lagertragkrafte sind P 
und Q. Wahrend des durch die Lagerkennzahl G gekennzeichneten Betriebes befinde sich der 
Zapfenmittelpunkt in A. Durch irgendeine Stérung werde er aus seiner Gleichgewichtslage A 
nach B’ ausgelenkt. Die dann von dem Schmiermittelfilm auf den Zapfen ausgeiibte Kraft Q’ 
soll nun dem Betrage nach gleich der Kraft Q sein, weil B und B’ die gleiche Exzentrizitat ¢ 
besitzen. Fiir die Bestimmung der Richtung der Lagerkraft in B’ wird angenommen, daB der 
Winkel zwischen Q und dem Fahrstrahl BO gleich dem zwischen Q’ und B’O ist. 

Unter der Voraussetzung kleiner Auslenkungen ergeben sich damit folgende Federkon- 
stanten fiir die x’- und z’-Richtungen: 


P ONE 
cosy, BM tet lh SM ie aye 
(22) 
A oP aes 
sin y , Cys’ = FZ cosy + —siny tg yp. 


1S. J. Needs, Trans. A.S.M.E., 56 (1934), APM — 56 — 16, S. 721. 
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Die entsprechenden Federkonstanten fir die x- und z-Richtungen sind 


Pre P 
Co =~ (sin y tg p + cos 8 es —~ (sin y — cos y tg Y) » 
(23) 


WeOP Ss), oP 
Gx => siny , Cz = 5 COSY. 
Zur numerischen Auswertung dieser Formeln mige die von Hagg auf Grund der Ergebnisse von 
Needs gefundene Beziehung zwischen der Lagerkennzahl G und der bezogenen Exzentrizitat 
benutzt werden: 
Caner x a= 1, 


Damit lassen sich P und 0 P/d¢ in (22) und (23) schreiben: 


Tae oP B 
— mya B Se 
Pa2 5 r(5) enacts und at eae 
Darin sind a und f dimensionslose Konstante, die lediglich von dem Verhaltnis 2 r/L = d/L 
abhangen. Zur Erleichterung der praktischen Anwendung sind die Federzahlen (23) als Funk- 
tion der Lagerkennzahl Gin Abb. 10 bis 13 dimensionslos dargestellt. 


Cre in 
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| et 
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1 
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Abb. 12. Abb. 13. 


4, Zusammenfassung. Die vorliegende Theorie erméglicht die rechnerische Ermittlung der 
hydrodynamischen Dampfungs- und Federeigenschaften von Gleitlagern. Die Ergebnisse 
dieser Theorie sind in den Abb. 8 und 10 bis 13 dimensionslos dargestellt. Diese Diagramme 
erméglichen die numerische Bestimmung des Schwingungsverhaltens von rotierenden Wellen 
unter Beriicksichtigung der Gleitlagereigenschaften, die besonders im Hinblick auf die Lauf- 
stabilitat von wesentlicher Bedeutung sind!. Auch die kritischen Drehzahlen und das Ver- 
halten der Wellen in der Umgebung einer solchen Drehzahl erfahren durch die Gleitlagereigen- 
schaften eine nennenswerte Beeinflussung. So zeigte z. B. die fiir einen Turbinengenerator- 
laufer durchgefithrte Schwingungsrechnung neben der zu erwartenden Herabsetzung und Auf- 
spaltung der Wellenkritischen eine betrachtliche Dampfungswirkung der Gleitlager, welche die 
Amplitude der zirkular-polarisierten Schwingung in der ersten kritischen Drehzahl auf das 
etwa sechsfache der Laufer-Exzentrizitat beschrankte. 


(Eingegangen am 21. Juli 1953). 
Anschrift des Verfassers: apl. Prof. Dr.-Ing. E. Pestel, Hildesheim, Weinberg 60. 


1 E. Pestel, Drehzahlfremde selbsterregte Schwingungen von Maschinenwellen, Verlag der Allianz- 
Maschinenversicherung, Wiesbaden 1953, 
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Die Wurzeln der charakteristischen Gleichung 
fiir das kreiszylindrische Rohr. 


Von E. R. Berger. 


Das System der Differentialgleichungen fiir die Verschiebungskomponenten eines biege- 
steifen Rohres wurde erstmalig von Miesel 1 aufgestellt. Spiater hat Fliigge *, von anderen Ver- 
schiebungsgréBen ausgehend, das entsprechende Gleichungssystem unter Beriicksichtigung — 
der Querdehnung aufgestellt. Es ist linear, mit konstanten Koeffizienten und, falls nur Rand- 
krafte angreifen, homogen. Man lést es durch den Ansatz 


exp (— Ax/a) - exp (im@) , 


wobei a der Rohrhalbmesser ist, x die Koordinate in Richtung der Erzeugenden, g das Azimut. 
Durch Einsetzen in die Fliiggeschen Gleichungen ergibt sich die charakteristische Gleichung 
[Fliigge, S. 125, Gl. (78)]: 


A8 — 4(m? — v/2) A8 4 ere + 6 m4 6m) | a) 
= 9 it? it (Fy) md No Gey 


Darin bedeutet » die Querdehnungszahl, k den Schalenparameter, der immer viel kleiner ist 
als Eins. Fir ein geschlossenes Rohr mu8B m ganzzahlig sein, A ist aus (1) zu bestimmen. 

Die charakteristische Gleichung von Miesel [S. 35, Gl. (27)] stimmt damit genau iberein, 
auBer daB neben 1/k noch der Summand 3 steht, ein Unterschied, der véllig bedeutungslos ist. 

Ausgehend von einem vereinfachten Ansatz hat Aas-Jakobsen® eine ahnliche Gleichung 
erhalten [S. 402, G1. (144)], die sich von (1) nur unwesentlich unterscheidet: anstelle der in (1) 
unterstrichenen Koeffizienten 6, 4, 2 erhalt. Aas-Jakobsen 4, 3, 1, ferner ist 1/k ersetzt durch 
1 + 1/k; beide Unterschiede sind fiir das Folgende ohne Belang. Auferdem fehlt bei Aas- 
Jakobsen die Querdehnung. 

Gleichung (1) ist als biquadratische Gleichung in A? mithsam aufzulésen. Miesel verwendet 
dazu die kubische Resolvente, und gibt daneben fiir kleine Werte von m Naherungsformeln fiir 


Aan [S. 37, Gl. (35)], die aber bis zu 5% Fehler ergeben, sobald 


1h =e 


Wire iiss 
m?(m 1) ar 


wird, also nur fiir ganz kleine Werte von m brauchbar sind. Fliigge [S. 126] gibt ebenfalls ein 
Naherungsverfahren fiir kleine m: er streicht einmal die hohen, einmal die niedrigen Potenzen 
von A und erhalt so zwei quadratische Gleichungen fiir 42. Das Verfahren versagt, sobald A in 
die Gré8enordnung von m kommt (ein Kriterium hierfiir ist nicht angegeben), auch erfordert 
es noch das unbequeme Ausziehen der Quadratwurzel aus dem komplexen Wert 42. Aas- 
Jakobsen |S.402—3] list die Gleichung fiir beliebige Werte m durch ein Iterationsverfahren, 
das aber immerhin noch einige Rechenschritte erfordert. 


Die folgende Methode gibt eine einfache Naherungsformel fiir 42, wobei nur vorausgesetzt 
wird, daB Vk < 1 ist, wahrend m beliebige Werte annehmen kann; der Fehler ist von der 
GréBenordnung Vk. Ist auBerdem noch 


Lo eeeyas 
salle 


so kann man Realteil und Imaginarteil von A nach Potenzen von mentwickeln und erhalt mit 
fiinf Gliedern eine Genauigkeit mit einem Fehler von unter 1%. 


1 K. Miesel, Ing.-Arch. 1 (1929—30), S. 22. 
2 Fliigge, Statik und Dynamik der Schalen. Berlin 1934. 
8 Aas-Jakobsen, Bauing. 20 (1930), S. 394. 
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Man setze 


(2) 
ym=n; yR=n, | 


dann wird (1) zu 


xt — 4 x8 (n—v y/2) 4 x?(1-+ 6 n?—6ny)—4«n(n—y)(n—y+ yy/2) + n2(n— y)?=0. (3) 


Setzt man weiter 


x= y|n(n— y), (4) 


so ergibt sich aus (3) nach einigen Umformungen 


if n — y/2 1 l—y» 1 1 
Perigo 4 i ae eis v= ate aie | Se Giese (Ve 5 
ee ae eS de eee Y 
Dies ist eine reziproke Gleichung fiir y, wenn man von dem Glied y — 1/y absieht; es ]aBt sich 


aber zeigen, das dieses Glied im Vergleich zu den iibrigen Gliedern nur von der Gré®enordnung 
y ist und daher vernachlassigt werden kann. 


SchlieBlich setzen wir noch 


1 z 
Nea = 6 
oY ain —y) ms 
und erhalten damit aus (5) 
a A ees Nh gc eat i eee (7) 
Hier kann y? = k neben 1 vernachlassigt werden; bezeichnet man noch 
n—y/2=n, also nm = (m*— 1/2)y, (8) 
so wird 
Z=itt2n. 
Nach (6) ist 
opened L1=0, 
20 foe 
oder (4) eingesetzt, 
x*—zx+tn(n— y)=0. (9) 


Daraus ergibt sich 


= =Slit2% +iyl —4in— 7). 
Wieder wird y? neben | vernachlassigt und wir erhalten 
x =—(1—2in+y1—4in). (10) 


Dann ist aber, wie man leicht bestatigt, 


je—E (1 + yi Ain). nb) 
Damit ergibt (2) fiir die gréBere Wurzel 
Vx, 
eos 
1 Vi 
= eet Ly Ain) (12a) 
2 Vk 
ee fie Van ai). (12) 
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Falls n klein ist, wird bei der Berechnung der kleineren Wurzel A, die Differenz in der 
Klammer von (11) sehr ungenau. Besser schlieBt man aus (9) 


te n(n —y) 


xX, 
2 
*y 


und daher fiir den Wert im ersten Quadranten 


is a m Vm? — 1 j (13) 


A, 


Solange n < 1/4 ist, kann man die Quadratwurzel in (12a) nach Potenzen von n entwickeln 
und erhalt damit 


peo! 
\2 Vk 
As — Tmt {((1—2n2+ 14n4) + (n— 573)]+ i[(1—2n?+ 14n4)—(n— 5n’)]}. (14b) 


Da die Teilreihen alternieren, ist der Fehler kleiner als das erste vernachlassigte Glied, also 


{((. + — 5 m4) + (w—29%)] + i[(1 + W—5H4)— (n— 27), (14a) 


feo { Realteil: negativ, <14n°(1+ 3n)< 25%/,, 
“' 1" | Tmaginarteil: negativ, <1475(1— 3%) < 3,5%y, 
. { Realteil: positiv, <42n*>(1— 3n)< 11/,, 

fir Ass 


| Imaginarteil: positiv, <4275(1+ 3 n)< 75%. 


Die oberen Schranken gelten fir n = 1/4, sie werden auch dort nicht erreicht. Zu diesen 


Fehlern kommt noch ein Fehler von der GréBenordnung Vk dazu, der von der Vernachlassigung 
in (5) herrithrt. 


Man kann ihn z. T. korrigieren, indem man den Wert A nach (12) oder (14a) mit 
Let ye (15) 
multipliziert, bzw. 4 nach (13) oder (14b) mit dem konjugierten Wert. 
Zahlenbeispiele: Fliigge [S. 126] berechnet die Wurzeln fiir 
b= 3,105 ~10-% “und” “mr 8; 
vist nicht angegeben, wir setzen es gleich Null. Damit wird 
y= Yk =1,925-10-3, n = 63,5-y = 0,1222, 
1/y2 y = 16,116, Jy/2-8- ¥63 = 1,9698. 
Die Reihenentwicklungen (14) ergeben 


1— 1,— 
m= 0,0149 i= 01299. |) 2272 = 0.0290 nm = 0,1222 
--57nt=-=0,0011 | —27? =—0,0036 1474— 0,0031 | —5n3 ——0,0091 
1,0138 0,1186 0,9732 0,1131 
Bie 324 _ {1,0863 
1,0138 + 0,1186 = | 0.8952 0,9732 + 0,1131 = | 0.8601 


also 
A, = 18,254 14,437, A, = 2,140 + 1,604 2 
(Fligge gibt an 18,22 + 14,427 bzw. 2,14 + 1,69.) 
Nimmt man dasselbe k, aber m = 40, so wird n = 3,079; hier sind die Reihen nicht mehr 
verwendbar, wir miissen nach Formel (12b) arbeiten. Man erhalt 
A, = 48,09 + 9.69%, 4, = 31,98 + 6,43 i 
(Fliigge gibt an 48,17 + 9,747 bzw. 31,89 + 6,51 1.) 


Aas- Jakobsen [S. 403, Zahlentafel VI] behandelt dasselbe Beispiel, jedoch ist dort durch 
einen Druckfehler k= 3,705 - 10-° angegeben; bei seinem ersten Beispiel (k= 2 -10-* m = 2) 
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ist A, falsch berechnet: in der Formel wurde irrtiimlich B anstatt /B eingesetzt. Bei richtiger 
Rechnung ergibt sich A, = 0,0926 + 0,0917 i. 


Miesel [S. 37, Tab. 1] gibt u. a. fiir k = 0,001 und m = 2 als strengen Wert an 
A, = 4,50 + 3,517, A, = 0,419 + 0,3547 
und nach seiner Naherung 
A, = 4,50 + 3,517, A, = 0,465 + 0,386 i. 
Unsere Formel (14) ergibt 
Ay = 4,45 + 3,597, A, = 0,471 + 0,381 i, 
und korrigiert nach (15): 


A, = 4,514 3,52i A, = 0,465 + 0,388 i. 
Rechnet man mit den Werten von Miesel die Kontrolle 
|A2| - |A3| in (a — |) 7127, 


so ergeben die genauen Werte 9,80, die Naherung 11,90; die Naherung diirfte also genauer sein, 
als die durch einen umstandlichen Rechnungsgang erhaltenen ,,genauen‘* Werte. 


(EHingegangen am 23. Juli 1953.) 
Anschrift des Verfassers: Dr.-Ing. E. R. Berger, Wien IV, Technische Hochschule. 
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Ein Beitrag zum Randstérungsproblem isotroper Kreiszylinderschalen. 
Von D. Riidiger. 
1. Einleitung. Im XX. Band dieser Zeitschrift hat Herr W. Zerna' fiir die Berechnung 


der Randstérungen kreiszylindrischer Tonnenschalen eine komplexe Differentialgleichung 
vierter Ordnung fiir eine Spannungsfunktion angegeben und mit deren Lésung den Spannungs- 
und Verschiebungszustand fiir Kreiszylinderschalen zwischen zwei Endbinderscheiben be- 
schrieben. Mit dem vorliegenden Beitrag werden ausgehend von den Beziehungen (13) und (14) 
der zitierten Arbeit die zur Lésung der Randaufgabe notwendigen Schnittkraft- und Verschie- 
bungsgleichungen in einer fiir die Berechnung zweckmaBigen Form so zusammengestellt, dah 
eine numerische Auswertung von Rekursionsformeln fortfallt. 


2. Geometrische Grundlagen, Definitionen. Diese sind in Abb. 1 und 2 zusammengestellt. 
Die Bezeichnung der Schnittgréfen ist so gewahlt, daB der erste Index den Ort, der zweite 
Index die Richtung des betreffenden Vektors angibt. 


z aS Ke 


‘ az @ i 
‘ ) fe No 
“eZ 


Abb. 1. Geometrie der Schalen- 
mittelflache und Koordinatensystem. Abb. 2. Definition der Schnittkrafte und Verschiebungen. 


3. Differentialgleichung. Aus den simultanen Differentialgleichungen (13) und (14) der 
Zernaschen Arbeit ergibt sich die Differentialgleichung des Problems als Determinante der 


0 per A Se eared SE eae 
peratoren = 9H T 38 un Sm 
AAAA® +*=* or" = 0. (1) 


Hierin ist » die Querdehnungszahl, k die ittbliche Abkiirzung fiir h?/12a?, und Striche bedeuten 
Ableitungen nach é. 


4, Lésung der Differentialgleichung. Als Liésungsansatz wird wegen den in x = + L/2 


geforderten Randbcdingungen der Ansatz = c- elm+149 cosh E mit A=na a/L verwendet. 
Die Einfithrung des Lésungsansatzes in die Differentialgleichung (1) liefert die biquadratische 
Gleichung 
m4 + ea = (), 
Mk 


Werden die Abkiirzungen 
x, =VV0 +e? 
ig = V0 See) ne = ee ee 


ae 
ss Moen aioe. eae 

1 W. Zerna, Ing.-Archiv. 20 (1952) S. 357. 
* Die Existenz der Differentialgleichung (1) fiir Kreiszylinderschalen hat bereits im Jahre 1952 
H. Neuber mit Hilfe der allgemeinen Schalentheorie unter Vernachlassigung von Gliedern, die gegentiber 
ahnlich gebauten mit dem Faktor k behaftet sind, in einer nicht veréffentlichten Arbeit nachgewiesen. 
Ihr theoretischer Teil ist im wesentlichen in der wissenschaftlichen Zeitschrift der T. H. Dresden 


Wea worden. H. Neuber, Z. angew. Math. Mech. 29 (1949) S. 97; Wiss. Z. T. H. Dresden 2 


| 


— 
—_ 
pe 
io) 
— 
sS 
H 
—— 
= 
Contant 
fom! 
ice} 
— 
in) 
fe 
pa, 
— 
icp} 
— 


i) 
io) 


XXII. Band 1954 Riidiger: Ein Beitrag zum Randstérungsproblem isotroper Kreiszylinderschalen. 16] 


definiert, so kann man die Spannungsfunktion @ in der folgenden Weise anschreiben: 
P= cosh § [e—%7 (c, sin uw, P + c, cos My P) + e~*? (cysin oP + cy cos My) l (3) 
+ e~*® (es sin uy, @ + C4 cos yw) +e? (C7 Sin ly + Cy COs fly@)]. J 
Die Schnittkrafte und Verschiebungen werden aus der Lésung (3) durch Differentiation ge- 


wonnen. Dabei lassen sich nach einem Vorschlag von U. Finsterwalder1 die Ausdriicke fiir dic 
Spannungs- und VerschiebungsgréBen simtlich auf die Form bringen: 


W= Aha 3 {+ Cm [(A, a, — B, B,) sin 4p + (A, 6, +B, %) cos yp | + 
ree = [(A, 2 — By By) sin up + (Ay Py + By ay) cos ppp] + (4) 
+ e~%°1(G, 0, — Hy B,) sin wo + (6, 6, + Hy «) cos pty] +- 
+ e~%°[(6, O, — Hy By) sin wa + (GC, By + Hy Op) COS fly ]} « 
Die Faktoren y und 9; , f; (i = 1,2) sind fiir » = 0 in der nachfolgenden Tabelle 1 eingetragen. 


Tabelle 1. 

Ww | yp | Oy By bs fs 
Nop ae ll aie 1 0 a 1 5 
Qos | + cE + [oy py (1—e)] | + [oy (1 — ©) py] | + Lema(1 + €) x09] | +09 (1 +e) tpg] 
ONo x mn 

Paul ee ae iy + og oy 
Mox tooo == 1 — (1 + ¢) ote il aol (pes 
Dees ule ee + fy ty] + [i — My] + [ty — m9] ee Cees 

1 
ae — (1+e) ss 1 Sea Gta a i 
a Wee + [ey (1—e) ay] | Ly (Le) to] | Dea (+e) tite] | Lite +8) — x9] 
2a 
E-w sac 
0 + 1 0 = 1 
a «2 
Mz + 222 
ay 
E-0 (4 
: - My =a “al ie Me ar He 
Moy + ed ee é 
Qe z Book “i 1 y 1 -- 1 + 1 
Ovs : a6 2 ee) (2 46) ae 2 sel ete) 


Die Tabelle 1 enthalt auBer den bereits definierten Schnittkraften und Verschiebungen, die 
Winkeldrehung # der Kreistangente (positiv mit dem positivem Wirkungssinn des Biege- 
moments M,.), sowie zur Belastung der Randglieder die mit den Drillmomenten verkniipften 
a 0 Ove = d der verein- 
uerkrafte Q,, =Q,.+0M,,/0% und Q,,=Q,2 + 0M,,/a op. Aufgrun 
ie Naherrae aye SchuitteroBeu sind die Schubkrafte N,, = Nyy und die Drillmomente 
M,, = M,x gleich. Damit sind alle Schnittkraft- und Formanderungsgleichungen gegeben, 


ee cur Lésung der Randaufgabe notwendig sind. 


1 U. Finsterwalder, Ing.-Archiv 4 (1933) S. 43. 
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Die vorliegenden Gleichungen lassen sich auch aus der von A. Aas Jakobsen1 angegebenen 
ersten Naherung der Differentialgleichung fiir die Verschiebung w berechnen, wenn man, wie 
auch U, Finsterwalder? vorschlagt, von der Higenschaft gedimpfter Schwingungen Gebrauch 
macht, daB die niederen Differentialquotienten kleiner als die héheren sind und fiir eine 
Naherung vernachlassigt werden kénnen. 


5. Zusammenfassung. Ausgehend von den simultanen Differentialgleichungen (13) und 
(14) der Zernaschen Arbeit wird fiir die Randbelastungen isotroper Kreiszylinderschalen die 
von H. Neuber auf anderem Wege aufgestellte partielle Differentialgleichung achter Ordnung 
fiir eine Spannungsfunktion abgeleitet. Die Lésung des Spannungs- und Verschiebungs- 
zustandes fiir Zylinderschalen zwischen zwei Endbinderscheiben kann auf eine einheitliche, 
von U. Finsterwalder vorgeschlagene, Form gebracht werden. Die numerische Behandlung der 
Randaufgabe fiir die Vielzahl von Schalentypen wird dadurch erleichtert, daB die Werte x;, 4; 
und o;, 6;(t = 1,2) in einfacher Form von nur einer FormgréBe abhangen. Fiir die Be- 
schreibung des Zusammenhanges der Schale mit Randgliedern sei auf die Arbeiten von 
U. Finsterwalder und von F. Dischinger® hingewiesen. 


(Eingegangen am 3. August 1953.) 
Anschrift des Verfassers: Dr.-Ing. D. Riidiger, Dresden N 6, Kénigsbriickerstr. 70. 


1 A, Jakobsen, Bauingenieur 20 (1939), S. 394. 
2 Siehe Note S. 161. 
8 F, Dischinger, Beton u. Eisen 34 (1935), S. 257. 
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Zur Biegungstheorie einer Vierpilzplatte mit rechteckigen Stiitzflaichen. 
Von W. Miller. 


1. Kinleitung. In einer friiheren Mitteilung ! haben wir die Durchbiegung einer querkraft- 
frei eingespannten, bzw. durchlaufenden Vierpilzplatte betrachtet, die gleichmaBig belastet 
und von vier unendlich diinnen, symmetrisch zu den Mittellinien gelegenen Siulen gestiitzt 
wird. Da wir bei der dort verwendeten Fourierschen Methode bereits von der Voraussetzung 
rechteckiger Stiitzflachen ausgegangen sind, den Grenziibergang zu unendlich kleinen Flachen 
aber schon bei der Aufstellung der Druckfunktion vorgenommen haben, kénnen wir un- 
mittelbar von dem ersten Ansatz wieder ausgehen und die Durchbiegungsfunktion w durch 
eine Summe dreier Reihen, zweier einfacher und einer Doppelreihe darstellen. Es wird dann die 
nachste und nicht ganz einfach zu lésende Aufgabe sein, diese Reihen teilweise zu summieren, 
oder in einfache, gut konvergierende Reihen umzuformen, wobei sich dann ganz zwanglos die 
Unterscheidung der verschiedenen Felder der Plattenebene ergibt, in denen die w-Funktion 
verschiedene Formen annimmt. Bei der Umrechnung werden wir einige Reihenformeln be- 
nutzen, die sich aus den frither gegebenen Formeln durch Integration gewinnen lassen ” 


2. Aufstellung zweier grundlegender Reihenformeln. Bei der Behandlung der Vierpilz- 
platte mit kleinen Stiitzflachen ist folgende Reihenforme! benutzt worden: 


4 


oz Sofoy+2 Ginw oy Sinw ee + Cof a baer 
> Cos my 2 2 2 a 
n 


2 12)2 wt 4 (1) 
(an2 1 /@)>) 16 w? Gin? oe 2@ 


(Oy <2), 
wobei n die Folge der geraden Zahlen durchlauft. Mit wm = 0 erhalten wir daraus 


cosny 1 /xA ; 
ni ales el aes 2a — yi (2) 


n=2, 4,6 


cos n x wt [1 5 
> Se @ rei Cen 0 o'). (2a) 


Integrieren wir jetzt beide Seiten von (1) nach y von 0 bis y, so erhalten wir die neue Formel 


oder mit y = 70 


sinny cd >: o ; Blo etn ( 
nina = rocinoy 2Sinw> 2Sin oF 7 
n=2,4,6... Loo Cin’ © =. (3) 
* toy Cojo(2 Nese paca fl 
oad, 2 y (" 2@8\2 oa We 
Ebenso folgt aus (2) 
BG ee aay! oe et as 
23 n> a ZU ty RY as (4) 
n=2,4... 
oder 
sinnaw@  m°[@ Ue Leer eee ; 
eee = alst gil oe one eee Oe Cy 
n=2,4... 


3. Aufstellung der Reihen, Wir wahlen das Bezugssystem wieder so (Abb. 1), daB der 
Anfangspunkt in der linken unteren Ecke des Grundrechtecks liegt, dessen Seiten a und 6 sind 
und die Achsen mit zwei Rechtecksseiten zusammenfallen. In den Rechnungen werden wir 
statt x und y meist die dimensionslosen Gréfen x/a = E, y/a =n verwenden, und das Seiten- 
verhiltnis b/a soll mit 2 bezeichnet werden. Ferner mégen die rechteckigen Stiitzflachen der 


1 W. Miiller, Ing.-Archiv 22/1954, S. 60; vgl. auch Oesterr. Ing.-Archiv, 8 (1954) S. 1. 
2 Vel. auch W. Miiller, Oesterr. Ing.-Archiv, 6 (1953), S. 404. 
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Sdulen mit den Mittelpunktskoordinaten &, 3; 1— &. moi» &o> A—%31 = eae A — No die 
Seitenlangen 2a und 2 af haben und dem Grundrechteck parallel liegen. Die gleichmaBige 
Belastung der Platte je Langeneinheit sei py. Bei einem Stiitzdruck q der Saulen haben wir 
die Gleichgewichtsbedingung 

Poth =1l6qga%xB oder’ pyA—1Ogap. 
Mit Hilfe der mehrfach angewendeten Fourierschen Methode haben wir fiir die Druckfunktion 
folgende Entwicklung gefunden: 


1a 
p(én)=—Po = : sin mz acos mz &,cos ma E+ > + sinna£ cos na cosn a | 
2A ae : B No a (5) 
+ or > > Sin maacos mn &,cosma é sin NT COSNT -— cosnr | 
(m,n = 24,05 4) ) 
Daraus leitet sich dann mit Hilfe der Differential gleichung 


aE w\? _ p(&n) ms Eh 
Adw == (Fe oa f (X= a7) (6) 


folgender Ausdruck fiir die Durchbiegung der Vierpilzplatte ab: 


yaar [ll Loy we Ibe B No n 
Ui = Wo — ee ls Sain MILK COS EGE DORI Gs ales > =; Sin nt -cos nm —- cos n—- 
™ 


m 


| 
=> l ae ‘0 ae 
sears a nUntAl f payin MEA cos mn c, con mi s sm M7 eee 


= Wy + We + W, + Wey (m:, 3 == 2,4, 62), 


dessen wesentliche Umgestaltung unsere nachste Aufgabe bilden soll. 


4, Umformung der Reihen im Gebiet I. Wir 
wenden uns zunachst der Doppelreihe w;,, zu 
und wollen versuchen, die auf die Variable 7 
und die Laufzahl n sich beziehende Teilreihe mit 
Hilfe der Formeln (3) und (4) zu summieren. 


XY alt Wie man bei der Anwendung dieser Formeln 
SG sofort erkennen wird, ist es erforderlich, wie ja 


IN 


auch aus der gewoéhnlichen Statik bekannt ist, 
verschiedene Gebiete oder Felder zu unterschei- 
den, wobei wir uns zunachst beschranken wollen 
auf das von achsenparallelen Geraden begrenzte 
Flachenstiick, welches die Stiitzflache um den 
Punkt &) 7) enthalt und bis zur angrenzenden 
Stiitzflache reicht (Abb.1). Wir haben dann 


Gebiet I: 0 <n <yH—f8; 
Gebiet IIT: nn—B <n <m+F6; 
Gebiet III: 4) + B <n <A—n— 8B. 
Im Gebiet J kann man die Reduktion auf 
die Forme] (3) durch die Darstellung 


sin nal cos nN Wo cos n ma 
A a A 
IL Wo ioae : 
= [in lo + B+) + sin" (m+ B—M) 
wae pee - nse - nn 
Abb. 1. Feldereinteilung der Vierpilzplatte. sin 7 (% B 4 ”) Noe (mB) (8) 


erreichen. Wenn man dann auch auf die Funktionen f(y) der rechten Seite der Formel (3) 
die Operation 


PAD =F SZ Omt B+] +45 (ot p—n|—S]z (o—B-+)) 415 (m—B—n)}} ) 
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anwendet, so gewinnt man mit @ = mA aus den Ausdriicken 


fi=no Sinwy, f= Sin os —7) ie Oy Cn) oF 7 


folgende neue Funktionen: 


i f= Dy(fi) 
= 7 eAm[Sin ma (ny + B +n) + Sin ma(H + P—n)— Sin ma(n—f + )— Sin ma(n—B—n)] 
=ahm Sinm xB Coj mx) Coj man, (10) 
=D) = — Sit max Eo} mal —n Co) man, (11) 


F,=D,(f,) = Gin map| macy Sin man ofma( 5 — no) + men, Sin mas = ng Co} m a] | 

s (12) 
+ ma B Co} max B Co mas — no} Col man. 

Der von 7 unabhangige Ausdruck les — 7) fiihrt dagegen auf 


1 It New, B I 
Daeala l= 2A (Am)i* We), 


Dann sieht man aber bei Multiplikation mit — sin maacos ma &é,cos ma € und dem kon- 
m 


stanten Faktor, da man wieder auf die Reihe — we kommt, die sich daher mit der ersten 
Reihe w; aufhebt. Nach Abspaltung dieser Reihe — we bleibt daher von w;, folgende Rest- 
reihe wbrig: 


ee ae Dy sin m 7% & Cos Be Sooke ES med Gin m x B Cof mx no Coj xy 


m=2, 4,6 m? Sin m x [ Sin mas 
ae Coj max(-F — ne) C0} may (2 Sinmazp —mxzf Coj mxp)— may Sinmap Sinman]) 


+2mxan Sinma p Sin ma — na} Co} mary} 
y 


a* 7), sin ma cos m m7 

Pins Soe a somee EE ¥i(t ) - (14) 
m=2, 4,6 m°? Sin ae 
Wenn man nun zur Umgestaltung der zweiten Reihe w, = w, die Gleichung (8) und dann die 
D,-Operation auf die rechte Seite der Formel (4) anwendet, so ergibt sich an Stelle der in w, 
auftretenden Summe eine ganze rationale Funktion von 7. Betrachten wir z. B. die dritte 
Potenz von o auf der rechten Seite von (4), so kommt durch Anwendung der D;-Operation 


BIE sae) ata) -ponre a rr. a 


Fiihit man diese Operation auf die iibrigen Potenzen 0, 0%, 0° aus, so entsteht schlieBlich fiir 
w, folgender Ausdruck: 


(10,)7 = 20% [gt + 22 9? — 6 ngh*(A— M0) + 18 (A—mo)® + 2 BCH? + 8 — Ano) = Sy Hel) (16) 


Das ist natiirlich eine partikulare Liésung unserer Differentialgleichung, ndmlich die Darstel- 
lung der Biegeflache eines in der x-Richtung unendlich langen Plattenstreifens von der Breite b, 
der gleichmaBig belastet und langs zweier symmetrisch zur Mittellinie gelegenen Streifen von 
der Breite 2a 8 gestiitzt wird (Abb. 2). w, kann auch als Ordinate der elastischen Linie eines 
entsprechend belasteten linearen Balkens gewonnen werden. Im Ganzen haben wir also folgende 
Durchbiegungsfunktion unserer Vierpilzplatte im Gebiet I (0 << <1 — B): 


12* 
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7} 
Wi = Wy 4 Pe [Wea Ono J?(A — 19) + nglA mo)? + 2 B?(n? + 76 —AMo)] | 
& | (17) 


Po ad ee ee 
I 0 . 


si 5 
ay m=2,4,6 m? Sin m a> J 


5. Die Durchbiegung in den Gebieten III und II. Betrachten wir zunachst das Gebiet IJT, 
so miissen wir hier ausgehen von der Darstellung 
No 


: B 7] 
sin ny Mec nu 
sin N71 -5- cos N-- cos 5 


i 4 4 ; . 
= |sin 2 (mo + B+) — sin (n— 140 — B)— sin (n— no + B)— sin (n+ No p)-| 


L260 ofa 


Abb. 2, Plattenstreifen. 


Man sieht aber, das der rechts stehende Ausdruck sich von dem entsprechenden in (8) nur 
durch Vertauschung der Groen 7 und 7 unterscheidet. Man schlieBt daraus, daB auch der 
vollstandige Ausdruck fiir die Durchbiegung, d.h. wzzz aus wz dadurch hervorgeht, daB man 7 
und 7) mit einander vertauscht. Wir werden den Wert in der Zusammenfassung angeben. 
Etwas schwieriger und weniger symmetrisch gestalten sich die Rechnungen fiir das die Stiitz- 
flache enthaltende Gebiet II (yy) — fh <9 <n)+). Man mu8 hier ausgehen von der 
Gleichung 


aos 
= 


% B 7 7 LS ieemererce 5 EE 
sin n 1 cos nT cos MMs = ‘sin 7 (74 sin (no + B — n) | 
m7 (19) 


r sin —- (1 — to + B) sin (1 t 9 — P)| = Drz (sin n V)-| 


Wendet man dann auch auf die rechte Seite von (3) die Djz-Operation an, so erhalt man nach 
etwas umstandlicher, aber nicht schwieriger Rechnung folgenden Ausdruck: 


(ws) r= 


Po aa sinmaacosma&cosma& 
— ENG B a! > ferns en, Yir (1 Mo P) 
mit m=2,4,6... m> Sin m a 


eo Sin m zB Coj ma Hy) Cof may 


2einma(y —"o~ B) Co) ma 7 


Sin m me 
=? Sint mar(ro— f) Go} ma(F ns m (H+ ) Co) ma(z—m—B) Cojmay (20) 
+ m z6(j — f) Co} m (mH — B) Co} ma(y—n 


_ mo n| in mala — m—P) Sin man — Gin ma(n — B) Sin mals _- ) : 

Subtrahiert man nun die beiden Funktionen Y7 und Y7, so ergibt sich folgende Beziehung: 
; A es 

Yu= Y;-+ ©in mT [m x(n — No + B) Sin ma(n — H + B)—2 Co] ma(n— y+ )]. (21) 


Ebenso ergibt sich, wenn man statt Y; den durch Vertauschung von » und 7) aus Yyz ab- 
geleiteten Ausdruck Yj77; verwendet, 


; A oer “ 
Yo = Yin + Sinm Uy [m x(n — Ho — B) Sit m x(n — Ho — B) — 2 Co} max(y— ny— B)]. (22) 
Daraus erkennt man dann, da an der Grenze zwischen I und II (yn = — P) 


Y77 = Y;— 2 Sin mae 
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und an der Grenze zwischen IJ und III fiir yn = yn) + f dieselbe Differenz entsteht. Um also 
einen stetigen Ubergang von w, an beiden Grenzen zu erzielen, mu man zu Yj; noch die von 7] 


unabhiangige GréBe 2 ©in ma. hinzufiigen. Welche Auswirkung das fiir die weitere Unter- 


scheidung innerhalb des Gebietes II hat, werden wir weiter unten besprechen. 


Zunachst miissen wir noch die Funktion w, den Bedingungen des Gebietes II anpassen, 
Durch Anwendung der D;7-Operation entsteht folgender Wert: 


a* { 


an aye G pA? nf? + 19 (A— No)? — 6 9 (A— Mo) + 2 B?(? + 76 An) A (1 No py 


(W) 77 = 


= anal? Ay(7 %0) s ica | Pe p)'| = — saa G Hyr1 (4 o) a (7— No a) : (23)) 


Der Klammerausdruck ist bereits so umgestaltet worden, da® der erste Teil die ent- 
sprechende, im Gebiet I, bzw. III giiltige Funktion H;, bzw. Hy enthalt. Aus dem Zusatz- 
glied sieht man dann ohne weiteres, dafs an beiden Grenzen des Gebietes H =n — B und 
N =m + f ein stetiger Ubergang sowohl der Biegewerte als auch ihrer ersten Ableitungen 
nach 7 stattfindet. 


6. Aufteilung des Gebietes II. Die in Ziff. 5 besprochene, aus Stetigkeitsgriinden geforderte 
Hinzufiigung der Konstanten 2 Gin mum , die auch im Nenner der Summe vorkommt, ist, 


wie man sieht, gleichbedeutend mit der Abspaltung der auf die Variable é sich beziehende Reihe 


Uy Sas 


Po aa sin m 4X COs m % Eo cos m mF 
INapa >, me 
m=2,4,6 


(24) 


Auf diese Reihe kann man aber wieder die Formel (4) anwenden und erhalt dann, wenn man die 
drei Gebiete unterscheidet 


TO oe dd ey eR ey 4 os IE ey pa f< 1 o£, 
in vollstandiger Analogie zu der Umwandlung der w,-Reihe w folgende Werte fiir w: 


(1)rr, = Boag [Et + EP 6 E082 (1— fo) + £8 (1— Fo)? +2 (E+ BN = SESE Gr, | 


(wn, = 85 [a7 — 6)? — 0 6 20 6) 4 2 (EF) 4 =o Gir,» } (25) 


Po ata 
(i) = G6 wpa |* oH 


1 Al eseaoy se! 1 A 
rE Sot 0] = a6 rag [2 On, — gE fo — a). 

7. Zusammenstellung der w-Funktionen in den verschiedenen Gebieten. Nach den bis- 
herigen Aufstellungen sind wir in der Lage, die Funktionen w fir die Durchbiegung einer Vier- 
pilzplatte in den fiinf Gebieten J, III I,, II,, I], ibersichtlich zusammenzustellen. Wir haben 


a*) sin maxcos maz &,cos ma € x;(m) 
Wy — Wo pare wn Ht (7 No) 1 yy eas > = xr (pe) 
ai ere: m®> Stim 2 — 
uP 2 
ey at) sin mwa cos mz &) cos mm § y/(m) Gehele 
= Deeded III 
WIT = Wo on ath No) 4 a Ba > eee 0 
Hi 2 
Poath Dye Po @as sin mmx cos ma fy) cos mam § y{7) ( ) 
wry = W, + G Hy > —— TF (Na) 
Sez 8 96.N 6. 1's "24. 4 Nap x? SM eta (26) 
i 2 
jaye (ee wor H pati sin maacos mz & cos m a & Yn No) s 
Wir, = ®o 96 NB Il, I 24.N II CNT, an eran Ls 
fh 2 
Po atd at Po ata sin m 7X Cos Mm 7% Fy Cos Mm 7X g YO) (4 
wir, = + 96 vB Het can ME ENB a 1 (7 No) 


mee A 
7 m°? SUL mm Uae 


(ni == 2, 4,70...) . 
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Die eingefiihrten Funktionen haben dabei folgende Bedeutung und Beziehungen zu einander: 


yin) _ maa Gin m x B Cof m x Hy Coj m aH , 
: 1 
Sin m a> 
+ 260] mal-5 ~ nn) (Gof man (2 Sin maz B — ma f Coj mz f) 
— ma Gittm af Sin marr} -+ 2 ma yy Sin mae B Sin max 5. — me) Cof mary. 
yin) ma Gin m x B Coj man, Coj may 
II 7 
Sin m x > 
+260) ma(3 —n) [Gof m a n9(2 Sin ma B — ma B Coj max f) 
: : A 
— mary Sin max Sin max yy) +2 mary Sin ma f Gof mary q Sint ma 5 — ): an 
| 
Yip =YP) + Sin maZ [ma (yo +B) Sin max (q—mo + f)— 2 Cf ma (n—Mo + BI 
= Vf] + Sin ma [nx (y—no— f) itt max (q— no— B)—2 Co} mae(n— no— BI 
Hy = 78 + A? 4? — 6 (A — 199) 4? + 9(A — No) + 2 BC? + 96 — A No) » 
Hy = 9 + 2? ng — 6 n(A— 1m) ng + 9? (A — 0)” + 2 BP (n? + 179 — An) ; 
A 
B Hie = B Hr “(1 Noe BP) B Ay — io =P) 
Grr, = 54 + 6? — 6 Eo(1 — So) & + S5(1 — So)? + 2 a? (6? + £5 — fo), 
Gy, = E54 Fe 6 S(1— FG et — 8)? 2 eee eg a) 
aGr, = 0 Grr, : (€ — & +a) = o Grr, a E&y— a). ; 


Aus Symmetriegriinden gelingt es nun leicht, auch fiir die noch fehlenden acht Felder (vgl. 
Abb. 1) die entsprechenden Funktionen fiir die Durchbiegung aufzustellen. So erhalt man wy 
fiir das dem Felde I gegeniiberliegende Feld V,ebenso wy aus wy; baw. aus wy durch Ersetzung 
von 7 durch A — yn. Ebenso ergibt sich w in den noch iibrigen beiden Feldern von II, namlich 
TI und IT, aus wy, und wz, dadurch, da®B man & durch 1 — £ ersetzt. Damit ist die Durch- 
biegung w in allen 13 Bente der Vierpilzplatte festgelegt und gleichzeitig ist auch der stetige 
Ubergang an den Grenzlinien zwischen diesen Feldern beriicksichtigt. 

Die Form der Darstellung kann noch dadurch etwas abgeandert werden, daB man statt m 
die Ganze Zahl 2 n einsetzt, wobei dann die Laufzahl n die ganze Zahlenreihe von 1 bis co 
durchlauft. So erhalt man z. B. fiir das Gebiet I die Durchbiegung 


i — 


Pot Int 4 22? — 6 (A — to) 9° + 2 B? (4? + m6 — Amo) + 6 (A — 10) 
_ poatd sin 2n 7a cos 2n xa &) cos 2nzé { nad Sin 2n mB Coj 2n xy Coj An aH 
aN cE n> Ginn wd \ Sinn wA 


n 


+ Goj anal = ny [Coj2nx7(2 Sin2nzp—2nxP Coj2nzp)—2nynSin2nzp Sin2naz7] 


+ 2n my Sin? nae B Sin2 na(— my) Col? nae (ie ase (28) 


Mit « =$ =0 geht dieser Ausdruck in Ubereinstimmung mit einer fritheren Aufstellung 
iiber in 


epee Fy hays ee Lae Sao Po ah cos2nz&)cos2nz& ) 
1 = (n! TA 19 Wag (A = Tho) Mg (A= Ne) tm ag es a 

Tein 

naAbojf2nanboj2nan , i 


Ginn wd Cof2ua(F —m] (Gof 2nay—2nay Sin 2 nxn) 


+2n2a) Sin2n «(> a) Coj2nzx z : 


— 
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Dieser Fall « = 6 = 0 gibt noch Anlaf zu einigen Bemerkungen. Man sieht, da® in den For- 
meln z. B. fiir das Feld IIa, das mit der Grenze zwischen I und III zusammenfiallt, eine schein- 
bare Unbestimmtheit auftritt, die sich aber auf Grund unserer bisherigen Entwickelungen 
leicht auflésen 1aBt. Bedenkt man namlich, daB fiir B = 0, Yyp = Y; — 2 Sin sek wird, 


2 
so sieht man, daB auBer der mit Y; geschriebenen Reihe und dem yenthaltenden Glied w, noch 


folgende Foden Glieder iibrig bleiben: 


at} at) cosmaécosmaé 
a Wp Me ONG se os (39) 
m=2,4.06 

die aber zusammen Null ergeben, wenn man die Formel (2) darauf anwendet, sodaB das Er- 
gebnis wieder mit wy fiir 7 = 7p iibereinstimmt. 

Unter den vielen weiteren Sonderfillen, die in den 
allgemeinen Formeln enthalten sind, erwahnen wir 
noch den Fall a = &,= 1/4, bei dem je zwei der von 
der x-Achse gleich weit abstehende Stiitzflachen zu- 
sammenriicken. Es bleiben dann zwei Stiitzstreifen 
iibrig, die sich in der x-Richtung iiber die ganze Platte 
erstrecken (Abb. 3). 

Da aber in diesem FalJle das Produkt 


eae 2 2 
fiir gerade m gleich Null wird, so verschwindet auch 
die gesamte Summe des letzten Gliedes in dem Aus- 
druck fiir w. Ferner reduziert sich aber die Funktion 
Gry, auf eine Konstante. Es ergibt sich also, daB nur 


die Lésung fiir einen Plattenstreifen iibrig bleibt, der 
die Breite b hat und in der x-Richtung ins Unendliche 
erstreckt. 


8. Die Berechnung der Momente. Mit Hilfe der ent - 


1 
wickelten Formel!n fiir die Durchbiegung kann man die Abb. 3. Der Fall & =a =~. 
Momentensumme ohne Schwierigkeiten ableiten. Es 

dirfte geniigen, wenn wir die Ergebnisse der Rechnung iibersichtlich zusammenstellen : 


2} sin maz a&cos ma &,cos mx & (mm) 
M; = CO eae : Py”) ( Mo) » 
a Bm m? Sin m 0 
2 SED iG sin maacosma &cosma€& y(m 
Myr =P tyr (n No) Pir (7 No) » 
: Dafa m® Sin m x 
a” Po aA sin maacosmaé,cosmaé& () ( 
Myr, = Be * gir,(6 £0) +2 har(n) — Pegs > - Pit'(n No): 


m m® Sin m x -- 


a?) nmazacosma &,cosma (m 
Mir, = ae ar £5) +P har 0) ee s 5 E pt )( Cone 


m?® Sin m 1 


7h a? py aA sinmaacosma&,cosmaé (mn) ( 
My =a air,(& £9) +2 hyx( No) ie : = Pir 


c 


(7 No) 
m> Sin m ms 


(Gis =7, 4, O57.) 
Dabei ist zu setzen 


; a 
PY} = 2 Sin Pe 
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! >: ‘, A 
P\?) — — Cof m nz —") Sin mx(m — B) — Sin mas —no—B) Co) may , 


5 , ap <a A 
= Pp”) — ©inm 0 Co) max(n— yy + B) = P= Gin m7 Co) mx(n — Hy — fA), 
(32) 
Ie a5 ese i 1 i ; 
hy =A5 6 A? - 3 P—=i—7 > hpr=A 6 A? 3 pire iG 
B hy = B hy — 3A(q — ny + B)? = 8B har — 3.A(y — Hy — 8)? 
> ee oma 
Sil, — & a 6 3 on eA Ge ) SII, = eS 6 3 a? es (ss ? 
agn,= + 8, — 3(§ — § + 4)? = 4% 8, — 3(F — &o — 4)? j 


Das Einspannmoment am Ende des Plattenstreifens, der auch als Gesamtlésung fiir £5 = a= 1/4 
erschejnt, erhalt man aus dem ersten Glied von M; mit 7 = 0. In gewohnlichen Koordinaten 
hat man, wenn man die halbe Breite des Stiitzstreifens f a = d setzt, zu 


Pp 1 tl 
My = = b ¥o 5 b? ; d? — y2l. (33) 
Zu dem Fall « = § = Oist noch zu bemerken, da hier wieder eine scheinbare Unbestimmtheit 
auftritt, deren Auflésung zugleich eine gewisse Bestatigung der Rechnung liefert. So haben wir 
im Gebiet II,, das sich natiirlich auf die Linie reduziert, auBer den in I auftretenden Gliedern 


wegen (Pr7)g-9 = — Sim ae + (Pr)g-o folgende Zusatzglieder : 
2} iL ay) cosmam & c¢ naeé 
a (fo : ee a L 1 > s oe mots (34) 


Setzt man cosma & cosa é& ae [cos ma (&)— €) + cos mz (&)-+ &)] und beachtet die 
leicht aus (2) durch Differentiation zu gewinnende Fouriersche Entwicklung 


cosmnag@ 2 
m? 24 


[l1— 60+ 6¢?], (35) 


m=2,4,6... 


so stellt man ohne weiteres fest, daB die Glieder (34) sich aufheben oder daB der Wert von M, 
" = % tibrig bleibt. 

Damit sind alle Voraussetzuggen entwickelt fiir eine weiter gehende Untersuchung der mit 
der Vierpilzplatte zusammenhangenden Fragen. Insbesondere soll an anderer Stelle gezeigt 
werden, da man auch in diesem Falle funktionentheoretische Methoden, z. B. die Theorie der 
Theta-Funktionen benutzen kann, um die Momente der Vierpilzplatte iibersichtlicher dar- 
zustellen. 


(Eingegangen am 3, August 1953.) 
Anschrift des Verfassers: Prof. Dr. W. Miiller, Miinchen 9, TaubenstraBe 2/III. 
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Schmiermitteldruck und Randverformungen des Rollenlagers*. 


Von J. Dorr. 


1. Aufgabenstellung. Die Frage nach der Lebensdauer von Rollenlagern setzt die Kenntnis 
der in den Rollen, der Lagerschale und der Welle wirkenden Schub- und Normalspannungen 
voraus. Zu ihrer Berechnung mu man die an den einzelnen Teilen des Lagers wirkenden 
Randkrafte kennen. Beim geschmierten Lager berithren sich die Oberflachen der Rollen nicht 
unmittelbar mit denen der Welle und Schale, wenn die Welle in Rotation ist. Dann zwangt 
sich namlich das Schmiermittel zwischen die Oberflaichen, so da® zwischen den Rollen und 
der Welle bzw. Schale Spalte entstehen, die in gewissen Breitenbereichen liickenlos vom 
Schmiermittel ausgefiillt sind. Die auf die Oberflachen der Lagerteile wirkenden Randkrifte 
sind dann also identisch mit den Randkraften des Schmiermittels. 

Bei ziemlich weitgehend giiltigen Voraussetzungen lassen sich die im Schmiermittel wirken- 
den Spannungen berechnen, wenn man in dem liickenlos vom Schmiermittel ausgefiillten 
Bereich die Spaltweite als Funktion des Ortes kennt. Nimmt man die Spaltrander als unver- 
formbar an, so ist die Spaltweite bekannt, sofern man die Mindestspaltweite vorgibt. 

In Wirklichkeit haben aber die Randkrafte Verformungen der Spaltrander zur Folge. Das 
heiBt, die Spaltform, aus der sich die Schmiermittelspannungen berechnen lassen, hangt ihrer- 
seits von der GréBe und Verteilung der Schmiermittelspannungen ab. Das fiithrt mathematisch 
dazu, daB die Spaltweite erst auf dem Umweg iiber eine nicht lineare Integralgleichung be- 
stimmt werden mu. Eine strenge Lésung dieser Integralgleichung diirfte selbst unter 
den starksten zulassigen Vereinfachungen aussichtslos sein. Auch der numerischen Behand- 
lung stellen sich noch so erhebliche Schwierigkeiten in den Weg, daB bei einem ersten Liésungs- 
versuch vereinfachende Annahmen geboten erscheinen. Fir die folgende Untersuchung 
werden unter anderem fo]gende Annahmen gemacht: 

a) In den Lagerelementen werden ebene Verformungszustande vorausgesetzt. 

b) Die Strémungsvorgange im Schmiermittelspalt werden als laminar und eben voraus- 
gesetzt. 

c) Die Zahigkeit des Schmiermittels wird als unabhangig von Druck und Temperatur 
angenommen. 

Die letzte Annahme ist fiir die Anwendbarkeit des hier mitgeteilten Verfahrens nicht 
unbedingt erforderlich, vereinfacht aber die numerischen Rechnungen. 


2. Die Berechnung des Schmiermitteldrucks bei vorgegebener Spaltweite. Strémungsvor- 
ginge von Flissigkeiten gehorchen der Navier-Stokesschen Gleichung. Betrachtet man die 
laminare Strémung eines Schmiermittels in einem engen Spalt, so lassen sich die Tragheits- 
krafte gegeniiber den Zahigkeitskraften vernachlassigen. Wenn wir auBerdem noch das 
Schmiermittel als inkompressibel voraussetzen, also div ) = 0 fordern, reduziert sich die 
Strémungsgleichung auf die Form 


gradg =7 Ab. (1) 


Hier ist q der Schmiermitteldruck, 7 der Zahigkeitskoeffizient und ) die Strémungsgeschwin- 
digkeit. 

Auch diese recht einfach erscheinende Beziehung zwischen Druck und Geschwindigkeit 
bietet der mathematischen Auswertung noch erhebliche Schwierigkeiten. Mit einigen zusatz- 
lichen Annahmen, die beim engen Spalt im allgemeinen erfillt sind, hat Reynolds von (1) 


* Die vorliegende Untersuchung wurde mir erméglicht durch ein fiir zwei Jahre gewahrtes Stipen- 
dium der Deutschen Forschungsgemeinschaft. Als sich herausstellte, daS allein mit formal mathemati- 
schen Methoden das Problem nicht zu bewiltigen war, erlaubte eine groBziigige Spende der Deutschen 
Shell A.-G. die Hinzunahme numerischer Untersuchungen. Die Deutsche Forschungsgemeinschaft 
gewahrte danach fiir die Zuendefiihrung der numerischen Rechnungen noch eine Sachbeihilfe. Auch die 
Pametrada Research Station in Wallsend/Northumberland férderte die Durchfiihrung der umfangreichen 
numerischen Arbeiten in wirksamer Weise. Fiir alle dargebotenen Unterstiitzungen danke ich herzlich, 
Mein Dank gebiihrt auch noch den Herrn G. Niemann, W. Peppler, G. Vogelpohl und A. Cameron, die mir 
bei Besprechung technischer und physikalischer Fragen vielfaltige Anregungen gaben, 
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ausgehend eine noch wesentlich einfachere Beziehung hergeleitet. Bei unserem zweidimen- 
sionalen Strémungsproblem kommt man aber auch auf sehr bequeme und anschauliche Weise 
direkt zu der gewiinschten Beziehung. Wir betrachten dazu den in Abb. 1 dargestellten Spalt. 

Der untere Spaltrand sei eben, der obere Spaltrand habe vom unteren den Abstand 
H = H(x). Das Koordinatensystem sei so gewahlt, da® H(0) = Hnin ist. Der Spaltbereich, 
der liickenlos vom Schmiermittel ausgefiillt ist, sei links durch x,, rechts durch 4%, begrenzt. 
An der Stelle x = 0 habe der untere Spaltrand die substantielle Geschwindigkeit V,, der obere 
Spaltrand die substantielle Geschwindigkeit V, nach links. Wenn also V, = FV, ist, rollt der 
obere Kérper schlupflos nach rechts auf dem unteren ab. 

In dem Spaltbereich x, <x < xs, der liickenlos vom Schmiermittel ausgefiillt ist, sei 


fH(x)|<1 und (m— %) > Hnin- (2a) 
Dann gilt fiir x, < x < % bis auf vernachlassighar kleine Bereiche 
0q oq |. 07x 0? vx 
Poles © dat | <| ay |: a 
Dadurch reduziert sich (1) auf die Naherungsbeziehung 
2 
dq ae 0 Ux (2c) 


dx " dy? * 
Fir laminare Strémungen im Parallelspalt gilt (2c) streng. Das Geschwindigkeitsfeld ]4Bt sich 
in diesem Sonderfall ebenfalls streng beschreiben und zwar durch einen quadratischen Ausdruck 
der Form 
v,=atby+tey’*. 

Bei Vorgabe der DurchfluBmenge 
lassen sich mit den Randbedingungen 
alle drei auftretenden Faktoren a, b, c 
bestimmen. Da unser Spalt wegen (2a) 
nur wenig vom Parallelspalt abweicht, 
wahlen wir fiir die Geschwindigkeits- 
Abb. 1. Geschmierte Rolle auf ebenem Untergrund. verteil ung den Ansatz 

v, = a(x) + y b(x) + y?e(x). (3a) 


Die Randbedingungen Jauten, wenn wir wegen (2a) quadratisch kleine Fehler zulassen, 


U,V tun yo ye 0, =), firsy= 02 (3b) 
Beachten wir noch, daB wir Inkompressibilitat, also 
divy = 0 (3c) 


vorausgesetzt haben, so lassen sich alle drei Funktionen a (x), b (x), c (x) bis auf einen offenen 
Parameter bestimmen. Dieser Parameter, den wir hy nennen, steht in Zusammenhang mit der 
DurchfluBmenge. Das Ergebnis der hier nicht wiedergegebenen elementaren Rechnung Jautet 


te =— Vat (Ve— VY gg + 8(V + VA Hla) — gly. 3d 
H (x) H3 (x) (3d) 
Damit nimmt (2c) die Form an 
d ides 
Te = 6n(Vi+ %) Fi (4) 


Negative Drucke sind im offenen Rollenlager nicht méglich, da der Schmiermittelfilm 
durch Unterdruck auseinander gerissen wird. Der Bereich des liickenlos vom Schmiermittel 
ausgefiillten Spaltes wird also durch die Randbedingungen 


q(%a) =q(%) =O (5a) 
festgelegt. Das ist aber nicht ausreichend zur eindeutigen Festlegung von q(x) bzw. v,x(x, ¥), 


da in den in Frage stehenden Beziehungen noch der offene Parameter hy auftritt. Druck- 
messungen an Rollenlagern ! zeigen, daB an der Hinterkante x, des Druckbereiches eine Art 


i Siehe z. B. W. Peppler, Druckiibertracung an seschmi j i itfli 
Forschungsheft 391 (1988), $2 und 8. ee g an geschmierten zylindrischen Gleitflachen. VDI- 
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AbfluShedingung zu beachten ist, die sich mathematisch durch die Forderung 


d 
= = fiir 55 c=) 6,0 (5b) 
ausdriickt. 
Die Randbedingungen (5a) und (5b) sind erfiillt, wenn man x, — -+ | x5| vorgibt, 
hy = H (xa) (6a) 
setzt und den Rand x, durch die Forderung 
_ 
H(xa)— H(z), , 
| ou = (6b) 
festlegt. Fiir den Druck findet man dann, wenn 7 als konstant angenommen wird, 
x5 
2 H(«a) — H 
EOS NU AN feels (6c) 


Die Lage des Randes x, wird durch die Starke der Schmierung, dem sogenannten Spalt- 
fillungsgrad festgelegt. Je starker die Schmierung ist, desto gréBer wird x, und auch %/|x,| . 
Beim unverformbaren Rollenlager ist H(x) bestimmt, wenn man H,,;, vorgibt. Dabei 


besteht ein eindeutiger Zusammenhang zwischen der Mindestspaltweite H,;, und der Gesamt- 
rollenlast 


P = f q(x) dx. (6d) 


Bei sonst gleichen Bedingungen ist 1 Pumgekehrt proportional zu H,,;,. Die Berechnung der 
Druckfunktion stellt beim unverformbaren Lager, weil die Spaltweite H (x) bekannt ist, 
lediglich eine Integrationsaufgabe dar. Nur die Bestimmung des hinteren Randpunktes 
Xq — f(x) bereitet gewisse Schwierigkeiten. W. Peppler? hat dieses spezielle Problem in 
seiner Dissertation unter anderem erschépfend behandelt. 


3. Berechnung der Normalverschiebung eines Rollenrandes bei vorgegebenen Randkraften. 
Die im Schmiermittel vorhandenen Druck- und Schubspannungen wirken als normale und 
tangentiale Randkrafte auf die Randflachen des Rollenlagers ein und rufen Normal- und 
Tangentialverschiebungen der Randflachen hervor. Die Spaltweite wird von den Tangential- 
verschiebungen der Rander nur quadratisch wenig gedndert. Jee 
Die Normalverschiebungen werden yon den tangentialen RSQ 


Randkraften nur quadratisch beeinfluBt. Deshalb brauchen 
wir nur den Einflu®B der normalen Randkrafte auf die nor- | 
malen Randverschiebungen zu beachten, um die elastische 


Verformung des Schmiermittelspaltes zu ermitteln. Zu ihrer RXQQAN 


Berechnung gehen wir von einem vereinfachten System (Abb. 2) GP bees Or ee 
aus. Wir nehmen an, der Rollenradius sei klein im Verhaltnis See ciache Walbeaaies 
zum Radius der Welle, was praktisch in den meisten Fallen 
zutrifft. Dann darf man die Kriimmung der Wellenoberflache und der Innenflache der Schale 
gegeniiber der Kriimmung der Rollenoberflache vernachlassigen. Wir betrachten deshalb zur 
Berechnung der Rollendeformation ein Ersatzsystem, bei dem die Rollen zwischen zwei eben 
begrenzten Halbraumen liegen. Die Rotation der Welle ahmen wir dadurch nach, da8 der obere 
Halbraum mit der Geschwindigkeit Vnach rechts bewegt wird. Dabei entstehen in den Spalten 
zwischen den Rollen und den Begrenzungsebenen der Halbraume Schmiermitteldriicke. Die 
auf eine Rolle einwirkenden Randkrafte sind aus Symmetriegriinden periodisch mit der 
Periode 7, wenn man der betrachteten Rolle Polarkoordinaten mit dem Rollenmittelpunkt 
als Ursprung gibt. 

Die Periodizitat der Randkrafte erlaubt die Einfiihrung einer speziellen Greenschen Funk- 
tion fiir die Radialauslenkung des Rollenrandes, mit deren Hilfe sich die Auslenkung fiir all- 


1 W. Peppler, a. a. O., Formel (17). 
2 A.a. O., Seite 2/4 und 12/13. 
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gemeine m-periodische Randkrafte berechnen 1a8t. Als Greensche Funktion wahlen wir die 
Radialauslenkung, die der Rollenrand erleidet, wenn zwei zum Mittelpunkt weisende, auf dem 
selben Durchmesser wirkende Einzelkrafte am Rollenrand angreifen. 

Die unter dieser Belastung auftretenden Spannungen in der Rolle sind bekannt. Sie lassen 
sich durch einfache Ausdriicke beschreiben, wenn man sie in Teilspannungen zerlegt, die auf 
drei verschiedene, in Abb. 3 erlduterte Koordinatensysteme bezogen werden. Diese Teil- 
spannungen sind ! 

il I 


2 / 2 ihy¢ Sp ee ee pe hy aie 
Oj. cos 0 , Cy = ee cos? . Ops Opa aan (7) 


y 


ay, 

Fiir unsere Aufgabenstellung miissen wir alle Teilspannungen auf ein gemeinsames Koordi- 
natensystem umrechnen. Wir wahlen dazu das System (r,y). In dem gemeinsamen System 
kénnen wir entsprechende Komponenten der Teilspannungen addieren zu 


Og — Ge ~) ? n= 0, (Tr, 2) 9 Tr, a Tals ) 


fz 
{0 

T,,9(R,v) =9 fiir alle g-Werte . 
Die Ausdriicke sind recht unférmig und werden deshalb nicht an- 
gegeben. Aus den Normalspannungen lassen sich die Radial verschie- 
bungen des Rollenrandes durch Integration gewinnen. Man findet 


R 
1— v 
G(y) = B | (—y2 5%) 
0 


Abb. 3. Koordinatensysteme Tie Integration, die miithsam, aber streng durchfihrbar ist, ergibt 


fiir eine durch zwei entgegen- 
gesetzte Randkrafte belastete 2 


Rolle. G(y) = 2 cos gin 


GAR, @) = 07 tur p+ 


' ee eek 

+ 2 ae sin gl (8) 

Die hiermit gefundene Greensche Funktion ist periodisch mit der Periode z, wie es sein mub. 
Die von einer Belegung normaler Randkrafte J/ (gp) = I] (p +z) hervorgerufene Radial- 


auslenkung des Rollenrandes ergibt sich zu 


u,(p) = | IT(s) Gp — a) do. (9) 


C—T 


YP 
ae wy 189 


Beim Rollenlager tritt der Schmiermitteldruck nur in zwei sehr kleinen Winkelbereichen 


PaSPSM, A+PSPSnt pm mit 9 —GK<n 
auf, und man braucht bei unseren Untersuchungen auch die Radialverschiebung u, nur in 
diesen kleinen Winkelbereichen zu kennen. Deshalb darf man in (8) |tg + 


£ und cosy=1 


setzen und den dritten Summanden in der Klammer ganz vernachlassigen. Wir verwenden 
also im Folgenden immer die Beziehung 


Pb 
u(p) =20 | IT (a) fin Meee 7 dx mit m<o<q- (10) 
Pa 


Den Spaltbereich 7+ @<ySa+y brauchen wir wegen der Periodizitat nicht zu 
betrachten. 


4. Die Integralgleichungen des geschmierten, elastischen Rollenlagers. Wir bleiben vorerst 
noch bei dem in Abb. 2 und 1 dargestellten Ersatzsystem, nehmen also weiterhin an, daGB der 
Untergrund eine ebene Grenzfliche hat. Weiterhin wollen wir aber noch annehmen, da8B nur 
die Rolle verformbar, der Untergrund also starr sei. Alle mit diesen Annahmen gefundenen 
Ergebnisse lassen sich nachtraglich sehr einfach auf den Fall umrechnen, daB auch der Unter- 
grund (d. h. die Schale oder Welle) gekriimmt und elastisch ist. 


' Vel. S. Timoshenko, Theory of Elasticity S. 104, New York und London 1934, 
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Wir gehen von den Winkelkoordinaten @ und « auf Lingenkoordinaten x und z itber. Fiir 
den schmalen Winkelbereich aS Py, bzw. Langenbereich x,<x< ~ kinnen wir 
auBerdem das Kreisbogenstiick der unverformten Rolle durch einen Parabelbogen ersetzen. 
So entsteht fiir die Spaltweite der Ausdruck 


x 


H(x)=s4 a | q(2) (in eat -{- 1 dz. (11a) 


Xa 


Bei der unverformten Rolle, also fiir Elastizitatsmodul E—> oo oder O (2) 0 iets sie Ta, 
Mit der bereits hergeleiteten Beziehung (6c) und der nicht holonomen Randbedingung 


Xp 
H( ta) — H(z) 
| OEE ie ==.) (11b) 


haben wir somit in (11a), (6c) und (11b) die formelmaBige Beschreibung der gestellten Aufgabe 


mVOr UNS. 


Wir haben also zwei gekoppelte Integralgleichungen, von denen die zweite nichtlinear ist. 
Durch Einfithrung von (11a) in (6c) kann man eine einzige nichtlineare Integralgleichung fiir 
die Druckfunktion q (x) erhalten. Umgekehrt ergibt die Einfithrung von (6c) in (11a) eine nicht 
lineare Integralgleichung fiir die Spaltweite H (x). Beide so gewonnenen Integralgleichungen 
sind aber so kompliziert, daf ihre direkte Untersuchung hoffnungslos erscheint. Die Auf- 
spaltung in zwei gekoppelte Integralgleichungen (lla) und (6c) ist also vorzuziehen, insbeson- 
dere deshalb, weil diese Darstellung sofort die Méglichkeit eines iterativen Lésungsverfahrens 
erkennen ]aBt. 


5. Losung der Integralgleichungen durch Iterationsyerfahren. Durch die Formulierung 
unserer Aufgabenstellung in Form der beiden gekoppelten Integral gleichungen (lla) und (6c) 
wird sofort folgendes iterative Lésungsverfahren nahe gelegt. Wir geben eine physikalisch 
sinnvolle Spaltweite H) (x) vor (zum Beispiel die Spaltweite des unverformten Rollenlagers) 
und berechnen hierfiir nach (6c) unter Beachtung der Randbedingung (11b) die zugeordnete 
Druckfunktion q(x). Mit dieser Druckfunktion berechnen wir nach (lla) eine verbesserte 
Spaltweite H, (x), wiederholen das Verfahren, um ein q, (x) und H, (x) zu erhalten und fahren 
in diesem Sinne fort, bis die Funktionen qg,_1(«) und H, (x) sich geniigend wenig von den 
Funktionen der vorhergehenden Iteration unterscheiden. Dabei brauchen wir noch einen 
Parameter, der die relative Starke der Verformung beschreibt. 

Bei einem ersten numerischen Versuch offenbaren sich drei Schwierigkeiten. Die Werte s 
und x, andern sich stark mit dem Verformungsparameter. Die Erfiillung der Randbedingung 
(11b) ist lastig. Die Berechnung des Integrals mit dem singuladren Integranden in Formel (11a) 
ist mithevoll und sowohl formelmaBig als auch numerisch sehr zeitraubend, da es in Abhangig- 
keit von x, das hei®t numerisch fiir viele x- Werte ermittelt werden muB. 

Alle diese Schwierigkeiten werden iiberwunden oder gemildert, 

1) indem man eine neue Koordinate € durch x =AC + x, einfiihrt, so da® der Abstand 
vom Beginn des Druckberges p(¢) = q(«(¢)) bis zu seinem Maximum immer den Wert 1 hat 
und der Beginn des Druckherges an der Stelle €¢ = 0 liegt, 

2) indem man zuerst h’(¢) berechnet und dann daraus streng oder angenahert h(¢) ermit- 
telt. Dabei ist h(¢) = H(x(¢)), 

3) indem man die Druckfunktion p(¢) durch Reihen von solchen Elementarfunktionen 
darstellt, daB die ihnen zugeordneten Integrale zur Ermittlung von h’ (¢) sich streng berechnen 
lassen. ; 

Wir wollen nun diese drei MaBnahmen in unserer formelmaBigen Darstellung beriick- 
sichtigen und setzen dazu 


pop, Fede ees, abt alepO©; HAC + =)=h6). (12a) 
Dabei soll also A so gewahit sein, dah 
p (0) =p'(1) =0 ist. (12b) 
AufSerdem fiihren wir noch die Bezeichnungen ein 


Vee Va Ves = w(E; PB, V, 7, ¥) (Verformungsparameter). (12c) 
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Mit diesen Transformationen und Bezeichnungen ]4Bt sich unser Problem wie folgt formulieren: 


‘ ‘ K 
Mah 2 4 L2¢+bu al (13a) 
c 
h2(0) __ [ [h®(0) 2 (0) 
Ogavi=| ee | Hes 13) 


Die durch das Integral 


ch 
Rie es dé (13c) 
3 


saa 


definierte Funktion bezeichnen wir als Verformungsfunktion. Der Verformungsparameter 
faBt alle GréBen zusammen, die auf die relative Spaltverformung Einflu8 haben. Wie sich 
u=u(E, P, Vin, v) aus den Betriebs- 
daten des Rollenlagers ermitteln 1aBt, 
wird spater erladutert. 

AuBer dem Parameter yw treten in 
(13a) noch die beiden neuen Parameter 
A und B auf. Diese lassen sich bei vor- 
gegebenem - Wert und vorgegebenem 
¢,- Wert aus den Bedingungen 


p'(0)=p(1)=0 und p(%) =0 
a ec ae el es festlegen. Das bedeutet, daB man A 


0 if z 
Abb, 4. Druckkurve iiber wirklicher und bezogener Langenkoordinate. und B aus den Bedingungen 
h(0) = h(1) und. p(G)—=—0- (lady 


bestimmen kann (wobei vorausgesetzt wird, daB h (¢) stetig ist). Uber das ,,StreckungsmaB“ 
2 braucht man sich erst Gedanken zu machen, wenn die numerischen Ergebnisse technisch 
ausgewertet werden. 


Die erste der beiden Bedingungen (13d) erlaubt die unmittelbare Berechnung von A, wenn 
der Verformungsparameter 4, und die Verformungsfunktion K (¢) vorgegeben sind, wie es bei 
dem Iterationsverfahren der Fall ist. Es ergibt sich 

1 


A=—1—4. Kw) dt. (14a) 


0 


Die zweite Bedingung p(¢,) = 0, die zur Festlegung des Parameters B fiithrt, 14Bt sich nicht 
unmittelbar, sondern nur durch Probieren bestimmen. Nachdem man A durch (14a) berechnet 
hat, wird aus (13a) mit einem geschatzten B- Wert 


ig 
h(¢) 
aes Bf E b2E + A K() dé (14b) 
und damit das Integral (13b) fiir den vorgegebenen Randwert ¢ = ¢; berechnet. Ergibt sich 
ob 

h?(0) —-h2(0) ia | 
| i) Wie al d¢=0, (14c) 
0 | 


so haben wir den richtigen B- Wert. 


Mit dem so gewonnenen A- und B-Wert kennt man nach (14b) die ,,bezogene Spaltweite‘ 
h (¢)/h (0) und kann damit nach (13b) die ,,bezogene Druckfunktion“ 


h?(0) 


PO) ETA (15a) | 


berechnen. Durch Normierung gewinnt man daraus 


P(S) = po(S) + Ap(2), (15b) 
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mit POF) and p'(1) = pp (1) = 0. 
Dabei ist py (¢) die normierte Druckfunktion der nicht verformten Rolle. Wir bezeichnen 
Ap (¢) als Druckdifferenz. 
Die der normierten Druckfunktion (15b) zugeordnete Verformungsfunktion bezeichnen 
wir mit 
fb 
K(¢) = [ PEs dg = Kyle) + AKC). (15c) 


0 


Wie sich dieses unangenehme Integral numerisch berechnen laft, wird spater beschrieben. 
Wenn K(C) in erster Naherung bekannt ist, wird das Verfahren in der angegebenen Weise 
wiederholt, um h(C)/h(0) , p() und K(é) zu verbessern. Wie die Funktionen h(¢), p(é), 
K(¢), das heiBt die Werte h(0), uw, A aus den Betriebsdaten der Rolle berechnet werden, 
wird ebenfalls spater angegeben. 


6. Numerische Durchfiihrung des Iterationsverfahrens fiir das stark geschmierte Rollenlager. 
Wir beschranken uns hier auf den wichtigsten Sonderfall, nimlich den des stark geschmierten 
Rollenlagers. Bei starker Schmierung und den technisch iiblichen Rollenlasten ist das Ver- 
haltnis x;/|x.| und somit auch der Wert & sehr groB. Dadurch ist es erlaubt, in den Bezichun- 
gen (13) und (14) den Fall starker Schmierung durch den mathematischen Grenzfall ¢,—» oo 
darzustellen 1. Nur dieser Grenzfall wird von jetzt ab behandelt. Dieser Grenzfall ist noch aus 
einem anderen Grund wichtig. Der EinfluB® des Spaltfwllungsgrades ist namlich in weiten 
Grenzen sehr gering, wenn man von extrem kleinen Spaltfiillungsgraden absieht. Das hat ins- 
besondere W. Peppler durch Rechnungen beim unverformbaren Rollenlager gezeigt (Seite 13 
der genannten Veréffentlichung). Das bedeutet bei uns, daB der fiir €¢,—>co berechnete Druck- 
verlauf auch fiir endliche ¢,- Werte (¢,> 1) eine gute Naherung darstellt. 

Mit den Ergebnissen von Peppler findet man fiir diesen Grenzfall als normierte Druck- 
funktion (bei 5-stelliger Genauigkeit) 
zt | 


il u 
Po(C) Be = 7,59566 (1 + u2)2 1+ u? 


; arc tg u (16) 
mit 
u = 0,95026 (€ — 0.5). 
Die zugeordnete Verformungsfunktion K,(¢) 14Bt sich nicht formelmaBig streng berechnen, 
wohl aber ihre erste Ableitung Ko(¢), aus der sich durch numerische Integration K,(¢) selbst 
berechnen 148t. Man kann auch direkt K,(¢) durch Reihenentwicklungen berechnen, die im 
gesamten Bereich 0 <¢ < cokonvergieren. Davon wird im folgenden Abschnitt die Rede sein. 
Uber den Variabilitatsbereich des Verformungsparameters «4 1aBt sich unmittelbar eine 
allgemeine Aussage machen. Bei unverformbarer Rolle ist ~ = 0. Die starkste relative Ver- 
formung tritt im Hertzschen Grenzfall auf. Dann beriihrt die Rolle den starren Untergrund 
auf einem ebenen Flachenstreifen. Die normale Randverschiebung ist demnach iiber dem 
Beriithrstreifen eine Parabel, die Verformungsfunktion K(¢) also eine durch den Punkt 
= 1 hindurchlaufende Gerade negativer Steigung. Da im Hertzschen Grenzfall auBerdem 
h (£)/h (0)—>1 geht, folgt aus (14b), daB uw—> 2 geben mufs. Es gilt also 


Qe (17) 


Das Iterationsverfahren beginnt naturgema8 mit einem relativ kleinen yw-Wert, z. B. 

= 0,5. Die numerische Berechnung von A auf Grund von (14a) ist problemlos. Bei uns 

wurde dazu K(¢) in den Punkten € = 0 (0,05) 1 berechnet und dann nach Simpson integriert. 

Bei der Ermittlung von B nach (14c) ist es lastig, daB ein von 0 bis oo gehender Integra- 

tionsweg vorliegt (wegen ¢, = co). Die Funktionen K(¢) und A(¢)/h(0) wurden numerisch 

nur fiir die Punkte 

€ = 0(0,05) 1(0,1) 2,5(0,5) 4 5 6 (18) 

berechnet, so da® nur in dem Bereich 0 <€ <6 nach Simpson integriert werden konnte, 

wahrend fiir den AuBenbereich € > 6 eine geniigend genaue Abschatzung vorgenommen 
werden muB8te. Dabei wurde folgender Weg gewahlt. 


1 W. Peppler, a.a.O., S. 3, 12, 13. 
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Die Verformungsfunktion hat das asymptotische Verhalten K(¢) = 90 (£1), wie aus (13c) 
folgt. Die Spaltweite hat deshalb nach (14b) das asymptotische Verhalten h(¢) = 0(¢?). Der 
Integrand von (14c), den wir hier voriitbergehend mit f(¢) bezeichnen wollen, hat demnach das 
asymptotische Verhalten f(¢) =0(C~*). Wir setzten deshalb 


FC) perp ear (19a) 


und bestimmten die 4 Koeffizienten durch die Forderung, daf die abbrechende Reihe die 
bekannten Funktionswerte f(¢,) an den Stellen ¢€, = 3, 4, 5, 6 annimmt. Setzt man 
T(G)G =slc,) und? g(C,) =g,, soergibt sich 

13 ] } 


1 J i 1 5 
He) ~( 6 o2 9 84 9 85 2 Bs) iG Bene Ba g 85 2 bs) | (9b 
a 63 47 1 
| ( 3 83 t 784 — 2185+ > b4) r (20 gz— 45 ga 3085 Wes) aac 


Die Giite dieser Darstellung 148t sich kontrollieren, indem man die Reihe auch fiir die Punkte 
€ =2,5 und € = 3,5 berechnet und mit den unmittelbar gegebenen Werten des Integranden 
vergleicht. Es ergab sich dabei ein relative: Fehler von héchstens 0,005, der als sehr klein zu 
bezeichnen ist. 

Im Bereich € > 6 wurde der Integrand von (14c) durch die Darstellung (19b) ersetzt, 
streng integriert und zu dem mit der Simpsonregel fiir den Bereich 0 <¢€ <6 ermittelten 
Integrationsbeitrag hinzuaddiert. Nach 3 bis 4 Versuchen mit geeigneter Interpolation wurden 
so fiir die verschiedenen y- Werte die richtigen B-Werte gefunden. 

Fir u = 0 ergibt sich auf Grund der Pepplerschen Rechnungen B, = 0,73669, so daB man 
einen Anhalt fiir die GréBenordnung von B = B (yu) bei kleinen yu-Werten hat. 

Mit der Bestimmung von B (uw) hat man auch fast alle Arbeit zur Berechnung der bezogenen 
Druckfunktion (15a) geleistet, aus der sich durch Normierung p(¢) = p,(¢) + Ap(€) ergibt. 

Geht man mit nicht zu groBem Verformungsparameter (etwa 0< 4 <1) bei der ersten 
Iteration von der Verformungsfunktion Ky, (C) aus, so erhalt man als erste Naherung fiir die 
Druckdifferenz Ap (¢) eine Funktion, die, abgesehen von einem konstanten Faktor, recht gut 
durch die analytische Funktion 


ee ae 
Pi(S) [Tt (e 0.55) (20a) 
angenahert wird. Diese Tatsache ist deshalb bedeutungsvoll, weil sich die zugeordnete Ver- 
formungsfunktion K, (¢) mit tragbarem Arbeitsaufwand streng berechnen und vertafeln ] aBt. 
Auf Grund dieser Tatsache wurden fiinf weitere Funktionen ahnlicher Art eingefiihrt und zwar 


_ eae oe 
Palo) =F G05 “ 
Ge 1 ee Ne a 2 il Ss \S: 
piQ)= se pal) =, (20c) 
C2 1 ean (72 i er 3 
pal) =» pol) == (20) 
Ihre zugeordneten Verformungsfunktionen 
K,(6) = | pl) at (20e) 
0 


lassen sich leicht aus K, (¢) herleiten. 


Die Funktionen p,(¢) und K,(¢) wurden fiir die €-Werte (18) d-stellig vertafelt. Das 
erlaubt eine wesentliche Vereinfachung und Beschleunigung des Iterationsverfahrens in fo] gen- 


der Weise. 


Die durch Iteration gefundene Druckdifferenz wird durch die ,,elementaren Druckfunk- 
tionen*‘ angenadhert 


6 
Ap(6) ~ = a, p, (6) - (21a) 
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Wenn die 6 Koeffizienten a, bekannt sind, hat man sofort die zugeordnete Verformungsfunktion 


6 
AK(¢) = 2% K, (¢) (2ib) 
und kann mit 
= 6 
K(¢) = K)(¢) + 2% K, (¢) (21e) 


die neue Iterationsrechnung beginnen. 

Die durch die abbrechende Reihe (21a) erzielbare Genauigkeit der Anniherung von Ap(C) 
ist natiirlich begrenzt. Eine weitere Wiederholung des Iterationsverfahrens wird sinnlos, wenn 
der Fehler in Ap(C) bereits so klein geworden ist, da® er durch die Reihe (21a) nicht mehr 
vermindert werden kann. Die numerische Rechnung zeigt, daf fiir 0 <p <1,6 mit diesem 
Verfahren nach 3 bis 4 Iterationsschritten ein Ergebnis erzielt wird, das allen praktisch erfor- 
derlichen Genauigkeitsanspriichen geniigt. Fir 1,6< 4 <2 waren wesentlich mehr elemen- 
tare Druckfunktionen nétig gewesen. Der dazu erforderliche Arbeitsaufwand zur Berechnung 
und Vertafelung der zusatzlichen Funktionen p,(¢) und K,(¢) erschien zu groB, so daB ein 
neues Verfahren gesucht und gefunden wurde, welches jede beliebige Genauigkeit zu erreichen 
gestattet. 


7. Die numerische Berechnung der Verformungsfunktion K (¢). Zunachst wird kurz an- 
gegeben, wie die Integrale K,(C),...,Kg(¢) berechnet wurden. Wegen 


po(0) =0 und — poco) =0 


erhalt man durch partielle Integration und anschlieBende Differentiation 


i) By oes pe ee ff 1 [9.67109 —7,88994 
aa o(5) =[ rane =| reer, la agua TE earn (22) 
0 0 
Mit den Abkiirzungen 
+ co 
i 1 
oe -| CoS «i tet 
a = 0,5 -0,95026 = 0,47513 , u = 2 a(¢ — 0,5), v =2 a(é — 0,5) (23b) 
gilt also 
< K,(6) = 9.67109 H,(w, a) — 7,88994 Hy(u, a) . (23c) 
Mit der Rekursionsfor mel 
1 1 ut 
Hy41 (us @) =~ M, (u, «) — ibe eae (24a) 
lassen sich miihelos alle H,, mit m = 2, 3, 4,... berechnen aus der Grundbeziehung 
u I 1 Jl + a? ; 24b 
HH, (u, a) = ce | are tg a | Teme lead (24b) 


Gebraucht wird nun noch der Wert von K, an der Stelle € = 0,5. Man findet mit der Be- 
ziehung (16) 


oo 


= are tgv 
K,(0.5) => | cee 1 i nS v |i : 
os ear +e 4 aretga)—F —aretga - = In |a|—a 4 aa ae oe 
K,(0,5) = 2,48687 . 
Damit erhalt man schlieBlich ; 
Ky(6) = Ko(0.5) + | £ Ky(é) dv, (25b) 
é 


13 
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wobei u und v durch (23b) erklart sind und der Integrand durch (23c), (24a), (24b) numerisch 
berechnet werden kann. Mit Hilfe numerischer Integration (z. B. Simpsonregel) 1aBt sich mit 
(25b) die Funktion K, (¢) beliebig genau berechnen. Spater wird noch ein Verfahren zur 
direkten Berechnung von K, (¢) angegeben, das uns als Kontrolle diente. : 

Mit Hilfe der Bezichungen (24a) und (24b) kénnen wir auch unmittelbar die Funktionen 
K, (¢) bis K, (¢) berechnen, Durch eine elementare Zerlegung von p, (¢) erhalt man 


oe 1 Sk iymeeees oe 
ENG, milberert ee rapa | Toa + oy ed) 
—0,5 
mit den Abkiirzungen 
x=&—0,5, Bat — 0,3. (26b) 
Verwenden wir wieder die durch (23a), (24a), (24b) erklarten Funktionen, so ergibt sich 
K,(¢) = H,(z, 0,5) — 2 H(z, 0,5) + 2 Hz 0,5) (26c) 


Alle itbrigen drei Funktionen kann man auf K,(¢) zuriickfiihren : 


K,(0) =(1—2) K,(0)— f (6) a = (1 — ¢) K,(¢) — 0,0944573 , (26d) 
0 
K.0) => a K© SOA ae) as] 
te as (26) 
= TEE aap | KO) — gop — 0.0472596 (¢ — 0,5)), 
K,(t) = (1 — ¢) K,(¢) — J ps(€) d& = (1 —¢) K,(¢) — 0,0472596 . (26f) 
0 
Mit entsprechenden Uberlegungen findet man 
Ks(0) =a gpa lKs(C) — 0,0349870 ¢ + 0,0089602] , (26g) 
Ki) 2 = ey RO) —0.0840870) (26h) 


Da bei diesem RekursionsprozeB zweifellos an Genauigkeit eingebiiBt wird, wurde die 
numerische Rechnung mit einer fiir unser Naherungsverfahren reichlich hohen Stellenzahl, 
namlich mit 5 geltenden Ziffern begonnen. 

Weil mit der Reihe (21a), wie schon erwahnt, nur im Bereich 0 <u <1,6 ausreichende 
Genauigkeit erzielt wurde und auch Hinzunahme weiterer Reihenglieder ahnlicher Art diesen 
Mangel nicht ausreichend behob, mufte im Bereich 1,6 < w~ <2 ein anderes Verfahren an- 
gewandt werden. Es erfordert zwar einen gréBeren Arbeitsaufwand, ]a8t aber durch Erhéhung 
der Anzahl] der Reihenglieder jede gewiinschte Genauigkeit zu. 

Wir betrachten das Integral 


K@) = | PO dE. (27a) 
0 


Hier gilt entsprechend unseren Abmachungen und Voraussetzungen 


p(é)=O(E*) fir &-> ©, ] 


2 
p(é)=0(@) fir ESO. | Oy 
Mit der Transformation 
ee ae 
C 1 + cosa AW ore acct ae (28a) 


und der Bezeichnung 


P(B) = p(E (6) (28b) 


nimmt (27a) die Form an 


re ‘ P() sin 
seh mee PU) (1 + cos B) (cos B — cos a) ap . (28c) 
0 
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Das hier auftretende Integral hat, als Cauchyscher Hauptwert verstanden, einen Sinn fiir 
beliebige Werte 0 <a <a, da P(f) als stetig und glatt vorausgesetzt ist und fiir b—> 0 und 
f > m wegen (27b) geniigend stark gegen Null strebt. 

Das Integral in (28c) 148t sich mit einer aus der Tragfliigeltheorie bekannten Formel 


rey eee A tars fit on 0,1, 2.64 


cos 6h — cos & sina 
0 


und der daraus folgenden Beziehung 


4 
eae peecer: 


cos fh — cosa 


dp =—acosna fiir Resa We eos. Ke 
0 


leicht berechnen, wenn man die Fourierdarstellung wahlt 


P(8) N . 
ape sin np. (28d) 
n=1 
Dann erhalt man namlich unmittelbar das gesuchte Integral in der Form 
K(¢) =a(1+ cosa) Dic, cosna. (28e) 
n=1 


Praktisch wird man so vorgehen, daf man K, (¢) fiir sich berechnet, damit man nach jedem 
Iterationsschritt nur noch die Druckdifferenz A p (¢) mit dem zugeordneten Integral AK (C) 
zu berechnen braucht. 

Die numerische Rechnung zeigt, da} man zur geniigend genauen Darstellung von AK (¢) 
leider sehr viel Reihenglieder benétigt und zwar wenigstens 23 fiir ~ < 1,7. Darin besteht der 
Nachteil dieses Verfahrens gegeniiber dem vorher besprochenen, bei dem man fiir 0 < yw < 1,6 
mit 4 bis 6 Reihengliedern auskam. Der Vorteil des neuen Verfahrens ist der, da® die Erhéhung 
der Anzahl] der Reihenglieder keine prinzipiellen Schwierigkeiten macht. AuBerdem kann man 
die umfangreichen Erfahrungen iiber Fourierreihen ausnutzen. 

Mit der nach (28) bestimmten K-Funktion wurde der Fall 4 = 1,7 behandelt. Der Fall 
po = 1,8 wurde ohne Erfolg in Angriff genommen, denn das Iterationsverfahren fihrte trotz 
zahlreicher Versuche nicht zur Konvergenz. Die Ursache dafiir liegt nicht in der Berechnung 
von K (C), sondern in der Tatsache, dafBs bei so starken relativen Verformungen sich die Spalt- 
weite als kleine Differenz groBer Zahlen ergibt. Fiir ~—>2 treten auBerdem an den Stellen 
€ = O0und€ = 2 Ecken in der Funktion K (¢) auf, was die Schwierigkeiten noch erhéht. Leider 
sind nun aber y-Werte nahe unterhalb 2, das heifBt Verformungszustande in der Nahe des 
Hertzschen Grenzfalles, fiir den Maschinenbau besonders interessant. Deshalb soll spater ver- 
sucht werden, wenigstens noch den Fall 4 = 1,8 numerisch zu bewaltigen. Der dazu erforder- 
liche Arbeitsaufwand diirfte betrachtlich sein. 


8. Die numerischen Ergebnisse des Iterationsverfahrens. Tabelle 1 gibt die numerischen 
Werte der elementaren Druckfunktionen p,(¢) wieder, von denen einige in Abb. 5, 6, 7 graphisch 
dargestellt sind. Tabelle 2 enthalt die elementaren Verformungsfunktionen K,(¢), von denen 
einige in den Abb. 8, 9, 10 graphisch dargestellt sind. Mit Hilfe dieser Elementarfunktionen 
wurden die Rechnungen fiir die Verformungsparameter w = 0,5 1,0 1,25 1,5 1,6 durchgefiihrt. 
Bei all diesen Parameterwerten wird bereits nach dem dritten Iterationsschritt Konvergenz im 
Rahmen der fiir technische Probleme erforderlichen Genauigkeit erzielt. Abb. 11 zeigt fiir 
fu = 1,6 die nach dem zweiten und dritten Iterationsschritt erhaltene Druckdifferenz Ap(¢), 
um die Gitte der Genauigkeit des Endergebnisses zu zeigen. Abb. 12 zeigt die durch die zweite 
Iteration ermittelte Ap-Kurve und ihre durch » a, p,(¢) angenaherte Ersatzfunktion fiir die 
dritte Iterationsrechnung. Der Vergleich beider Kui ven gibt ebenfalls einen Eindruck von der 
GroBenordnung der erzielbaren Genauigkeit. 

Bei der ersten Iteration wurde als Grundlage die Druckdifferenz des nachstbenachbarten 
p-Wertes gewahlt, fiir den das Iterationsverfahren bereits vollstandig durchgefiihrt war. 
AuBerdem wurde bei den ersten Iterationen zur Berechnung der Verformungsfunktion A K (¢) 
nach (21b) die Druckdifferenz Ap(¢) nur durch die ersten 3 bis 4 Elementarfunktionen p, (C) 
angendhert. Lediglich bei der letzten Iteration wurden — falls nétig — alle 6 Elementar- 


13* 
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Tabelle 1. Elementare Druckfunktionen p, (¢) fiir vy = 0-- > 6. 


| 
| 


| 
Po (6) Pils) p2(6) Ps (6) ps (6) ps () Po (6) 


Tay 


0 0 0 0 0 0 0 
0,019 73 0,004 474 0,004026 | 0,003 857 0,003 471 | 0,003 325 0,002 992 
0,082 80 0,018 136 0,014509 | 0,016 639 0,013 311 | 0,015 265 0,012 212 
0,190 15 0,037 697 0,026 388) 0,036 247 0,025 372 | 0,034 853 0,024 396 
0,334 48 0,055 352 0,033 212 | 0,054 804 0,032 882 | 0,054 261 0,032 556 
0,500 00 0,062 500 0,031 250] 0,062 500 0,031 250 | 0,062 500 0,031 250 
0,665 53 0,055 352 0,022141 |] 0,054 804 0,021 922 | 0,054 261 0,021 705 
0,809 86 0,037 697 0,011 309} 0,036 247 0,010 874 | 0,034 853 0,010 456 
0,917 21 0,018 136 0,003 627 | 0,016 639 0,003 328 | 0,015 265 0,003 053 
0,980 28 0,004 474 0,000 447) 0,003857 0,000 386 | 0,003 325 0,000 333 
1,000 00 0 0 0 0 0 0 
0,983 19 0,003 537 | —0,000 354 | 0,002601 | —0,000 260 | 0,001912 | —0,000 191 
0,939 30 0,011686 | —0,002 337] 0,007843 | —0,001569 | 0,005 264 | —0,001 053 
0,877 94 0,021026 | —0,006 308} 0,012821 | —0,003 846 | 0,007818 | —0,002 345 
0,807 30 0,029 219 | —0,011688 | 0,016143 | —0,006 457 | 0,008919 | —0,003 567 
0,733 59 0,035 156 | —0,017578 | 0,017578 | —0,008 789] 0,008789 | —0,004 395 
0,661 08 0,038 634 | —0,023181] 0,017481 | —0,010489 | 0,007910 | —0,004 746 
0,592 40 0,039 952 | —0,027966 | 0,016374 | —0,011 462) 0,006711 | —0,004 698 
0,529 02 0,039602 | —0,031 682} 0,014722 | —0,011778| 0,005473 | —0,004378 
0,471 54 0,038091 | —0,034282 |) 0,012 869 | —0,011582) 0,004348 | —0,003 913 
0,420 02 0,035 853 | —0,035 853 | 0,011032 | —0,011032 | 0,003394 | —0,003 394 
0,333 74 0,030 438 | —0,036525| 0,007825 | —0,009390 | 0,002012 | —0,002 414 
0,239 01 0,022500 | —0,033 750] 0,004500 | —0,006 750) 0,000900 | —0,001 350 
0,143 13 0,013030 | —0,026060; 0,001 797 | —0,003594} 0,000 248 | —0,000 496 
0,090 68 0,007656 | —0,019141 | 0,000 766 | —0,001914]| 0,000077 | —0,000191 


s 


SOSONOUNSDOONIANEBWNHHSOBNIAUNBWHH 


’ 


| 0,06047 | 0,004672 | —0,014016) 0,000353 | —0,001058}| 0,000027 | —0,000 080 
, 0,030 34 | 0,001962 | —0,007 847} 0,000092 | —0,000 369} 0,000004 | —0,000 017 
; 0,016 23 0,000 944 | —0,004719 | 0,000030 | —0,000151 | 0,000001 | —0,000005 
Peo 

40 

08 

06 - 

OF 

| 

| ¢ 
0 7 a 3 4 
Abb. 5. Elementare Druckfunktion p,(é) = normierte Druckfunktion des unverformten Rollenlagers. 
‘a 
0,08 
fe 
GO8 T 


(ar (a c)2 Abb. 7. Elementare Druckfunktionen 
1 j= — ed tj) = —_ re Se er(1 —_ fe 
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Tabelle 2. Verformungsfunktionen K,(C) fir y =0---6. 
ee ee CANES 8 

é Ky (6) K, () K, (¢) K,(¢) K,(¢) K; (6) K, (6) 

0 1,497 09 0,100 69 0,006 23 0,078 12 0.030 86 0,069 66 0,034 6 

0.05 1,623 92 0,116 04 0,015 78 0,092 01 0.040 15 0,082 51 0,043 40 
(eal 1,768 54 0,133 01 0,025 25 0,107 97 0,049 91 0,097 78 0,053 02 
0,15 1,919 59 0,148 18 0,031 50 0,122 99 0,057 28 0,112 87 0,060 95 
O22 2,068 51 0,158 54 0,032 38 0,133 99 0,059 93 0,124 72 0,064 79 
0,25 2,206 97 0,161 46 0,026 64 0,137 98 0,056 23 0,130 06 0,062 56 
0,3 2,326 89 0,155 02 0,014 06 0,132 50 0,045 49 0,125 92 0,053 16 
0,35 2,420 72 0,138 47 | —0,004 45 0,116 28 0,028 32 0,110 50 0,036 84 
0,4 2,481 81 0,112 47 | —0,026 98 0,089 63 0,006 52 0,083 76 0,015 27 
0.45 2,504 99 0,079 18 | —0,050 94 0,054 74 | —0,017 16 0,047 84 | —0,008 68 
0,5 2,486 87 0,041 96 | —0,073 48 0,015 30 | —0,039 61 0,006 77 | —0,031 60 
0,55 2,426 14 0,004 90 | —0,092 25 | —0,02407 | —0,05809 | —0,034 26 | —0,050 40 
0,6 2,323 70 | —0,02795 | —0,105 64 | —0,058 75 | —0,070 76 | —0,070 08 | —0,063 02 
0,65 2,182 47 | —0,053 17 | —0,113 07 | —0,085 01 | —0,07701 | —0,096 62 | —0,068 80 
0,7 2,007 17 | —0,06860 | —0,115 04 | —0,100 69 | —0,077 47 | —0,111 75 | —0,068 51 
0,75 1,803 84 | —0,073 52 | —0,112 84 | —0,105 41 | —0,073 61 | —0,115 47 | —0,063 86 
0.8 1,579 35 | —0,069 39 | —0,108 34 | —0,10114 | —0,06749 | —0,110 25 | —0,057 04 
0.85 1,340 88 | —0,055 69 | —0,102 81 | —0,08811 | —0,060 48 | —0,09701 | —0,049 54 
0,9 1,095 42 | —0,03719 | —0,09818 | —0,071 34 | —0,054 39 | —0,080 92 | —0,043 08 
0,95 0,849 39 | —0,015 81 | —0,095 25 | —0,05301 | —0,049 91 | —0,06427 | —0,038 20 
1,0 0,608 38 0,005 99 | —0,094 46 | —0,035 45 | —0,047 26 | —0,04918 | —0,034 99 
I] 0,158 70 0,043 27 | —0,098 79 | —0,008 64 | —0,046 40 | —0,02807 | —0,032 18 
iy? —0,228 96 0,065 40 | —0,107 54 0,003 79 | —0,048 02 | —0,019 62 | —0,031 06 
3 —0,545 49 0,071 13 | —0,115 80 0,00405 | —0,048 48 | —0,019 80 | —0,029 05 
1,4 —0,791 24 0,063 95 | —0,12004 | —0,002 90 | —0,04610 | —0,023 71 | —0,025 50 
158 —0;973 15 0,048 81 | —0,118 87 | —0,012 56 | —0,040 98 | —0,028 04 | —0,020 97 
1,6 —1,100 93 0,030 13 | —0,11254 | —0,021 95 | —0,03409 | —0,031 21 | —0,016 26 
1.7 | —1,184 79 0,011 12 | —0,102 24 | —0,02961 | —0,026 53 | —0,032 84 | —0,012 00 
1,8 —1.234 24 | —0,006 31 | —0,089 41 | —0,035 10 | —0,01918 | —0,03313 | —0,008 48 
1,9 —1,257 41 | —0,021 24 | —0,075 34 | —0,038 54 | —0,01257 | —0,032 45 | —0,005 78 
2,0 —1,261 03 | —0,033 37 | —0,06109 | —0,040 29 | —0,006 97 | —0,03117 | —0,003 82 
Zell —1,250 45 | —0,042 80 | —0,047 38 | —0,040 75 | —0,002 43 | —0,02957 | —0,002 46 
252 —1,229 78 | —0,049 83 | —0,03466 | —0,040 32 0,001 12 | —0,027 85 | —0,001 57 
7a —1,202 21 | —0,054 84 | —0,023 16 | —0,039 29 0,003 81 | —0,026 13 | —0,001 02 
2,4 —1,170 13 | —0,05818 | —0,01300 | —0,037 88 0,005 78 | —0,02449 | —0,000 70 
J) —1,135 30 | —0,060 20 | —0,00416 | —0,036 28 0,007 16 | —0,022 96 | —0,000 55 
3.0 —0,954 55 | —0,059 00 0,023 54 | —0,028 11 0,008 96 | —0,01712 | —0,000 75 
a5 —0,798 46 | —0,051 22 0,033 59 | —0,021 97 0,007 67 | —0,013 55 | —0,001 11 
4,0 —0,675 54 | —0,043 36 0,035 61 | —0,017 77 0,006 05 | —0,011 23 | —0,001 30 
5,0 —0,579 82 | —0.031 60 0,031 96 | —0,012 75 0,003 74 | —0,008 41 | —0,001 35 
6,0 —0,504 85 | —0,024 22 0,026 64 | —0,009 95 0,002 49 | —0,006 75 | —0,001 24 


’ 


K 


ou “ 
Abb. 8. Elementare Verformungsfunktion K,(¢) =4h S © dé. 
0 


funktionen im Reihenausdruck (21a) herangezogen. In diesem Fall wurden die sechs Koeffi- 
zienten a, so bestimmt, daB die Funktionen in 4 Punkten tibereinstimmten und daB die beiden 


Stellen ihrer Maxima sich deckten, 
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Abb. 9. Elementare Verformungsfunktionen K,(¢) und K,(¢) . Abb. 10. Elementare Verformungsfunktionen K,(¢) und K,(¢). 


G7 


Abb. 11. Durch Iteration gefundene Druckdifferenz. — Ap, (€) (2. Iteration); - - - 4p; (C) (3. Iteration); yw = 1,6. 


07 


-O7 


6 
Abb. 12, — Durch 2. Iteration gefundene Druckdifferenz Ap, (6); --- Ersatzfunktion 4p(¢) = Dd’ a Pp, (2) 
als Grundlage fiir die 3. Iteration; = 1,6. ges 


Aye 


Abb. 13. — Durch 4. Iteration gefundene Druckdifferenz Ap, (€) 3 


23 
--- Ersatzfunktion 4p [¢ (A)] = (1 + cos f) Dc, sinn fp; p = 1,7. 
n=1 
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Fir uw = 1,7 wurde K (€) fiir jede folgende Iteration mit Hilfe der Beziehungen (28) berech- 
net. Mit diesem y-Wert treten beim Iterieren zum ersten Mal Schwierigkeiten auf. Die 
Ap-Kurven erster und zweiter Naherung weisen groBe Unterschiede auf. Die Ap-Kurve 
dritter Naherung liegt wieder nahe bei der erster Naherung, so daB das iibliche Verfahren als 
vierte Naherung wieder ungefahr das Ergebnis der zweiten ergeben hatte. Um diesem Pendel- 
vorgang auszuweichen, wurde die Druckdifferenz erster und zweiter Naherung gemittelt, so daB 
bei der vierten Naherung die Druckdifferenz 


= {4p (Oh + [4P(O))} 


zugrunde lag. Das so erhaltene [Ap], ist bereits eine sehr gute Naherungslisung, denn das 
damit berechnete [Ap]; unterscheidet sich im Rahmen der Zeichengenauigkeit kaum von [Ap],. 
Abb. 13 zeigt die Ap-Kurve vierter Naherung und ihre Approximation durch 


(1 + Sao sin nf 
n=1 


fiir die fiinfte Iterationsrechnung. 

Fir « = 1,8 wird das bei ~ = 1,7 bereits beobachtete Pendeln der Iterationslésungen 
starker. Auch Mittelwertbildungen zwischen den Ergebnissen zweier aufeinander folgenden 
[terationen fiithrt nicht zur befriedigenden Konvergenz. Trotz langwieriger Versuche verschiede- 
ner Art ist noch keine geniigend genaue Lésung bekannt. Die Zusammenstellung der Ergeb- 
nisse zeigt aber spater, daB die Ubersicht iitber den Parameterbereich 0 <u <1,7 bereits 
ziemlich klare Schliisse dariiber erlaubt, was in dem nicht untersuchten Restbereich leat Se 
eintreten wird. 

Alle jetzt folgenden Abbildungen und Erlaute- I -A 
rungen beziehen sich auf den jeweilsletztenIterations- 2 + 
schritt der Rechnungen. ame 

Abb. 14 zeigt die Parameter A und 1/[(2 — yw) B?] 
iiber w aufgetragen. Die erste Kurve ist fast eine 
Gerade, die zweite weicht nicht viel von einer Kon- 
stanten ab. Diese Tatsache erleichtert die Vorher- 
bestimmung yon A und B bei Iterationsrechnungen 
fiir neue y- Werte. 

Abb. 15 zeigt die Druckdifferenzen Ap und Abb.16 Wait 
die normierten Druckfunktionen P= py+ Ap iiber 
€ aufgetragen mit uals Parameter. Man erkennt, wie 
mit ~—>2 die Druckkurven der Form der Hertzschen 
Halbellipse zustreben. 

Abb. 17 gibt die Verformungsfunktionen iiber ¢ ji ; 3 te 
aufgetragen wieder mit yw als Parameter. Hier er- : 
kennt man besonders deutlich, wie mit ~—>2 die ABN 1S ser arartet ol Shey ae ar nee 2 
Kurven im Bereich 0<€ <2 der Geraden 2 (1—€) aufgetragen. 
des Hertzschen Grenzfalles zustreben. 

Abb. 18 zeigt noch fir ~ = 0 und 1,0 und 1,7 die bezogene Spaltweite iiber € auf- 
getragen. Die Kurve fiir w = 1,7 1aBt wieder deutlich erkennen, wie im Bereich 0 <¢€ <2 
der Rollenrand durch den Einflu® der Schmiermittelkrafte sich einer horizontalen Geraden 
angeglichen hat. 


9. Ubertragung der Ergebnisse auf das wirkliche Rollenlager. a) Effektiver Rol- 
lenhalbmesserundeffektiver Elastizitatsmodul. In unseren Unter- 
suchungen wurde bisher der Untergrund als starr und eben angenommen. Was hier als Unter- 
grund bezeichnet wird, ist aber in Wirklichkeit die Lagerschale oder die Lagerwelle, die beide 
im unverformten Zustand gleichmabig gekriimmt sind. 

Bei der Herleitung der Schmiermittelkrafte bleiben alle Beziehungen (2) bis (4) mit dem- 
selben Naherungsgrad giiltig, wenn man an Stelle des ebenen Untergrundes einen gekriimmten 
wablt, sofern nur im Bereich, wo der Schmiermitteldruck nicht nahe bei Null ist, die Kriim- 
mungsradien der Randflachen groB sind in Bezug auf die Spaltweite H (x). 

Der Rollenradius sei R,, der Wellenradius R,,, der innere Schalenradius R,, dann sind alle 
bisherigen Ergebnisse giiltig, wenn man in den Formeln fiir R folgenden effektiven Wert 
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einsetzt: 


Rr Rw rue Ry Rs 29 
R= St Re bzw. Rak . (29) 


Der linke Ausdruck gilt fiir den Spalt zwischen Rolle und Welle, der rechte fiir den Spalt 
zwischen Rolle und Schale. Die Beziehungen (29) gibt auch W. Peppler * fiir das nichtverform- 
bare Rollenlager an. 

Schale und Welle haben nicht nur gekriimmte Oberflachen, sondern sind auch elastisch. 
Das unseren Untersuchungen zugrunde gelegte Verformungsgesetz in der differenzierten Form 


Abb. 16. Durch letzte Iteration gefundene normierte Druckfunktion p(¢) = po(¢) +4 p(¢); uw als Parameter. 


(13c) ist wegen der in (10) vorgenommenen und dort begriindeten Linearisierung vom Rollen- 
radius ganzlich unabhangig. Es gilt deshalb, wie man auch direkt zeigen kann, fiir beliebige 
positive und negative Rollenradien. Das bedeutet, da es auch fiir die Oberflache von zylin- 
drischen Bohrungen in einem unendlich ausgedehnten Kérper und fiir den Grenzfall der ebenen 
Randflache eines unendlichen Halbkérpers gilt. Naherungsweise gilt das Verformungsgesetz 
aber auch noch fiir die beiden Oberflachen von Kreisréhren und zwar umso genauer, je dicker 
die Réhrenwand bezogen auf den Radius der Bohrung ist. Die Bohrung braucht nicht einmal 
zentrisch zum AuBenzylinder sein, wenn iiberall die Wandstarke gro8 gegeniiber dem Bohrungs- 
radius ist. 

Diese Tatsache erlaubt eine wichtige Folgerung. Es seien E,, E,,, E, die Elastizitatsmoduln 
von Rolle, Welle und Schale. Dann bleiben alle bisher mitgeteilten Ergebnisse giiltig, wenn 


19 Hy Elo Dan Lioiereaall ()). 
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man fiir den in den Formeln auf- 
tretenden E-Modul folgende 
effektiven Werte einsetzt: 


E; Ep 7 
ees se | 
bzw. ( (30) 
_ EE Es 
Ey -b Es j 


Der linke Ausdruck bezieht sich 
auf den Spalt zwischen Rolle 
und Welle, der rechte auf den 
Spalt zwischen Rolle und Schale. 
Fir den Spalt zwischen Rolle 
und Schale ist die Ubertragung 
der Ergebnisse nur naherungs- 
weise erlaubt. Wahrscheinlich 
ist die Naherung aber nicht nur 
dann gut, wenn die Lagerschale 
dick zum Innenradiusist,sondern 
auch dann noch, wenn die Scha- 
lendicke wenigstens von der- 
selben GréBe wie der Rollen- 
radius ist. In unseren Herleitun- 
gen und Formeln wird namlich 
nur ein Spaltbereich betrachtet, 
der schmal ist in Bezug auf den 
Rollenradius. Wenn dieser Spalt- 
bereich auch schmal ist in Bezug 
auf die Schalendicke, so wird 
unser Verformungsgesetz in die- 
sem schmalen Bereich auch fiir 
die Lagerschale naherungsweise 


ho 


EA 41__ 


Abb. 17. Durch letzte Iteration gefundene Verformungsfunktion 
co 
pO cee (2 Par 

Kmax Kmax, c 


0 


pw als Parameter. 
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Abb. 18. Durch letzte Iteration gefundene bezogene Spaltweite h (é)/h (0); als Parameter. 
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giiltig sein. Eine genaue Untersuchung dieses elastizitatstheoretischen Problems diirfte 
schwierig sein. 


b) Berechnung des StreckungsmaB8es 4. Alle berechneten Druck- und 
Verformungsfunktionen beziehen sich auf eine normierte Ortskoordinate C derart, dah P (0) 
= p’(1) =0 ist. Wir wollen nun von der dimensionslosen Koordinate ¢ auf die wirkliche 
Lingenkoordinate x iibergehen und miissen dazu A berechnen. 

Zur Herleitung der asymptotischen Darstellung (19) wurde bereits erwahnt, daB K (C) fiir 
eroBe C-Werte wie k/¢ gegen Null strebt. Daraus folgt aus (13a) fiir groBbe ¢-Werte, das heibt 


fir ¢ SJ 


h''(C\ eo B hi(0)- (31a) 
Mit x =AC + », und der Abkiirzung h(0)= A, erhalt man daraus die asymptotische Beziehung 
Cplal ON Tashi s 
a fiir KS>A+ %, (31b) 
wobei H(x)=h[€(x)]ist. Andererseits gilt im Bereich 
L+ m<x<R (32a) 
die bekannte Naherungsformel 
ae (32b) 
dx” R 


Die Bedingung x < R ist bei allen unseren Untersuchungen implizit vorausgesetzt worden. 
Wir haben namlich immer angenommen, daB8 sich im Bereich des Druckberges und einer im 
Vergleich zu (A+ x,) groBen Umgebung der unverformte Rollenrand durch eine Parabel 
anndhern la8t. Nur mit dieser Voraussetzung hat der mathematische Grenziibergang ¢, —> 00 
(sehr starke Schmierung) einen Sinn. Dann gibt es aber einen Bereich (32a), fiir den (31b) und 
(32b) in guter Naherung giiltig sind, und es folgt 
Ae = 2 Bh he (33) 

Diese Beziehung gilt, solange unsere hydrodynamische Grundbeziehung (2c) giiltig ist, umso 
genauer, je kleiner hy ist. Auch W. Peppler macht bei seinen Untersuchungen fiir das unver- 
formbare Rollenlager implizit dieselben Voraussetzungen. 

In (33) ist der Rollenradius R vorgegeben, B = B (uw) durch die numerischen Rechnungen 
ermittelt, hy aber vorlaufig noch unbekannt. Zur Berechnung von hy miissen wir auch noch die 
Gesamt-Rollenlast P je Einheitsbreite vorgeben. Es ist 


P = [ q(x) dx =A f pe) ae. (34a) 


Die wirkliche Druckfunktion q(x) bzw. p(¢) kennen wir nicht, sondern nur die bezogene 
Druckfunktion (13b). Wir kénnen also das Integral 


0= [lenin] epi (34b) 


unmittelbar aus den Ergebnissen unserer Iterationsrechnung ermitteln. Die Werte Q = Q (12) 
sind dimensionslos und hangen nur von yu ab. 

Berechnen wir hy aus (34b) und setzen diesen Ausdruck in (33) ein, so erhalten wir schlieBlich 
die gesuchte Beziehung 


Q(u) 
A 2A Bu) ’ (35) 
In (35) sind die Werte R, 7, V, Pals technische Daten vorgegeben, wahrend B und Q durch die 


Iterationsrechnung geliefert werden, wenn der Verformungsparameter ut bekannt ist. Formel 
(35) ist also nur dann auswertbar, wenn wir noch jz durch die technischen Daten des Rollen- 
lagers ausdriicken kénnen., 


c) Bestimmung des Verformungsparameters faustechnischen 
Daten des Rollenlagers. In Formel (13a) wird die erste Ableitung der normalen 
Randverschiebung, die voritbergehend mit o'(¢) bezeichnet wird, durch den Ausdruck be- 
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schrieben 
eC) =n Bie mit K(t) = | pe. dé. (36a) 
= 
; i 
Auf Grund von (11a) folgt durch Differentiation und Ubergang auf die Koordinaten (12a) 
ey) 20 =) fre) 
e'(2) = E [eee #: (36b) 


i 
In unseren numerischen Rechnungen wurde (36a) fiir die normierts Druckfunktion p (¢) 


berechnet. Numerisch liegen namlich die Funktionen p(¢) und K (¢) vor. 
Die beiden Beziehungen (36a) und (36b) haben eine weitere fiir uns niitzliche Beziehung zur 


Folge: 


2(1—v)A Hie 
aE PC) == Be). 
Durch Integration von € = 0 bis € = o erhalt man daraus 
tH wu Bh, 
(1— Des rea {PO p()d 


Hier stért der zunachst noch unbekannte Term hy. Ersetzen wir in (35) die linke Seite durch 
die rechte Seite von (33), so finden wir 


hy = 12 Ry VB(w) SY. (37) 


Setzt man das in die vorhergehende Beziehung ein und bringt alle nur von ~ abhangigen 
Glieder auf die eine und alle nicht von u abhangigen auf die andere Seite, so erhalt man 

Gh Ge xs 2 ped 

RAVE =e eB) Ou) | BE) a, (38) 


max i 
0 


Rechts steht ein Ausdruck, der eine eindeutige Funktion von y darstellt und durch die Erebnisse 
der Iterationsrechnungen bekannt ist. Wenn also die linke Seite von (38) durch die Betriebs- 
daten des Rollenlagers vorgegeben ist, kann der zugehérige w- Wert unmittelbar bestimmt wer- 
den. Die Werte Rund E sind gema8 (29) und (30) einzusetzen. 

d) Bestimmung der Mindestspaltweite und des gréften 
Schmiermitteldruckes. Durch die numerischen Rechnungen ist die bezogene 
Spaltweite und somit auch der Wert Amin/hy bekannt. Deshalb folgt aus (37) unmittelbar 

hmin = 12 Ry VB (yu) 20) Pi (39) 
‘0 

Fiir die Technologie der Schmierung ist neben der Mindestdicke h;, des Schmierfilms auch 

noch der maximale Schmiermitteldruck Png, von Bedeutung. Seine Berechnung ist sehr ein- 


fach. Aus 


P=A | p(c) dt =2 Prax | B(G) a 


folgt 


ws (40a) 


Pmax = ao 


Af p(s) de 


Die normierte Druckfunktion P (¢) und somit auch ihr von 0 bis co genommenes Integral sind 
bekannt. Das StreckungsmaB A ist durch (35) gegeben. Setzt man den aus (35) sich ergebenden 
Ausdruck ein, so erhalt man: 


Pmax mae: Le ; Ss . (40b) 
4 a) f plo) a 
0 
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10. ZweckmaBige Fassung der Ubertragungsbeziehungen und allgemeine Folgerungen. Wir 
kennen jetzt alle Beziehungen, die zur Ubertragung unserer numerischen Ergebnisse auf ein 
wirkliches Rollenlager erforderlich sind. Die Umrechnung 1a8t sich aber wesentlich verein- 
fachen, wenn man einige von 1 abhangige dimensionslose Funktionen einfithrt. . 

Zur technischen Auswertung der vorliegenden Ergebnisse mu8 man als erstes den zu den 
Betriebsdaten des Rollenlagers gehérigen u-Wert kennen. Seine Bestimmung ist am einfach- 
sten, wenn man (38) in folgender Form schreibt: 


lai. 
PY ape = 7) (41a) 


mit 


6 a uw Q(u) SO dé 
0 1 ng (41b) 


p(w) = B(w) fe 


7 (UH) (eee 


iiber yw dargestellt, da sich herausstellte, daB sie in dem untersuchten Bereich 0 < w <1,7 fast 
konstant ist. Sie weicht nur um héchstens + 3,0% vom Mittelwert 2,235 ab. Es gilt also in 
guter Naherung 


In Abb. 19 ist die Funktion 


(wu) 2 235 )/-* 42a 
p(w) ~ 2, fo (42a) a 
010 

| f- 
ee ee 

0 et) a 

[=u Abb Berech et 
Abb, 19. o Berechnete Werte von  (y) | —— - 20. © Berechnete Werte von (2 nu) y(n) 
. be eal j 
— Mittelwert 2,235. Se irsatzfaukction “st , 


Das erlaubt eine numerische Berechnung des Parameters w mit einer Genauigkeit, die allen 
technischen Anspriichen geniigt. Man erhalt, wenn man (42a) in (41a) einfiihrt, 


2 

> Ry7V E 
1+ (1—2?) P? 
Diese Naherungsbeziehung erleichtert die Auswertung unserer numerischen Ergebnisse in 
auBerdordentlicher Weise. Man muB8 es deshalb als Gliicksfall ansehen, daB der Verformungs- 
parameter win der Form (13a) definiert wurde. Zu Beginn der Untersuchung! wurde eine 
andere Definition des Verformungsparameters gewahlt, der nicht das einfache Gesetz (42b) 
ergeben hatte. Diese alte Definition wurde fallen gelassen, weil sie zu mehreren Unannehm- 
lichkeiten anderer Art fiihrte. 

Nachdem der w-Wert bekannt ist, kann das StreckungsmaB 2 bestimmt werden. Wir 
schreiben zu diesem Zweck (35) in folgender Form: 


2 = RV ww (43a) 


bX (42b) 


mit 
y(u) = B(u) $24 Q(u) . (43b) 


Kine fiir manche Fragestellungen zweckmaBigere Form der Beziehung erhilt man, wenn man 
hier P mit Hilfe von (41a) eliminiert: 


g- VER vw) 
) E Vp(u) 


1 J. Dorr, Z, angew. Math. Mech, 32 (1952) 268. 


(43c) 
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Die durch (43b) definierte Hilfsfunktion y (u) hat eine einfache Form, die man am besten 
erkennt, wenn man 
1 
y?(u) (2 — a) 


itber ws auftragt. Das ist in Abb. 20 geschehen. Die Kurve ist fast linear und strebt monoton 
fallend vom Wert 0,113 zum Wert 0,093, wenn fu von 0 bis 1,7 zunimmt. Es gilt also in sehr 
guter Naherung 

9,4 


DW Th) Vw aa . 4Aa 
Ee Veme aie a ae 
Damit erhalten wir anstelle von (43a) und (43c) die Naiherungsbeziehungen 
~94R]/" oa Wee , 
Wee V2 (2— pw) (10 — p) ’ ey) 
EPS a —)7VR i 
A 6.27]/ E p=" ee) 
0,80 
Os 
80 
“ ie a 
0 7 —— 0 7 Z 
Abb. 21. o Berechnete Werte von = Abb. 22. o Berechnete Werte von a 
’ — Interpolationskurve. — Mittelwert 0,154. 
Die Beziehung (39) fiir die Mindestspaltweite hy;, schreiben wir in der Form 
| Ce = ATF o(u) (45a) 
mit 
nae 
o (11) = 12 Bn) Q(q) (45b) 
oder 
E hmin 
o(u) = vi(u) : (45c) 


a 2 B(u) ho 


Der Quotient h,,;,/hy, der in Abb. (21) dargestellt ist, nimmt von 4 =0 bis u =1,7 von 
0,8157 auf 0,7866 monoton ab, andert sich also in diesem Bereich nur um -+ 1,8%, wenn man 
auf einen Mittelwert bezieht. Das bedeutet, daB das Verhaltnis zwischen Mindestspaltweite 
und der Spaltweite am Beginn des Druckberges nur wenig von der Spaltverformung beeinfluBbt 
wird. 
In Abb. 22 ist 
1 


o(u) (2,7 — p) 


dargestellt. Diese Funktion ist fast konstant und schwankt im Bereich 0 Sw <1,7 nur um 
+ 1,8% um den Mittelwert 0,1540. Es gilt also fiir 0 <j <1,7 in sehr guter Naherung 


6,5 


o(M) =57 29 (45d) 
und 
_RnV 6,5 
nin = lo eres (45e) 


Die Beziehung fiir den maximalen Schmiermitteldruck schreiben wir in der Form 


] ips 
Dice =R | ay On) (46a) 


192 Dérr: Schmiermitteldruck und Randverformungen des Rollenlagers. Ingenieur-Archiv 


mit 


— 


iu) =<. (46b) 
y(u) f(s) de 


Trigt man die numerisch vorliegenden Werte von 


2—p 


8 (u) 


iiber us auf, wie es in Abb. 23 geschehen ist, so streuen die Punkte mit einem Gré8tfehler von 


+ 1,9% um die Gerade 4,78 (uw + 5,3). Daraus folgt fir0 Su S 1,7 die Naherungsbeziehung 


(yu) = 0,457 Va (46) 


mit einem Fehler von héchstens + 1%. 

Damit haben wir fiir alle Funktionen, die zur Ubertragung unserer numerischen Ergebnisse 
auf wirkliche Rollenlager erforderlich sind, handliche Naherungsformeln entwickelt. Diese 
Naherungsformeln sind etwa ebenso genau, wie die Eregbnisse des Iterationsverfahrens selber. 
Man darf annehmen, daB sie nicht nur fiir den untersuchten Parameterbereich, sondern auch 
noch fiir etwas gréBere w Werte, etwa bis uw = 1,85 giltig sind. 


a as tee 
0 7 vA 
Cre 
Abb. 23. o Berechnete Werte von res » — Ersatzfunktion 4,78 (u + 5,3). 


#(u) 


11. Auswertung und Zusammenfassung der Er- 
gebnisse. Fiir das verformbare, geschmierte Rollen- 
Jager werden Schmiermitteldruck und Randverfor- 
mung numerisch bestimmt, indem die nichtlineare 
Integralgleichung des Problems durch ein Iterations- 
verfahren mit technisch ausreichender Genauigkeit 
gelést wird. Die Schmiermittelstrémung wird dabei als laminar, die Viskositat als un- 
abhangig von Druck und Temperatur angenommen. Die relative Spaltverformung wird 
durch einen Verformungsparameter «u beschrieben, wobei der Bereich 0 <u <2 alle még- 
lichen Verformungszustande erfaBt. Mit dem angewandten Verfahren gelingt die Lésung 
des Problems fiir den Bereich 0 Sw <1,7, also fiir 85% des Gesamthereiches. 


Durch das Zusammenwirken von hydrodynamischen und elastizitatstheoretischen Erschei- 
nungen wirken beim Rollenlagerproblem zahlreiche physikalische Gré®en ineinander. Das 
erschwert die numerische Rechnung und auch das anschauliche Verstehen der formel maBigen 
Zusammenhange. Bei den mathematischen Herleitungen und numerischen Auswertungen 
werden deshalb die Langen, Driicke und Verformungen durch dimensionslose Variablen und 
Parameter ausgedriickt. Die so gefundenen Ergebnisse miissen fiir die Anwendung auf ein 
wirkliches Rollenlager durch gewisse Ubertragungsbeziehungen umgerechnet werden. Die 
gewahlte Definition des Verformungsparameters erlaubt es, fiir diese Ubertragungsbeziehungen 
einfache Naherungsformeln anzugeben, die allen technischen Genauigkeitsanspriichen geniigen. 
Die Anwendung dieser Ubertragungsbeziehungen, das hei®t der Gang der technischen Aus- 
wertung wird im Folgenden zusammengestellt. 


Wir betrachten die beiden Schmiermittelspalte zwischen Rolle und Schale sowie zwischen 
Rolle und Welle des in Abb. 24 dargestellten Systems. Die Elastizitatsmoduln und Radien von 
Rolle, Welle und Schale werden durch die Indizes r, w und s gekennzeichnet. Fiir den effek- 
tiven Elastizitatsmodul E und den effektiven Radius R gilt dann beim Spalt zwischen 


Abb, 24. Rolle zwischen Welle und Schale. 
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Rolle und Welle: Rolle und Schale: 
Rr Ru lisa AN 
R — uA x: S: 
Ry + Rw x Rs — R; 
Tre Ee Be 
th + Be Ek, 


Die Welle drehe sich mit der Winkelgeschwindigkeit m. Bei geringem Schlupf zwischen 
Rolle und Welle sowie zwischen Rolle und Schale gilt fiir die Summe der auf die Stelle kleinster 
Spaltweite bezogenen Randgeschwindigkeiten eines Spaltes 

V=V,+ ),.2o R, 


Die Welle driicke mit der Kraft P pro Einheitslinge auf die betrachtete Rolle. Dann ergibt 
sich fiir den Verformungsparameter der Wert 


wobei die Zahigkeit des Schmiermittels und y die Querdehnungszahl ist, die hier fiir Rolle, 
Welle und Schale als gleich angenommen wird. 

Als Mindestdicke des Schmiermittelfilms findet man 
— RY 6S 


min P oa — pb 


Der gréBte Schmiermitteldruck, das hei®t das Maximum des Druckherges ist 


204574? 2= a 
Poke SOR aan 
Der Abstand zwischen dem Orte des Druckmaximums und dem Ende des Druckherges auf 


der Spaltseite, wo sich der Rollenrand substantiell vom Rand der Welle bzw. Schale abhebt, 
ergibt sich zu 


V 1 
A=9,4 Ry” : 
P (2— ¢) (10 —#) 
Der Fehler dieser Naherungsbeziehungen ist héchstens etwa + 3% fiir den Parameter- 


bereich 0 Sw <1,7. Genauere Ergebnisse erhalt man, wenn man anstelle der Naherungs- 
forme] fiir ~ die in Abb. 19 eingetragenen genaueren numerischen Werte zugrunde legt. 

Mit A und pmax 14Bt sich die wirkliche Druckkurve q(x) aus den in Abb. 16 dargestellten 
normierten Druckkurven p(¢) fiir jeden w-Wert im Bereich 0 <y <1,7 bestimmen. Eine 
gegebenenfalls notwendige Interpolation bringt dabei nur geringfiigige Fehler, da sich die 
Kurven p (¢) nur wenig mit 4 andern. 

Auch die wirkliche Spaltform H(x)148t sich aus Abb. 18 leicht bestimmen, wenn / und hyj, 
bekannt sind. 

Die hier behandelten Probleme wurden fiir das Rollenlager mit kleiner Verformung, das 
heif®t fiir kleine p- Werte bereits von A. Meldahl, C. Weber und H. Blok mit Naherungsansatzen 
untersucht. Herr Prof. H. Blok (Techn. Hochsch. Delft) hatte die Freundlichkeit, einige der 
hier mitgeteilten genaueren Ergebnisse mit denen der drei genannten Autoren zu vergleichen. 
Die von A. Meldahl und C. Weber fiir u< 0,5 gefundenen Ergebnisse erwiesen sich dabei in 
leidlich guter Ubereinstimmung mit meinen. Die Ergebnisse von H. Blok stimmen bei kleinen 
p- Werten fast genau mit meinen iiberein, wahrend sie fiir wachsende y-Werte (u > 0,5) eine 
gréRer werdende systematische Abweichung aufweisen. 


(Eingegangen am 6. August 1953.) 


Anschrift des Verfassers: Dr.-Ing. Johannes Dérr, Institut fiir Praktische Mathematik der Technischen 
Hochschule (16) Darmstadt. 
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Rechnerischer Ausgleich von MeSkurvenscharen. 


Von H. Freytag. 


1. Einfithrung in die Aufgabenstellung. Bei der experimentellen Untersuchung der Ab- 
hangigkeit einer physikalischen MeBgréBe m von zwei von einander unabhangigen EinfluB- 
gréBen (z. B. von der Frequenz fund der Temperatur t) geht man gewohnlich so vor, da man 
die eine EKinflu&BgréBe jeweils wahrend einer MeBreihe konstant halt und die andere variiert. 
Welche EinfluBgré8e dabei innerhalb der MeBreihe konstant gehalten, und welche variiert 
wird, hingt von der Aufgabenstellung und der Einstellbarkeit und Reproduzierbarkeit der 
EinfluBgréBen ab. Auch bei der graphischen Darstellung der Funktion in Form einer Kurven- 
schar ist es im allgemeinen ohne Belang, welche EinfluBgréBe man als Abszisse wahlt und 
welche als Parameter. Diese Freiziigig- 
keit ist jedoch durchaus nicht immer 
gegeben. Wenn sich die zu unter- 
suchende GréBe min den betrachteten 
Bereichen relativ zu ihrem Anfangs- 
wert nur wenig andert, so wird es der 
meBtechnisch Erfahrene vermeiden, 
die Anderungen von m aus der Diffe- 
renz zweier ungefahr gleich grofer 
MeBwerte zu errechnen; er wird viel- 
mehr zur Steigerung der MeBgenauig- 
keit ein Verfahren wahlen, das die 
unmittelbare Messung der Differenz 
selbst gestattet. Es ist dann bei ge- 
schickter Versuchsanordnung meist 
méglich, den Gang der MeBgréBe m 
bei Variation einer EinfluBgréBe mit 
wesentlich gréBerer absoluter Genauig- 
keit zu bestimmen als den Anfangs- 
) f wert von mselbst. In diesem Fall ist 
i 73 “man in der Wah! des Parameters nicht 
mehr frei. Hat man beispielsweise den 
Frequenzgang1 der MeBgréBe bei 
mehreren konstanten Temperaturen aufgenommen, so mufs man auch beim Auftragen die 
Temperatur als Parameter wahlen, weil von jeder Kurve m = F, (f) die Gestalt (die durch 
die Ordinatendifferenzen der MeSpunkte definiert ist) wesentlich genauer bekannt 
ist als ihre Héhenlage; denn die Ordinaten jeder Kurve sind samtlich um den gleichen Betrag 
gefalscht, den (nach unserer Voraussetzung wesentlich gréReren) Me®fehler des Anfangswertes 
von m (vergl. hierzu Abb. 1). Dieser Fehler hat fiir jede Kurve einen anderen Wert. Es ware 
daher verfehlt, aus den MeBwerten des Frequenzganges bei verschiedenen festen Temperaturen 
Kurven des Temperaturganges bei verschiedenen festen Frequenzen bilden zu wollen, da die 
zur gleichen Frequenz gehérenden MeSpunkte untereinander zu sehr streuen. Man muB viel- 
mehr den Temperaturgang von m bei verschiedenen festen Frequenzen gesondert aufnehmen. 
Tragt man die Ergebnisse mit der Frequenz als Parameter in Kurvenform auf, so ist auch hier 
die Gestalt der Kurven genauer bekannt als ihre Hohenlage (vgl. Abb. 2). 

Wenn man nun derartige mittels eines Differenz-MeSverfahrens ermittelte Kurvenscharen 
des Frequenzganges und Temperaturganges vorliegen hat (die zweckmaBig unter Benutzung 
stets derselben MeBfrequenzen und MeStemperaturen aufgenommen werden), so sollte man 
annehmen, es kénne nicht schwer fallen, aus diesen Kurvenscharen wie bei rechnerisch ermit- 
telten Kurven das Liniennetz der gesamten Flache m = F(f, t) mittels gegenseitiger graphischer 


‘ Wir verwenden im folgenden der Kiirze halber als Beispiele der EinfluBgréBRen Frequenz und 
.Temperatur, ohne damit die Allgemeinheit der Gedankenginge einschranken zu wollen. 


4m 


= Mebwerte 
—-o — fehfertrele Kurve 


Abb. 1. Frequenzgang von m. 
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Interpolation geniigend genau zu bestimmen und darzustellen. Ein Versuch zeigt jedoch, daB 
dies nicht ohne weiteres méglich ist, weil die MeBwerte des Frequenzgangs und des Tempe- 
raturgangs einander zum Teil widersprechen. Man sieht dies am deutlichsten, wenn man die 
MeSpunkte des Temperaturgangs in das Diagramm des Frequenzgangs mit eintragt. Dies ist 
in Abb. 3 fiir 5x 4 MeBwerte ausgefiihrt. In dieser Abbildung und im folgenden Text sind zur 


leichteren Unterscheidung die aus Mes- 
sungen des Frequenzganges stammen- 
den MeBwerte von m mit q, die aus 
Messungen des Temperaturganges 
stammenden mit Tt bezeichnet. Der 
erste Index bezieht sich auf die Fre- 
quenz, der zweite auf die Temperatur. 

Wir stellen uns nun die Aufgabe, 
aus der Uberlagerung zweier offenbar 
fehlerbehafteter Netze von je 5x 4 
Me8Bpunkten das ausgeglichene, d.h. 
von Fehlern so gut wie moéglich befreite 
Netz abzuleiten. Dies erscheint auf 
den ersten Blick recht einfach. Da 
fiir jeden Zustand f =f;, t = zwei 
Me8ergebnisse vorliegen, liegt es nahe, 
jeweils deren arithmetisches Mittel als 
den wahrscheinlichsten Wert anzu- 
sehen. Bei einem derartigen Vorgehen 
wiirde aber die héhere Genauigkeit, 
mit der die Differenzen Ay und Ar ge- 
messen wurden, neben den gréBeren 
Fehlern der Anfangswerte von g und 
t nicht zur Geltung kommen kénnen; 
Fre quenz- und Temperaturgang von m 
wiirden auf diese Weise gréblich ver- 
falscht werden. Ein anderer nahe- 
liegender Gedanke ware der, die ein- 
zelnen Kurven des Frequenzgangs und 
des Temperaturgangs ohne Gestalt- 
anderung parallel zu sich selbst in 
Richtung der Ordinatenachse solange 
zu verschieben, bis g-Netz und t- Netz 
in Abb. 3 zur Deckung kommen. Ein 
derartiger Versuch mite zum Ziel 
fiihren, wenn die Differenzen Ay und 
Ar vollkommen fehlerfrei gemessen 
waren. Da jedoch auch diese Diffe- 
renzen gewisse (nach unserer Voraus- 
setzung allerdings kleine) Fehler auf- 
weisen, 14Bt sich der Ausgleich nicht 
ohne kleine Gestaltanderungen der 
Kurven durchfithren. Es leuchtet ein, 
daf diejenige Lésung die gréBte Wahr- 
scheinlichkeit fiir sich hat, bei der die 
Korrekturen an der Gestalt der Kur- 
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Abb. 3. Gegeniiberstellung der g- und t-MeS werte. 


ven am kleinsten sind. Ein solcher Ausgleich laBt sich freilich durch Probieren nicht ein- 
wandfrei durchfiihren. Deswegen wurde nach der Gaufschen Methode der kleinsten Quadrat- 
summe ein Rechenverfahren entwickelt, das den Ausgleich im obigen Sinne, frei von Willkiir 
und ohne zeitraubendes Probieren erméglicht. 

Es ist bekannt, daB der in der Industrie Tatige Ausgleichsrechnungen scheut, da sie in der 
Regel langwierige Rechnungen erfordern. Deswegen sei im Voraus bemerkt, dai das im folgen- 


14 
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den beschriebene Verfahren soweit durchgebildet wurde, daB es nicht mehr notwendig ist, fiir 
jeden Einzelfall die Normalgleichungen von neuem aufzustellen und in bekannter Weise durch 
schrittweise Ausmerzung der Unbekannten zu lésen. Es werden vielmehr fiir die an jedem 
MeBpunkt anzubringenden Verbesserungen einfache Formeln angegeben, mit denen sich die 
Verbesserungen schnell berechnen lassen. Das Rechenverfahren ist nicht auf 5x 4 Me8punkte 
beschrankt, sondern ohne weiteres auf beliebig viele MeBpunkte anwendbar. Abgesehen von 
der oben genannten Voraussetzung einer wesentlich héheren MeBgenauigkeit bei der Bestim- 
mung des Ganges der MeBgréBe wird an die Kurven nur die eine Bedingung gestellt, daB sie 
stetig verlaufen miissen. 


2. Aufstellung der Fehlergleichungen. Das von Gauf angegebene Verfahren zum rechne- 
rischen Ausgleich unsystematischer Fehler geht von der Forderung aus, da8 die Summe der 
Quadrate der anzubringenden Verbesserungen ein Minimum werden soll. Da es in unserem 
Falle darauf ankommt, die Gestalt der Kurven méglichst wenig zu andern, haben wir als 
Verbesserungen v im Gaufschen Sinne die Unterschiede zwischen den Ordinatendifferenzen 
der ausgeglichenen Kuive und der ihr entsprechenden g- bzw. t-Kurve anzusetzen. Bezeich- 
nen wir die ausgeglichenen MeBgréRen mit r und benummern die g-, t- und r-Werte wie in 
Abb. 1 und 2, so ist demnach zu fordern, daB folgende Differenzen paarweise méglichst genau 
iibereinstimmen: 


Tro To Mit Py —Poo 3 T20—Too Mit Px9—Poo USW- 

aber auch Toi —T on Wit Toy—To 3 Tox Top WMIt T.—To) USW. 
Wenn die y- und t-Messungen bei a MeBfrequenzen und b MeBtemperaturen ausgefiihrt 
wurden, so lassen sich auf diese Weise an Hand der 2 a b Punkte zah!reiche Fehlergleichungen 


aufstellen, und zwar ergeben sich b (a — 1) Gleichungen aus den g-Differenzen und a (b — 1) 
Gleichungen aus den t-Differenzen. 


Fehlergleichungen fiir a = 4 MeBfrequenzen und b = 5 MefBtemperaturen 
(E20 3) clea Ope el 


p10 = (Tro Too) (Pio — Poo) > ror = (To1 Too) (To Too) » 

Up 29 = (T20 — Too) (P20 — Poo) > roe = (Toe Too) (To2 — Too) 
Up 30 = (T30— Too) — (P30 — Yoo) > Yos = (To3 — Too) (To3 — Too) » 
Voir = (T11 — M1) (P11 — Yor) > eos = (Toa — Too) (Tos — Too) > 
Vpor = (To To1) (G1 — Pu)> Yer = (M11 — No) (11 — Ty) 
Uys = (T31— To) — (P31 — Yor) > x12 = (M12 — No) (T12 — To) 
Upi2 = (T12 — To2) (P12 — or) > %13 = (13 — Mo) — (%13 — To) > 
: : Ura = (Tyra — M10) — (T1a — To) > 
(T21 == Too) — (oq — Teo) 


Usa = (34 — Toa) — (P34 — Poa) » 
(15 Gleichungen) (16 Gleichungen) 
Bei den Klammerausdriicken in den obigen Fehlergleichungen handelt es sich um Differen- 
zen nahezu gleich groBer Zahlen. Um bei den spateren Zahlenrechnungen mit einer geringen 


Stellenzahl] auskommen zu kénnen, fiihren wir andere, kleinere Gré®en ein. Wir gehen dabei 
von der halben Summe und der halben Differenz der Werte 1;, und Yi, aus. ndmlich von 


Sih = 0,5 (tin + Qin) und. dj, = 0,5 (ti — on), 
und erhalten, indem wir 


Tih = Sik + Up , 
setzen, die Beziehungen 


Tk — Dik = Uy + dj, , Nik — Gk = Uj, — dyy . 
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Dann bekommen die Fehlergleichungen folgende Form: 
Vo19 = (Uyo Ugo) + (dio doo) » U01 = (Uo, Ugo) (do d 
Vp 09 = (Ug9 — Ugo) + (449 — do), og = (tog — Upo) (dy. — doo) 
p30 = (39 Ugo) + (439 doo) + %r03 = (Ups Ugo) (dj; — d 


Vp 34 = (Uzq — Uy) + (ds4— dyy), 2, me ae (ws, U4o) Ca d4). 


Die Koeffizienten-Matrizen fiir u;;, und d;, der Fehlergleichungen sind in den Tabellen 1 und 2 
wiedergegeben. 


Uos Uog Uo 
} 


Wir haben nun nach Gauf zu fordern, da8 Lv? ein Minimum wird. Wenn die MeBgenauig- 
keit bei der Bestimmung der Differenzen Ag und At verschieden gro ist, muB man den Ver- 
besserungen v, und v, verschiedenes ,,Gewicht“ beilegen. Das Gewicht eines Fehlers ist dem 
Quadrat des mittleren Fehlers (der aus Messungen ermittelt oder auch geschatzt werden kann) 
umgekehrt proportional. Ordnen wir den v,-Werten das Gewicht 1 zu und sind uw, und yu, die 
mittleren Fehler bei der Messung von Ap bzw. Art, so ergibt sich als Gewicht der v,- Werte 


P= Pat Mg - 
Wir haben demnach zu fordern 
S=2Lpwuy+ Lv = Minimum . 
3. Aufstellung der Normalgleichungen. Die Fehlergleichungen enthalten einerseits die zu 
ermittelnden Zuschlage u, die an den mittleren Ordinaten s anzubringen sind, andererseits die 


bekannten halben Differenzen d zwischen je zwei zum gleichen Zustand gehérenden g- und t- 
Werten. Als Verdnderliche der Minimumsaufgabe haben wir die u anzusehen. Aus der Be- 


14* 
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Tabelle 2. Koeffizienten- Matrix (h) in den Fehlergleichungen. 
EN eS 
dy, d,, dys di, yg dy, dy» dys | dy, | dy ds, dy dy, ds, 
(eee 


Ay, | dio 


V943 
Voogs 
V033 


Vpy4 —l ] 1 
Vopg4 —l | 


Vp34 


dingung, daB die partiellen Ableitungen der Quadratsumme S nach allen u;; gleich Null werden 
miissen, erhalten wir a - b sog. ,, Normalgleichungen“ fiir die a - b Unbekannten u;,. 
Bei der partiellen Ableitung benutzt man zweckmafig die Umformung 


os A(X pve) 0(= vz) c ave dur | a %a)| 
TT LP ep 


In den Tabellen 3 und 4 sind die Koeffizienten-Matrizen der Unbekannten u;, und der 
Konstanten d;;, fiir alle a: b Normalgleichungen wiedergegeben. Das Bildungsgesetz in Ab- 
hangigkeit von aund bist deutlich zuerkennen. AuBerdem ist zu ersehen, daB die Summe aller 
Koeffizienten in jeder Spalte u;, und d;; gleich Null ist. Das bedeutet aber, da8 nur (a b — 1) 
Gleichungen von einander unabhangig sind. Denn wir kénnen jede der a- 6 Gleichungen auch 
dadurch erhalten, daB wir alle anderen addieren. Wir haben somit gerade eine Gleichung zu 
wenig. Dies ist nicht verwunderlich, wenn man die u-Koeffizienten-Matrix der Fehlergleichun- 
gen betrachtet. In allen Fehlergleichungen kommen die Unbekannten u immer nur in Form 
von Differenzen vor, niemals in anderer Kombination. Da die Normalgleichungen lineare 
Kombinationen der Fehlergleichungen sind (die partiellen Ableitungen liefern nur Konstante), 
liefern auch die Normalgleichungen nur Bedingungen fiir die Differenzen je zweier Unbekann- 
ten, niemals fiir die Unbekannten selbst. Wir miissen daher noch eine Bedingung hinzufiigen, 
durch die die Héhenlage des ausgeglichenen Punkthaufens festgelegt wird. Wir wahlen die 
Bedingung 2 u = 0, die besagt, da der Schwerpunkt des ausgeglichenen Punkthaufens 7;; 
mit dem Schwerpunkt der Mittelwerte s;, der y- und t-Werte zusammenfallen soll. Damit 
schaffen wir gleichzeitig die Voraussetzung dafiir, daB auch XY u? ein Minimum wird. 
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4, Allgemeine Auflésung der Normalgleichungen. Bei Ausgleichsrechnungen ist es iiblich, 
das System der Normalgleichungen nach dem Einsetzen der Zahlenwerte fiir p und alle d von 
Fall zu Fall numerisch aufzulésen. Da man es bei nur drei g- und drei t-Kurven schon mit 
9 Gleichungen, bei fiinf y- und vier t-Kurven aber mit 20 Gleichungen zu tun hat, ist dies eine 
ziemlich langwierige und Ubung voraussetzende Arbeit. Um diese immer wiederkehrenden 
Rechnungen zu vermeiden, wurde das System von a:b Gleichungen generell algebraisch auf- 
gelést. Im folgenden werden die Ergebnisse mitgeteilt. Von einer auch nur auszugsweisen 
Wiedergabe der durchaus elementaren, aber ziemlich umfangreichen Berechnungen muB aus 
Raumgriinden abgesehen werden. 


Die an den Mittelwerten s;, anzubringenden Zuschlage u;,, durch die das Netz entsprechend 
unseren Forderungen ausgeglichen wird, sind, wie zu erwarten war, lineare Funktionen der 
halben Differenzen d;,, und zwar tragt im allgemeinsten Fal] jede dieser halben Differenzen 
nach MaBgabe eines Koeffizienten c zu jedem Zuschlag u etwas bei. An dem Mittelwert sp, , 
der der Frequenz f und der Temperatur t zugeordnet ist, ist ein Zuschlag 


anzubringen. Die Koeffizienten elt) sind algebraische Funktionen von a, b, und p. Das 


zweite Indexpaar f,t weist darauf hin, daB sich c nicht nur nach der jeweiligen Differenz dj, 
richtet, sondern auch nach der Verbesserung uf,, fiir die es gelten soll. Die Mannigfaltigkeit ist 
jedoch keineswegs so groB, wie es zundchst scheint. Unter den a: b Koeffizienten ¢ gibt es im 
Héchstfall 16 der Formel und dem Zahlenwert nach verschiedene. Die Formeln fiir die 16 
Koeffizienten, die zur leichteren Unterscheidung mit Buchstaben des groBen Alphabetes 
(A bis R, unter Auslassung von O) bezeichnet werden, sind in Tabelle 5 zusammengestellt. 
Die Zahlenwerte von A... Rsind fiir einige Werte von a, b, und p der Tabelle 6 zu entnehmen. 


Tabelle 5, Formeln fiir die Koeffizienten A... R. 


VA =—(p+1)(bp +1) (a + p)[(1 + ab) (e —bp) + 2 4b(p—1)], 
VB=—(p+1)(bp+1)(a+ p)(2abp +a—bp), 
VC =—(p+1)(bp + 1) [ab(a + bp) (a—p—2) + 2abp(p +1) +(a+4+ p)(a—bp)], 
VD =—(p+1)(bp +1) [2abp(p+ 1) + (a+ p)(a—bp)], 
VE =-+ (p+1)(bp+1)(a+p)(2ab—a-+ bp), 
VF =—(p+1)(bp +1) (a+ p)(a—bp), 
VG =+(p+1)(bp +1) [2ab(a + bp) —bp(a—p) —a(a + p)], 
VH =—(p+1)(bp+1)[(e+ p)(a—bp) +2abp], 
VJ = + (p +1) (a + p) [ab(a + bp) (bp —1) —bp(a+ bp) (2a—1) + 2ab—a+ bp], 
VK =—(p+1)(a+ p)[bp(a+bp)(2a—1) +a—bp], 
i 


Ser ens 


[ab(p —1) (a+ p) (a + bp) —(p + 1) (a —bp) (a + p) + 2ab(a + bp —p)] 
= 24D p? (o> p)(ar bh p= i). 


VM =— 2abp[bp+1+p(a+ p)(a+bp—1)]—(p +1) (bp +1) (a4 p) (a — bp), 
VN =+ (p+1)(a+ p)[2a(b—1) + (bp +1) (a + bp)], 

VP =—(p+1)(a+p)[(6p + 1) (a —b p) —2abp}, 

VQ=-+ 2ab[(bp +1) (a+ bp—p) + p?(a+ p)] —(p + 1) (bp + 1) (a + p) (a — bp), 
VR=—2abp[bp+1—p(a+ p)]—(p+1)(bp + 1) (a+ p) (a —bp). 


V=ab(p+1)(bp +1) (a+ p)(a + 5p). 


5. Die Struktur der d;;,- Koeffizienten-Matrix. Die Zuordnung der 16 Koeffizienten 4A...R 
zu den halben Differenzen d;, und den zu berechnenden Zuschl agen u zeigt die Tabelle 7, die 
fiir den Fall von a = 5 Mef®frequenzen und b = 6 MeBtemperaturen aufgestellt ist. Die Koeffi- 
zienten-Matrix ist quadratisch und symmetrisch zur Hauptdiagonale. Sie weist auBerdem 
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einige GesetzmaBigkeiten des Aufbaus auf, 
deren Kenntnis fiir die Umstellung auf andere 
Werte von aund 6 von Wichtigkeit ist. Sie zer- 
fallt namlich in a* zum Teil gleiche Teil quadrate 
von der Seitenlinge b, die ebenfalls symmetrisch 
zur Hauptdiagonale aufgebaut sind. Bei der 
gesamten Matrix wie bei den Teil quadraten kann 
man einige Teilflachen unterscheiden, die in 
ihrem Koeffizienten-Inhalt einheitlich auf- 
gebaut sind. Diese sind: ein quadratisches 
Eckfeld (links oben), ein Streifen in der Haupt- 
diagonale, zwei Randleisten (oben und links) 
sowie zwei Dreieckflachen. Dabei hat man 
als Struktureinheit ein Quadrat anzusehen, 
das bei der gesamten Matrix, die Kantenlange 
b, beim Teilquadrat die Kantenlange 1 hat. 


6. Praktische Durchfiihrung der Ausgleichs- 
rechnung. Gegeben seien aus Messungen! der 
Frequenzgang der MeBgriBe m = F, (f) bei b 
konstanten Temperaturen t)...t,_, (b Kurven 
des y-Diagramms, aus denen bei den Abszissen 
fo---+ fay insgesamt a-b MeBwerte g;, ent- 
nommen werden), der Temperaturgang der 
MeBBgréBe m = F,(t) bei a konstanten Fre- 
amenzen jf,...f,, (a Kurven des 1-Dia- 
gramms, aus denen bei den Abszissen ft)... th_4 
insgesamt a:b MeBwerte 1;, entnommen 
werden). 

Ferner sei (aus Messungen oder Schatzun- 
gen) die mittleren Fehler uw, und pw, bei der 
Ermittlung des Frequenzganges und des 
Temperaturgangs bekannt, woraus p = [dz : Mg 
errechnet werden kann. Man geht dann wie 
folgt vor: Aus den 2x a-b MefSwerten 9, 
und 7;; werden die a-b Mittelwerte 


Sik = 0,5 (tie + Qix) 
und ebensoviele Differenzen 


d;, = 0,5 (Tit = Pik) 
berechnet. 

Dann stellt man an Hand der Regeln des 
Abschnittes 5 und der Beispiel matrix Tabelle 7 
die d-Koeffizienten- Matrix fiir die betreffenden 
Werte von aund b auf, und zwar zunachst nur 
mit Buchstaben (A... R). Nunmehr berechnet 
man diese Koeffizienten A... R entweder 
aus der Tabelle 7 (aus den Zahlenwerten 


U-(A... R) und U) 


oder aus der Formelzusammenstellung Tabelle 5. 


Die einzelnen Summanden, aus denen sich 
die us, zusammensetzen, ergeben sich dann 
aus der Multiplikation der Zahlenwerte der in 


1 Entsprechend dem oben gewahlten Beispiel 
sind als EinfluBgréBen wieder Frequenz und Tempe- 
ratur gewahlt, 


TRY. 


Zahlenwerte der Koeffizienten A... 


Tabelle 6. 
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hens zu addieren. 
Rkann dienen, daB 


... R und des dariiber im Kopf der Tabelle ge- 
) simtliche Zuschlage u verschwinden miissen. 


sind dann unter Beachtung ihres Vorzeic 


, stehenden Koeffizienten A 


a 
o 
I 
= 
oc 
= 
q 
= 
WN 
gl 
4 
ar 
ts: oS 
= 4 
a8 
N 3 
wa 
o 
ae 


ee 


e Berechnung der Koeffizienten Ax 
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Als Kontrolle fir die richtig 


fiir lauter gleiche d (z. B. 


Tabelle 7. 
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Beitrag zur Potentialstr6mung durch radiale Schaufelgitter. 
Von W.-H. Isay. 


1, Einleitung. In neuerer Zeit erfolgt die Berechnung der Potentialstrémung durch axiale 
Schaufelgitter immer mehr durch die Anwendung der direkten Singularitatenmethode, indem 
die von Birnbaum! seinerzeit fiir den einzelnen Tragfliigel geschaffene Theorie der Pageuden 
Wirbelflache auf ebene Schaufelgitter ausgedehnt wurde. Diese Theorie ist sehr iibersichtlich 
und es ist méglich, mit ihr Kinzelgitter sowie Gitterstufen mit beliebigen in der Praxis vor: 
kommenden Profilformen zu behandeln 2; ferner kann sie in beschranktem Umfange auch bei 
kompressiblen Strémungen angewendet werden 3, >, : 

Ks ist somit naheliegend, diese Theorie auch zur Berechnung der Potentialstrémung durch 
radiale Schaufelgitter (Schaufelsterne) heranzuziehen. Dabei wird es sich zeigen, da8 alle in 
Arbeit I fiir axiale Gitter gewonnenen Ergebnisse auch auf Radialgitter tibertragen werden 

énnen, 


2. Das Einzelgitter mit linienhaften Profilen. Wir betrachten in der komplexen z-Ebene 
ein Radialgitter mit N linienhaften (,,unendlich diinnen‘) Profilen. Wir wollen annehmen, es 
sei ein Laufradgitter (Abb. 1). Die Strémung durch das mit der Winkelgeschwindigkeit w 
rotierende Schaufelgitter ist nur im Absolutsystem eine Potentialstrémung, wahrend die 
Relativstrémung einen konstanten Wirbel (— @) besitzt 4. Wir beziehen uns daher auf das 
Absolutsystem. Dabei denken wir uns die Strémung durch das Gitter entstanden durch 
Uberlagerung der Strémung (u,, v,), die von den Zirkulationsverteilungen auf den Profilkon- 
turen induziert wird und der Strémung (up, vp) einer im Mittelpunkt des Gitters (Nullpunkt) 
angebrachten Quelle (oder Senke) Q) und einer Zirkulation J) (innere Zirkulation). Es ist 
also § 


. 1 . 1 
Co(2) = Up — tM =5—(Q+ily)—. (1) 
ae 
Nun ist das von N kreisférmig in den Punkten €,¢ e ae 
ary Nhe Y 
§22—_ 
Ce S angeordneten Wirbelelementen y ds induzierte Stré- Ss Ys, 2a Ps 
mungsfeld im Punkte z gegeben durch? Sa Grit 
A t N2zN-1 
de, = du, — i dv, Fe! ds EN 
Damit ist das resultierende von der auf der Profilkontur an- 
gebrachten Zirkulationsverteilung y(s) (s = Bogenlange) indu- 
zierte Geschwindigkeitsfeld 
: i N 2N-1 
¢,(z) = u, — ty = | v(s) NEN ds. (2) Ret 


Wir setzen z = x+iy =re'? und €=&+in = ge'” und geben die Kurve der Schaufel- 
kontur in der Form r (¢) bzw. 0 (y) vor. Es sei r (py) ein Pol ynom bzw. ein Fourierpol ynom in 
, jedenfalls jedoch eine eindeutige Funktion von 9. 

Da wir y nur auf der Kurve o (y) anbringen, kénnen wir statt y(s) schreiben y (y). Ferner ist 


ds =o? -+ ody. 


1 W. Birnbaum u. W. Ackermann, Z. angew. Math. Mech, 3 (1923) S. 290. 

2 W. H. Isay, Beitrag zur Potentialstrémung durch axiale Schaufelgitter. Diss. T. U. Berlin 1952, 
Z. angew. Math. Mech. 33 (1953) 5S. 397, sowie die hier verzeichnete Literatur. Diese Arbeit wird als Arbeit I 
bezeichnet. 

3G. V. R. Rao, Journ. of the Aeronautical Sciences 1953 S. 217. 

4 Vel. z. B. W. Spannhake, Hydraulische Probleme, S. 180, Berlin VDI 1926. 

5 P. F. Byrd u. M. T. Huggins, Ing.-Arch, 21 (1953) S. 191. 

6 Allgemein bedeuten hier u und v die x- und y-Komponente der Absolutgeschwindigkeit ¢(c? = u?-+ v?), 

7 Vgl. z. B. W. Miiller, Mathematische Strémungslebre, S. 97, Berlin 1928. 
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Damit ergeben sich die x- und die y-Komponente des von der Zirkulationsverteilung y () 
(die auf der Kontur 0(y) angebracht ist) induzierten Geschwindigkeitsfeldes zu 


N 
) sin(Ngy —Ny —p) +singp 


N 72 
u(t, 9) =s—-— | y(y) = sw Vo? t ody, (3) 
By Tee cos N(p —y) + (20) 
N 
N a ee cos (Ny — Ny —g) —cos 9 wa = 4) 
Vy Ley =) YY N y)\\2 N Q @ Y: 
“ eat cos Nip —w) +(20) 
Ferner ist nach (1) 
Lsj 
Ta es (5) 
i le 
NO epee (6 


Es ist jetzt die folgende Aufgabe gestellt: Auf der Schaufelkontur r(q) soll eine Zirkulations- 
verteilung y (y) so bestimmt werden, daB r (y) Stromlinie der Relativstrémung wird. 

Da wir uns auf das Absolutsystem beziehen, bedeutet das fiir unser Problem die folgende 
Randbedingung (vgl. Spannhake): An jeder Stelle der Profilkontur muB die Normalkompo- 
nente der Geschwindigkeit c) + ¢, itbereinstimmen mit der Ge- 
schwindigkeitskomponente des Schaufelelementes normal zu 
diesem selbst. Ist 


(p) sin p + r(p) cos Y 
(p) cos p — r(y) sin g 


der Anstieg der Profilkontur, so fiihrt die Randbedingung auf 
die Relation (Abb. 2) 


w rcos (8 — ~) = (v% + »,) cos i — (uy + u,) sin #, 


r’ 
tg 0 oO 


Abb, 2, 


oder 
wr(p) cos y {r'(~) cos p — r(y) sin gy} + w r(¢) sin p {r'(y) sin p + r(p) cos p} 
= [v0 (r(9), 9) + % (119). PI Lr'(p) cos » — r(g) sin G7] 
— [wo (r(¢), 9) + % (r(p) 91 [F'(y) sin @ + r(g) 008g] - 
Diese Relation stellt eine Integralgleichung fiir y (y) dar. Sie lautet, wenn wir noch zur Ab- 
kiirzung 


v(y) Vo? + 0? = a(y) 


setzen, 


— a r(g) r’ (9) 


Ae An r’(#) (ow) \™ 
ate Gag, sin N (py — yp) — A a cos N (py — yp) 


dy. 
sy Fn) N(p —y) + ae 2 4 
1) r(@) | 


Dies ist die Integralgleichung unseres Problems. Die einzige Singularitat ihres Kernes liegt an 


der Stelle m = y. 
Bezeichnen wir den Kern von (7) mit K (g, y), so wird 


cit 


r’ rN_oN 0 N 
lim ( ) K Ni ese ee a ee 
im — = j ee 
se ee (py) =5_ lim yea = 


6 2 
p>y mai : 
( 2) a) ar rN oN N? 
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Demnach kénnen wir K (9, y) auf die Form 


Kip.) =, : + Hp.) 


ese 
bringen, wobei H (¢, y) im ganzen Bereich in beiden Variablen stetig ist. Speziell ergibt sich, 
wenn man den Grenziibergang durchfihrt, 


, ,\3 (oa av} 
NU 4+(N »(") pit 


1] r 
H(p, 9) =5 SO - (8) 
1+ () 
- 
Mit 
Ne ae A Ns I, r’(~) Qy 
DD gee OTT AG) a 
nimmt (7) die Form an 
p) — I A ) L 4) 
fo) =z | ow E>, + He 9] ay. (9) 


Die Auflésung dieser Integralgleichung erfolgt nach der von W. Schmeidler in seinem 
Buch fiir Integralgleichungen erster Art angegebenen Methode der Kernmatrix und wurde 
bereits in der Arbeit I behandelt, so da® wir hier nicht mehr darauf eingehen!. Die in der 
Lésung o (wy) enthaltene freie Konstante 


iP = J ow) dy = [yy Vo + ody 


wird durch die AbfluBbedingung 

Via) = 0 
bestimmt. Die sogenannte innere Zirkulation J kann durch die Bedingung des stoBfreien 
Eintrittes 


ya) == 0 
festgelegt werden. Die auBere Zirkulation des Gitters ist dann 
Mtn. 


Damit 14Bt sich aus (3), (4), (5), (6) das Geschwindigkeitsfeld der Gitterstrémung in jedem 
Punkt (r, y) der Gitterebene berechnen. 

Fir die Berechnung der Druckverteilung am Profil (hier natiirlich nach der Bernoullischen 
Gleichung im Relativsystem) und des Drucksprunges von der Druckseite zur Saugseite des 
Profiles, mu man noch die resultierende Relativgeschwindigkeit w an diesen beiden Seiten 
ermitteln. Es ist 


w? = (uy + u, + orsing + Au,)? + (v) + v4, —wreosp + Av,)*, 
und zwar (++) an der Saugseite und (—) an der Druckseite. Dabei ist 


_ __—s¥ ~ cos Y—rsing _ __—«y rsing +7 cos p 
AOS SO pa rare ae aa a re 


Liegt nicht ein Laufradgitter sondern ein Leitradgitter vor, so behalten alle Formeln unver- 
andert ihre Giltigkeit, sofern nur iiberall m = 0 gesetzt wird. 


3. Das Einzelgitter mit endlich dicken Profilen. Wir betrachten ein einzelnes Schaufel- 
gitter mit endlich dicken Profilen (Abb. 3). Es mége sich wieder um ein Laufradgitter handeln. 
Der Fall des Leitrades 1a8t sich dann sofort entnehmen, indem man w = 0 setzt. Die Kurve 
der Profilkontur, auf der wir die Zirkulationsverteilung anbringen, sei in diesem Fall in der 
Parameterdarstellung 


M 
y(t) = — acost , rt), 2. (= 7 t= a) 
v= —M 


1 Der singulare Hauptteil des Kernes ist wohlbekannt. (Vgl. W. Schmeidler, Integralgleichungen I, 
Leipzig 1950, S. 51). Der stetige Zusatzkern wird durch ein Fourierpolynom (py = —a@ cost gesetzt) 
approximiert, so da die Berechnung und Auswertung der Kernmatrix méglich ist. 
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gegeben. Die c, hingen von der speziellen vorliegenden Form des Profils ab. Bei Profilen mit 
spitzer Hinterkante (Abb. 3) mu8 7 (0) = 0 sein. Bei abgerundeter Hinterkante ist r (0) += 0. 
Im iibrigen gelten die Betrachtungen von Ziff. 2 sinngem4B auch hier, insbesondere kénnen wir 
die Formulierung der Randbedingung wortlich iibernehmen. Mit 


__7(t) sin p(t) + r(t) p(t) cos ¢(t) 
8 OY) =F) cos lt) —r(t) ot) sin 92) 


als Anstieg der Profilkontur lautet sie 
o» r(t) cos (O(t) — g(t) = Lv (r(0), g(t) + 2, (x(t), g(t) ] e08 8(2) 
— [ug (r(t) . (0) + u, (x(t) « p(d)) sin O(0) 


Indem wir 


olny). = E 
sin (Ng — Ny(t) —¢/) + sin yp 


r r 


N r . ° 
u(r, 9) = y(t) onl ory? + edz, (10) 
as td pe 2 (Ar) 20s N (vp — (t)) + ae 
K N 
N eel cos (Ng — N y(t) —¢) — cos p / 
ar, @) = y(t) + ory? + or dr (11) 
: a il —2(242))" cos N (vy — y(z)) ae 


— 


einsetzen und noch zur Abkiirzung 


y(t) Vo2p? + 2? = a(t) 


schreiben, erhalten wir fiir unser Problem die folgende Integralgleichung: 


c eerie 
— o r(t) F(t) —52 ua — 2 5(2) 
; Reps F(t) (o(t)\N 
nN a + H(t) i) sin N (p(t) — ¥(2)) ae cos N (p(t) — ¥(2)) (i) 
ez aD 1—2 BT cost pe) = (eee Gah a 
av r(t) Ce as r(1) 
Der Kern von (12) ist genau dann singular, wenn zugleich gp = y und r = gist. Das ist immer 
der Fall, wenn t = 7 ist; auBerdem wenn t = 2 und t = — a oder umgekehrt ist. 
Bezeichnen wir den Kern von (12) mit K (t, 7), so wird 
N—oN 
a ¢pn-17 8 4 pw on NGL —* 
lim tg’ * K(t, 7) =lim> 1 prs cat ae t—t ae 
fees top 2 t —T 2 =e Sorts Sel ree 4x 
t—T t—t 
Ferner ist 
: Co é t—T 1 
Le 3 5 T(E) ae rer K(t, 7) Sores © 


da r(--z) = r(x) und y(— 2) = 7 (az) ist. Wir kénnen also 
7 K(t,t) auf die Form bringen: 


K(t,7) = etgt=* ce H(.}. 


An 


Dabei ist H (t, 7) im ganzen Bereich in beiden Variablen stetig. 
Die Durchfithrung des Grenziiberganges liefert noch 


* ry \3 GQ, oA 
vp +n toes v—a(F) spas 
(et) : 


cawanmeeie 
Abb. 3, t alg! 


(13) 
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Da der Kern von (12) eine vollstandige Ableitung nach t darstellt, gilt 
fKe.2) dt—0, 
Daraus folgt, daB auch 7: 
Ec) dt 0 (14) 


sein mu, da von selbst 


ist. Mit 


1aBt sich somit (12) auf die Form bringen: 


fi) = [ow jets" + H(t,7)| dr. (15) 
Dabei ist H (t,t) an die Bedingung (14) gebunden. 
Fiir die Auflésung dieser Integralgleichung verweisen wir auf Arbeit I1. Die in der Lisung 
enthaltene freie Konstante 


i = foe) dt = frole wy? + 0% dr 


wird durch die AbfluBbedingung y (0) = 0 bestimmt. Die innere Zirkulation [, kann z. B. 
festgelegt werden durch die Wah] einer Anstrémrichtung. 

Damit 148t sich aus (3), (4), (10), (11) das Geschwindigkeitsfeld der Gitterstrémung in 
jedem Punkt (r,q) der Gitterebene.berechnen. Fir die Berechnung der Druckverteilung 
auBen am Profi] mu’ man die resultierende Relativgeschwindigkeit an der AuBenseite des 
Profils kennen. Ihr Quadrat im Punkte t ist gegeben durch 


w(t) = (up(t) + u,(t) + @ r(t) sin y(t) + Au,(t))? + (vo(t) + 0,(t) — @ r(t) cos p(t) + Av,(t))? 


mit 


_ _ «yy ~reosp—roepsing 2 y rsing +rgceos@p 
ae 2 7? + 7? g : Av, 2 Vr? + rp : 

Y r: 7 
Der Drucksprung von der Druckseite zur Saugseite des Profils ergibt sich hier automatisch fiir 
die entsprechenden (voneinander verschiedenen) t- Werte. 


4. Die Gitterstufe mit linienhaften Profilen. Wir kommen jetzt zur Behandlung einer 
Gitterstufe bestehend aus einem Leitradgitter und einem Laufradgitter (Abb. 4). Wir beziehen 
uns natiirlich auf das Absolutsystem und betrachten eine ,,Momentaufnahme“ der Stufe. Die 
Konturkurven der Schaufeln seien in der Form r, (¢) 3 T2 (P) als stetige, eindeutige Funktionen 
gegeben. Bei allgemeinster Stellung von Laufrad und Leitrad gegeneinander ist 


— 4 SH S44, b—a, SY, Sb—a,+ 24). 


Wegen 6 > 0 ist ferner stets 1,}< 1. 
Auf der Profilkurve r, (y,) sei die Zirkulationsverteilung y, (y,) und auf der Kurve 1, (~,) 
die Zirkulationsverteilung y, (y,) angebracht. Es ist dann z. B. 


y f Siero: 
Uyq(Tp> Pp) mets ¥q(Ya) og\N oq\2 NV Voy + 0? ty, 
, La (A) cos N (yr — vo) + ( 


N 
(22 sin (NV pp— N yq — #p) + 8i0 yp 
) 


Tp 


Tp 
(9) f 


1 Vgl. auch FuBnote von S. 205 und in diesem Falle: W. Schmeidler, Integralgleichungen I, S. 37, 
Leipzig 1950. 
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die x-Komponente der von der auf der Kurve 1(9q) angebrachten Zirkulationsverteilung Y in 
Punkten der Kurve 1,(y,) induzierte Geschwindigkeit (p,q = 1,2). Unsere Aufgabe ist es 
nun, die Zirkulationsverteilungen y, und 7, so zu bestimmen, daB r, (~,) Stromlinie in Bezug 
auf das absolute Geschwindigkeitsfeld und r,(p,) Stromlinie in Bezug auf das relative Ge- 
schwindigkeitsfeld wird (vgl. auch die Formulierung der Randbedingung in Ziff. 2). Diese 
Aufgabe fiihrt auf folgende Bedingungsgleichungen, die ein System von Integralgleichungen 
fiir y, und y, darstellen: 


= [v (71(71) 2 (1) + %1 CAC) > $1) oe Uy o(T1(1) ; 71)] cos A, 
aa [uo (71(1) ’ (1) + yy (r1(91) ’ 1) + Uy» (r1(91) 7 ~1)) sims 
( Ty(P2) C08 (Fy — Po) = [% ( P2) » 2) + Yy1 (72(P2) ? $2) + Ye (Tp ie ? eo cos B, 
up) 


(16) 
( 
a [Up (T2(@2) ’ 2) + Uy, (r3(P2) ° 2) a Uy » ( Qo) » 2) | sin 0. 


Leitred Le/trad 


Lautrad 


Laufrad ww =fe/ativgeschw. 


w= Relativgeschw. 
c=Absolurgeschn. 


¢ =Absolutgeschw. 


Abb. 4. Abb. 5. 


Eine Betrachtung ganz analog zu der in Ziff. 2 durchgefiihrten zeigt, da® (16) auf die Form 
gebracht werden kann. 


a, b—a,+2a, 
1 1 1 
fl) == [a gay + Balm wh an +f only) Halos +9) ays. 
—a, b—a, 
a, b—a,+2a (17) 
Ge eee H Appts Bl bes | \ 
2(P2) = es | O1(Y1) Hoy (P2 > Yr) dp, + oa i J(Po) | 2 — 2 + Ay(Q2 ; Ye) | dy, . 
—4a, b—a, 


Alle H,, sind stetig in beiden Variablen (p, q = 1,2). Fir H,, und H,, gilt Formel (8). Fiir die 
Auflésung des Systems (17) sei auf Arbeit I verwiesen!. Die in den Lésungen o, und o, enthal- 
tenen freien Konstanten 


ay b—a,+2a, 
ry = ACH dy, , r, = I O2(P2) dps 


werden durch die AbfluBbedingungen y, (a,) = 0 und y, (b — a,) = 0 festgelegt. Damit ist die 
Berechnung des Geschwindigkeitsfeldes méglich. 


Die Ermittelung der Druckverteilung an den Profilen erfolgt prinzipiell wie beim Einzel- 
gitter. (Fiir das Laufrad ist die Bernoullische Gleichung im Realtivsystem zu verwenden.) Das 


1 Die Auflésung erfolgt, indem zur Auflésung der einzelnen Integralgleichung (9) angewendete 
Methode ganz formal auf das System (17) iibertragen wird. 
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Quadrat der resultierenden Absolutgeschwindigkeit an der Leitradschaufel ist gegeben durch 


ce? = (uy + Uyy + Uyy + Au,,)? ee ie Cer Urs oe Av,,) > 


wihrend sich fiir das Quadrat der resultierenden Relativgeschwindigkeit an der Laufrad- 
schaufel ergibt 


pie (iy PY, 4 Uys —F © Ty Sin Po + Au»)? + (v9 + ty + ty. — @ Fy COS My + Avis Ness 
(+) Saugseite, (—) Druckseite. 


5. Die Gitterstufe mit endlich dicken Profilen. Unter Beriicksichtigung der in Ziff. 2 bis 4 
gewonnenen Ergebnisse ist es nun nicht mehr schwer, auch noch den Fall einer Gitterstufe mit 
endlich dicken Profilen zu behandeln (Abb. 5). 

Die Profilkonturkurven, auf denen die Zirkulationsverteilungen y, und y, angebracht 
sind, geben wir in der Parameterdarstellung vor 


M, 
Pil) =a, cost,, (4) = Y ey etre, 
y= —M, 
Pa(t,) = — a, cost, + a, + b—a,, ro(tz) = LY cy” ei? 
y= —M, 
(uh es Moby FE) 


Die auch hier vorliegende bereits in Ziff. 4 formulierte Aufgabe fiihrt auf folgende Bedingungs- 
gleichungen, die ein System von Integralgleichungen fir y, und y, darstellen: 


0 = [v% (r,(t,) 2 Yi (t,)) + %y4 (r,(t) ’ Pr(t,)) 
+ ye (r1(t), x(t) ] cos 3,(t,) — [Uy (r,(t,) ’ Pr(ty)) 
+ Uy, (71(ty) ’ Pr(t,)) + Uys (ry(t), Py (t,))] sin 3, (t,), 


(18) 
@ Tp(ty) cos (9,(t,) — Palts)) = [U% (r(t.) ’ Po(to)) ae 04 (75(ts) ’ Pa(to)) 
+ Yye (79(ty) , Pa(to)) J cos @,(t,) — [Up (75(t.) > Po(to)) 
+ Uy (7,(t.) > P2(ty)) + Uy» (r2(t2) ’ Po(to)) ] sin ,(t,) . 
Wie eine Untersuchung analog zu Ziff. 3 zeigt, kann (18) auf die Form gebracht werden 
n = ; EG 
fit) = = i 9, (7) jets — + Ay, (t, » n) dt, + ane | O,(T) Hyy (ty » T) At 
a . (19) 
ci 7 es 
Sots) = 35 i 01(7) Hy (ty ,%) dt + 5 | 0,(T») |etg ~ 5} Ty Fz» (ty - T)| AT « 
Alle H,, sind in beiden Variablen stetig, und fiir sie gilt die Relation 
f Apaltp » 1%) dtp = 0 (p,q =1,2). (20) 


Fiir H,, gilt auBerdem Forme] (13). 
Fiir die Auflésung des Systems (19) mit der Zusatzbedingung (20) wird wieder auf Arbeit I 
verwiesen. Die in den Lésungen o, und o, enthaltenen freien Konstanten 


T= | OG) at. f, = | 6, ,) 4%; 


werden durch die AbfluBbedingungen y, (0) = 0 und y, (0) = 0 festgelegt. Die Berechnung des 
Geschwindigkeitsfeldes ist nun méglich. ; 
Die Ermittelung der Druckverteilung an den Profilen erfolgt prinzipiell in gleicher Weise wie 


vorher (vgl. Ziff. 3 u. 4). 
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6. Der Sonderfall schwach gekriimmter Profile. Liegt ein Gitter mit linienhaften Profilen 
vor, die so schwach gekriimmt sind, daB man in guter Annaherung die Zirkulationsverteilung 
auf der Profilsehne anbringen und die dort induzierte Geschwindigkeit auf die Skelettlinie des 
Profils ibertragen kann, so ergibt sich aus den in Ziff. 2 aufgestellten Formeln an Stelle von (7) 
eine Integralgleichung, die bereits von Kucharski 1 gelést worden ist. 

Im Falle der Gitterstufe wird man mit der Naherungstheorie fiir schwach gekriimmte 
Profile auf das Integralgleichungssystem von Séhngen ? gefiihrt, das geschlossen auflésbar ist. 


7. Zusammenfassung: Es wird die Strémung durch radiale Schaufelgitter nach der Sin- 
gularitatenmethode behandelt. Dabei werden die Zirkulationsverteilungen auf den Profilkon- 
turen selbst angeordnet. Es ist méglich, sowohl Einzelgitter als auch Gitterstufen mit linien- 
haften und mit endlich dicken Profilen nach einer einheitlichen Methode zu behandeln. 


(Eingegangen am 7. September 1953.) 
Anschrift des Verfassers: Dr.-Ing. W. Isay, Berlin-Dahlem, Schweinfurthstr. 90. 


* W. Kucharski, Z. angew. Math. Mech. 21 (1941) S. 65. 
* Vgl. W. Schmeidler, Integralgleichungen Bd. I, S. 56, Leipzig 1950. 
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Die Berechnung der Krafteumlagerungen bei einem einfach 
symmetrischen, in mehreren Lagen vorgespannten Betonbalken 
infolge Kriechen und Schwinden des Betons. 


Von H. Bechert. 


1. Die plastischen Verformungen des Betons. Nach Dischinger verlaufen die plastischen 
Verformungen einer Betonfaser, das von der Belastung hervorgerufene Kriechen und die 
Higenverkiirzung des Betons, gen. Schwinden, nach der Gleichung: 


Oe, __ dp | &scn dp 
Rane a ead” () 


Darin bedeuten ¢ die Betondeformation zur Zeit t, &,., die Schwindverkiirzung des Be- 
tons, p= (t) das Verhaltnis der plastischen zur elastischen Verformung des Betons, yg 
dieses Verhaltnis nach AbschluB der plastischen Verformungen, t die Zeit nach dem Ein- 
leiten der Krafte auf den Balken. 

Setzt man nun die Betondeformation zusammen aus der elastischen Verformung é, zur 
Zeit t = t) und dem dazutretenden plastischen Verformungsbeitrag €, dann ergibt sich aus (1) 
bei Annahme eines in der Zeit und iiber den Querschnitt konstanten Elastizitatsmoduls des 
Betons die Gleichung: 

de, de 


de, de = \dp | Fen dp | 
dt aa dt + (26, Se &b) a Si PE dt Q (Ia) 


2. Herleitung der Gleichungen. Es wird nun angenommen, alle Stahleinlagen seien un- 
mittelbar vor oder nach Einleitung der standigen Lasten, dem Zeitpunkt t = t), mit dem 
Beton in Verbund gebracht worden, so dai bei den nunmehr folgenden Deformationen die 


20) 616" oe. 6 6 (6 6 6: 


4,4 


as 


Ses ky SI 


+@ 


Abb. 1. Verformungen und Krafte am Element infolge der Plastizitat. 


Stahleinlagen jeweils dieselben Verformungen wie die im selben Abstand von der Schwer- 
linie des Betonquerschnitts liegende Betonfaser mitmachen miissen, d.h. es muf sein 


(2) = Sy o 


wenn ¢,; die Stahldeformation der Lage i und (e); die Betondeformation der entsprechenden 
Faser ist. Fir jede Lage haben wir eine solche Gleichung, also insgesamt n Gleichungen. 
Aus (1) mit (Ia) ergibt sich nach einer Umformung 


= Ge (ss, 4 =e 0 (isan (2) 

In der Gleichung (2) sind nun die bekannten GréSen von den unbekannten getrennt. Der 

Term €,.4/Yx ergibt sich aus Erfahrungswerten und wird in der Regel fiir alle Fasern gleich 

groB angenommen. Die (e,); ergeben sich durch Berechnung der Gleichgewichtsgruppe im 

Zeitpunkt t), wobei wir Stauchungen als positiv einfiihren wollen. Dieser Gleichgewichts- 

gruppe wird nun die durch die plastischen Verformungen ausgeléste Kraftegruppe tiberlagert. 
15 
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Die Deformation einer Betonfaser i im Abstand a; ergibt sich zu 
£5; = €s + a; Vb» (3) 


wenn ¢, die Verformung in der Betonschwerachse und y, der Verdrehungswinkel des Be- 
tonquerschnittes bedeuten. Wir definieren weiterhin die Steifigkeiten mit 


K, = E, F, (E, = Betonelastizitaétsmodul), 
S, = By Ji (F, = Betonflache), 
(J, = Betontragheitsmoment), (4) 


(K,,); = (Eu F,); mit E,,; als Stahlelastizitatsmodul der Lage iund F,; als Flache 
derselben Lage. 


Dann fiihren wir die durch die plastischen Verformungen in den Stahlen hervorgerufene 
Krafte X; als Druckkrafte ein. Die Gleichgewichtsbedingungen im Querschnitt liefern dann 


No — M, = — S a, X,. (5) 
k=1 k= 
Setzt man die Gleichungen (4) und (5) in (3) ein, so kommt zusammengefaBt 
dy z it a; a), 
Ce (x, + ea Nee (6) 
k=1 


Definieren wir nun die folgenden Matrizen: 


ie [ Sea + Esch|PE 
X, Ebo2 + Esch| PE 
re aa (7) ae (8) 
Xn \ Ebon i Esch/PE 
1 1 qe @, a, . 1 | 4a \ 
Ka eke ian Ky Ca sian ea Sp 
1 a, A, | 1 1 as J aan 
a Kiem Sir Ky ke ea ee Sb (9) 
1 a, An on Ay Gn } 
Ky ‘Sie Ks ‘Sie 
il aj I a, a, i a, An 
Ks x S§? Ky fs Sb Ks + So 
1 a, Ay 1 a5 1 @, An 
ae Ke Ss ” Kross rat Sp p (10) 
1 a, an 1 Qn, 1 an 
Ria are ea ‘Sie “Kes: 
dann Jautet das System der Differential gleichungen 
d 
Ra 4+ or=3, (II) 


wobei die Matrizen 8 und & symmetrisch sind. ¢ und 8 sind Spaltenvektoren. 


3. Lésung des Systems. Die exakte Auflésung des Systems bereitet theoretisch keine 
Schwierigkeiten. Dagegen ist schon bei mehr als zwei Lagen die Rechenarbeit zur Erlangung 
zahlenmaBiger Werte recht erheblich. Von den drei Schritten zur Lésung des Systems 

a) Berechnung des expliziten Systems, 

b) Berechnung der Eigenwerte der Matrix des expliziten Systems, 

c) Lésung des Anfangswertproblems 


ist durch das folgende, rasch konvergierende Iterationsverfahren nur der erste notig. 
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Die explizite Form von (II) lautet 
dy = aia 
dp —= & 1 8 —— & 1 o t 4 (11) 


wobei die reziproke Matrix St! leicht nach schon bekannten Verfahren ermittelt werden 
kann. Definiert man die Matrizenprodukte mit 


p= R18, (12) 
R21 G = B, (13) 
dann lJautet (11) 
dy aoe 
= ng Ot (11a) 


Die Matrix § hat nur reelle negative Eigenwerte J;, wie aus den technischen Gegebenheiten 
folgt. Ferner folgt daraus 


al< 1. (14) 
Das erweiterte Iterationsverfahren von Picard-Lindeléf liefert dann meist schon nach we- 
nigen Schritten brauchbare Werte. Die Iterationsvorschrift lautet 


Pp 


Emi = | (3 — B tm) dp. (111) 
0 
Das fiihrt auf die Reihe 
12 3 m 
Im = 59 — B35, + B95 —4+..-+ (171 B42. (15) 


Setzt man die Reihe beliebig weit fort, so ergibt sich nach Erweiterung mit § und beider- 
seitiger Subtraktion von € - 4, wobei & die Einheitsmatrix darstellt, die Lésung des Systems 
von Differentialgleichungen zu 


t= (6 — e_ **) 3. (IV) 


Zu dieser Lésung kommt man noch schneller durch formale Integration von (II). Fir die 
praktische Anwendung geniigt bereits nach (15) 


ee po (15a) 
ts ra 67 3 9! aie 3 ont u 2 
Der exakte mathematische Beweis des Verfahrens ist in der Literatur! zu finden. 


4, Beispiel. Gegeben ist der in Abb. 2 skizzierte Plattenbalkenquerschnitt mit den Werten 
F, = 33040 cm*, J, = 76 440 000 cm4, E, = 300 000 kg/cm?, Eo 2,1 ° 10? kelem?, 
Fy = Ferg = 63,6 cm?: F,,3 = 42,4 cm’, a, = 82 cm:a, = 63 cm:a, = 47 cm. 

Ferner betragen die Spannungen in den ein- 
zelnen Fasern bei den einzelnen Lastzusténden 
zur Zeit t= ty 


ee Js 
fea] 


f 

ie 

8S 

a) Eigengewicht : 051 = — 69; O42 > — 535 | 
043 = — 39 kg/cm? (Zug) & 

b) Vorspannung: 0,4; = +149; o43;=—+ 124; t 


053 = +103 kg/cm? (Druck). 
Ferner sei @z = 2,0 und e,,, = 0,0002. 
Daraus ergeben sich dann die 10%-fachen Werte der Matrizen nach 
—230 +497 +100 
3, = (=n) ; Bye (+213 A Spa tae (4-100) ; (8) 
.— 130 +-343 +100 


(i 879; 0,000 334; 0,000 a) 
(ee 


nach (9) 0,000 334; 0,007 759; 0,000 238 


0,000277; 0,000 238; 0,011 435 


1 K. Sattler, Theorie der Verbundkonstruktionen, Berlin 1952; E. Kamke, Differentialgleichungen 
Bd. I, Leipzig 1942, R. Zurmiihl, Matrizen, Berlin 1950. 


By 
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12725203 —5,3867; — 2,9704 
und gat — [—sas 129,1930; — 2358) : 
—=2, 97043 —2,5584; +87,5761 
0,000 392; 0,000 334; 0,000 277 
nach (10) a [0.00 334; 0,000272; 0,000 238 ) 
0,000 277; 0,000238; 0,000 205 


Fiir die explizite Form nach (12) und (13) und (lla) ergibt sich 
—27 928 +60 000 +11 889 
dg = (=a 205 A dv se +-49 802 ° dsch a +12 425 i 
—10 248 271506 + 8204 
0,047 26; 0,040 33; 0,033 36 
a | 


0,040 333, 0,032 73; 0,028 73 
0,022 24; 0,019 15; 0,016 52 


0,004 60; 0,003 86; 0,003 28 
4 == (oo 86; 0,003 25; 0,002 76). 
0,002 19; 0,001 84; 0,001 56 


Nach (15) ergibt sich dann, wenn man nach dem dritten Glied abbricht, da die dazukom- 
menden Glieder schon kleiner sind als 0,5 % des ersten, 


— 55 856 — 5 041 — 324,9 — 51 140 
t= (= 42 1 — (- 4 235) + | — 273,0] = | — 38 628] , 


— 20 498 — 2 396 Le — 18 256 
+ 120 002 +. 11523 + 667,5 + 109 146 
t =(+ 99 603 -(¢ 9 680) + (+ 624.3) ={(+ 90548), 
+ 55011 + 5485 +. 353,7 + 49 880/ 
4099779 +. 2673,4 + 166,6 21273 
Cech =| 24-849) — 4 19045, 8) 4 ao ee eos ae 
+ 16 409 + 1275,8 1 BONO + 15 215,6 


Die einzelnen Anteile am Spannungsabfall im Spannstahl errechnen sich dann zu 


Ao = 6, + dy + Ozer» 


Ao, = — 804 + 1716 + 334 = + 1246 kg/cm? , 
Ao, = — 607 + 1424 + 357 = + 1174 kg/em? , 
Ao, = — 430 + 1176 + 359 = + 1105 kg/cm? . 


Das Vorteilhafte an dieser Berechnungsart ist, da®B man standig mit dem Rechenschieber und 
der Rechenmaschine die einze]lnen Operationen ausfithren kann. 


5. Ubersicht. Von den von Dischinger entwickelten Ansitzen fiir die plastische Ver- 
formung des Beton ausgehend, werden die Gleichungen fiir die Umlagerungskrifte eines vor- 
gespannten, von mehreren Lagen Stahl durchsetzten Betonbalkens angegeben. Sie stellen sich 
in Form eines gekoppelten Systems von Differentialgleichungen dar. Zur Lésung der expli- 
ziten Form dieses Systems wird das duferst rasch konvergierende Naherungsverfahren von 
Picard-Lindeléf beniitzt. 


Fiir die Anregung zu dieser Arbeit méchte ich Herrn Prof. Dr.-Ing. B. Fritz, T.H. Karls- 
ruhe herzlich danken. 


(Eingegangen am 16. September 1953.) 


Anschrift des Verfassers, Dipl.-Ing. Hch. Bechert, Karlsruhe, Lehrstuhl fiir Baustatik der Techn. Hoch- 
schule. 
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Zur kinetischen Bestimmung der Kennlinie eines nichtlinearen 
freien Schwingers. 


Von H. Kauderer. 


1. Aufgabenstellung und Ubersicht. LaGt sich fiir einen freien ungedampften nichtlinearen 
Schwinger auf Grund statischer Messungen der Verlauf der Kennlinie als Funktion des Aus- 
schlags ermitteln, so kann man bekanntlich die bei nichtlinearen Schwingern von der Ampli- 
tude Q abhangige Schwingungsdauer T durch ein Integral darstellen, das sich stets naherungs- 
weise und in gewissen Sonderfallen sogar exakt auswerten 1a8t. Nun gibt es aber auch zahl- 
-reiche Falle (z. B. in der Uhrentechnik), wo man die Kennlinie nur sehr ungenau oder iiber- 
haupt nicht mit statischen Methoden aufnehmen, wohl aber dank der technisch hochentwickel- 
ten Zeitmefigerate die Schwingungsdauer in Abhangigkeit von der Amplitude mit jeder ge- 
wiinschten Genauigkeit messen kann. Méchte mannunaus dieser Ausmessung eines kinetischen 
Vorgangs den Verlauf der Kennlinie bestimmen, so ist das eingangs genannte Problem umzu- 
kehren: Vorgegeben ist die Schwingungsdauer T eines einfachen freien Schwingers als wili- 
kiirliche Funktion seiner Amplitude Q; gesucht ist seine (nichtlineare) Kennlinie. Es soll hier 
versucht werden, einen Beitrag zur allgemeinen Liésung dieses Umkehrproblems zu liefern. 

Wir behandeln zuerst (Ziff. 2) eine fiir unser Problem notwendige Hilfsaufgabe, die auf 
eine Abelsche Integralgleichung fiihrt. Sodann zeigen wir (Ziff. 3), wie man fiir den Fall, daB 
die Schwingungsdauer T als Potenzreihe in Q? vorgeschrieben ist, aus den Koeffizienten dieser 
Potenzreihe die entsprechenden Koeffizienten einer nach Potenzen des Ausschlags fortschrei- 
tenden Reihenentwicklung fiir die Kennlinie gewinnt. Dasselbe Problem fiir den allgemeineren 
Fall, da8 T als Potenzreihe in Q vorgeschrieben ist, 1aBt sich in ahnlicher Weise lésen, webei 
sich fiir die Kennlinie eine Parameterstellung ergibt (Ziff. 4). Es folgen einige Anwendungs- 
beispiele (Ziff. 5) und ein Verfahren zur Nachpriifung der Genauigkeit der Rechenergebnisse, 
die man erhalt, wenn man bei der praktischen Auswertung die Potenzreihen nach einer end- 
lichen Anzahl von Gliedern abbricht (Ziff. 6). Zum SchluB soll noch fiir den Sonderfall, daB 
die Schwingungsdauer als lineare Funktion der Amplitude vorgeschrieben ist, eine geschlossene 
Lésung fiir die Kennlinie in Parameterform hergeleitet werden, die bis jetzt noch nicht be- 
kannt zu sein scheint (Ziff. 7). 

2. Die Abelsche Hilfsaufgabe. Wir gehen davon aus, daB der Ausschlag q des Schwingers 
eine Bewegungsgleichung von der Form 


aq+ R(q)=0 (1) 
befolgen muB, in der a eine positive Konstante (Masse) und R (q) die Riickstellkraft bedeuten. 
Thr Kurvenhbild (Kennlinie) verlaufe symmetrisch zum Null punkt [d. h.es sei R(q) = — R(-q), 
also speziell R(0) = 0], und sie sei fiir positive q selbst positiv. Dann ist ihr Potential 

js gq 
U(q) = {R(q dq (2) 
0 


fiir ¢>0 monoton wachsend, und es existiert somit eine eindeutige positive Umkehrung 
q= 7 (U) dieser Funktion. Mit Hilfe der aus (1) folgenden Energiegleichung 
SF +U(Q)=h, 


in der h =U (Q) das Potential fiir die Amplitude q= Q (die sog. Schwingerenergie) bedeutet, 
14Bt sich dann die Schwingungsdauer T als Funktion der Amplitude durch das uneigentliche 


Integral 


Q 
d 
TQ) = 2720 | pe (3) 
0 


darstellen. Wenn wir in dieser Gleichung mittels der Funktion g = y (U) an Stelle von q das 
Potential U als Integrationsveranderliche cinfiihren, so erhalten wir mit dem nunmehr der 
oberen Grenze Q des Integrals entsprechenden Wert h die Schwingungsdauer als Funktion der 
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Schwingerenergie in der Form 


h dp 


Tih) a dU. 
0 


Denken wir uns nun die Energie h nicht mehr als konstant, sondern als veranderlichen Para- 
meter, und die Schwingungsdauer T als Funktion dieses Parameters vorgeschrieben, so stellt 
die Beziehung eine Abelsche Integralgleichung fiir die Ableitung dg/dU dar, die sich in 
bekannter Weise! nach y (U) auflésen 1aBt. Da q = p (U) wegen (2) fiir U = 0 verschwinden 
mu, erhalt man hierbei den Ausschlag q als Funktion des Potentials U in der Form 


U 
1 T (h) 
Pane a 4) 
0 
Mit der aus (2) folgenden Beziehung 
dU 
= (5) 
wird dann die Riickstellkraft als Funktion des Potentials 
1 
R ( U) 7 dg/dU } (6) 


Die Gleichungen (4) und (6) setzen uns instand, mit Hilfe des Potentials U, das hier die Rolle 
eines Parameters spielt, die Beziehung zwischen Riickstell kraft Rund Ausschlag q herzustellen, 
wenn die Schwingungsdauer T als Funktion der jeweiligen Schwingerenergie h vorgegeben 
ist. 

Dieses im wesentlichen schon auf Abel zuriickgehende Ergebnis soll uns dazu dienen, unsere 
eigentliche Aufgabe, bei der T nicht als Funktion von A sondern von der Amplitude Q vor- 
geschrieben ist, zu lésen. Wir werden diese Lésung nun zuerst fiir den Fall erbringen, da®B T 
durch eine nach ganzen Potenzen von Q? aufsteigende Potenzreihe darstellbar ist. 


3. Die Schwingungsdauer ist durch eine Potenzreihe in Q* vorgeschrieben. Wir haben hier 
vorauszusetzen, dafs sich die Schwingungsdauer T als Funktion von h in einer Taylorschen 


Reihe 


co 1 4 
T (h) = > = FW) (0) h (7) 
v=0 
mit 
Bas 
dh’ 


Te) (0) = ( 


darstellen ]aBt. Gehen wir hiermit in Gleichung (4) ein und integrieren gliedweise, so ergibt 
sich »(U) ebenfalls in Form einer Reihe zu 


p (U) 


, TO (0) = T(0) 


1 


~~ 1 
aaa | 


S12 7a (0) 0 (8) 
Ali 


q| 


1 
Belt 
a}! 
Wir fithren nun eine Funktion 

C(O) =e Oe 


ein und erhalten so mit den ersten fiinf Gliedern der Potenzreihe (8) 


o ~ ee 4 / 8 Vote 16 ‘tt 32 Vtit 2 
CU) = pag UPTON + FTOUT ET’ (OU + 8 re (us4 Zr ousry.]?*, 
Dies stellt, wenn wir ausquadrieren, eine Potenzreihe fir ¢ (U) dar, die sich in der Form 
CU) CU" (9) 
v=1 


1 Vgl. etwa F. Schwank, Randwertprobleme, S. 291, Leipzig 1951. 
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schreiben 148t, wobei die ersten fiinf Koeffizienten C, die Werte 


4 TION, | 
10) T'(0), 

== i FH Ty tee To} | i 
a a 1(0) T'"(0) +42 T(0) To) | 

bs = gra [ons TO) TO) + 35 TO) TO) tae TO} | | 

annehmen. 


Ist nun die Schwingungsdauer T in Abhangigkeit vom Ausschlag durch eine Potenzreihe 
in der Form 


T= % + 120? + ty Q*+- te QO + 4, OF +: =D Tn 02” (11) 
n=0 
vorgegeben, so fithren wir die Rechnung folgendermafen weiter: Da zu q = Q der Wert U=h 
gehort, ersetzen wir in (11) die Potenzen von Q? durch diejenigen von ¢ (h) = Se und 
1 


erhalten so nach Umordnung der Reihenglieder fiir T eine Reihe nach Potenzen von h, aus der 
sich durch Vergleich mit der entsprechenden Reihe (7) die Werte der Ableitungen T (0) in 
den ¢; und den 7; ausdriicken lassen. Wenn wir mit diesen Ableitungen in die Gleichungen (10) 
eingehen, bekommen wir Rekursionsformel]n fiir die €;, aus denen sich diese der Reihe nach in 
den 7; ausdriicken lassen. Hiermit sind die Koeffizienten der Potenzreihe (9) in den gegebenen 
t; dargestellt. Nun kehren wir diese Potenzreihe um und erhalten U(¢) in Form einer nach 
Potenzen von ¢ = q? fortschreitenden Reihe, deren Koeffizienten in den 7; ausgedriickt werden. 
Die Ableitung dieser Potenzreihe nach qliefert dann schlieBlich gemaB (5) die gesuchte Reihen- 
darstellung fiir den Verlauf der Riickstellkraft R(q) = dU/dq. 
Mit Q? —¢ (h) erhalten wir aus (11) 


T (hi) = 3 r4p C(I" = 3 tan( 3 6,¥°) 
n=0 n=0 =) 
oder, wenn wir potenzieren und das Exgebnis nach Potcnzen von h ordnen, 
PCy) "7 Ta Gah (tC go) Cg 2 ty Ce tg Sh) AP A: 
=p [tas + Va (2 510s + C8) 4-3 Teli Co TyCi] A ee, 

woraus sich durch Vergleich mit (7) 

TO) 4,5 

TO) =a, 

T’'(0) = 2 (w0, + 7407); 

T’"(0) = 6 (1203 + 2 T0102 + T C3), 

T'"(0) = 24 [7,0, + ty (20103 + C3) + 8 tell + 567) 


ergibt. Gehen wir hiermit in (10) ein, so folgen die Rekursionsformeln fiir die ¢;: 


se Ae 

b= agg tot 

Ca = sag |ig to 2 (Habe + ta68) + HCH, 

Ca = aa ligg tO (tata + 2 ralsb + 662) + 756, °2 (taba + 7462) ; 
ie sia gas 7 24 [taba + te (2016p +03) + 3 6030, + teC8] + 


— ere OTs Ga ct 2g Sy Go ren) Hae A(t 03 - 2 ty T4070. + gh) 
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deren Auswertung der Reihe nach 


OU=a oy 
1 1 
Co Gaai' 3 70%: 
meet ea : 
b= ap go (T+ 470%) (12) 
i 2 


eyes aay qa 15 (392 t3 + 609 t% tT, + 135 19 Te) » 


1 2 
Gs (x2 a)® 1488 375 


18 (88 396 t4 + 252 672 ty 12% + 107 055 13 7%) % + 
+ 46 935 7277 +18 900 73 7) 


liefert. Wenn wir die Potenzreihe (9) umkehren, so erhalten wir mit € = q 


= 1 q— st sai oon sOUS: g® — 5 68 —5C,6,0,+ C2C, g 
a Te — tt, qe a tease 
‘ 
cy 


Hieraus folgt durch Ableitung nach q die Relgncnrwich inns der Riickstellkraft R(q) in der 


Form 


dU } 
BQ) RG ae Bs tah iis ice eee ta (14) 
wobei die Koeffizienten R,, R;, ...nach (13) in den ¢; und diese nach (12) in den vorgegebenen 
Koeffizienten Ty, T,... der Reihenentwicklung (11) auszudriicken sind. Auf diesem Weg 


ergibt sich nach einiger Zwischenrechnung die Darstellung der ersten fiinf Koeffizienten von 
R(q), die wohl in den meisten Fallen ausreichen werden: 


4x a 
tae 2 2 
Soe : 
he 348 T2 9 | 
2 
Rye 2" (107? 18 tt) 
Rie 2 003 20s 90 
eae a Gee 7575 78 | 3 ToT Tq + 90 Th) » 
647 a : 
Ree 207 615 14 (227 801 t§ — 534 408 t) 1374 + 252180 72 t,t, + 129780 72 72 — 
— 15600 13 7.) . 
Praktisch sind also einfach die gegebenen Koeffizienten 1, tT), tT, ... in die Gleichungen (15) 
einzusetzen, und wir erhalten hiermit unmittelbar die Koeffizienten R,, R;, ... der gesuchten 


Reihenentwicklung (14) der Rickstellkraft nach Potenzen des Ausschlags q. 


_ 4, Die Schwingungsdauer ist durch eine Potenzreihe in Q vorgeschrieben. Um diesen gegen- 
iitber demjenigen von Ziff. 3 al}gemeineren Fall behandeln zu kénnen, miissen wir voraus- 


setzen, daB sich T in Form einer nach Potenzen von Vh ansteigenden Reihe mit den Koeffi- 
zienten t, darstellen lat: 


T (h) =2 ie (16) 


Gehen wir hiermit in Gleichung (4) ein, so ergibt sich g(U) durch gliedweise Integration in 
der Form 


(17) 


TAU Set Lue, 


24200, eal 
9 ! 
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und es folgt gema8B (6) die Riickstellkraft R(U) aus 
; v 


vy—l 
2 


—|! 
Meese ~ (5 
R(U) Toe la (18) 


Nun sei die Schwingungsdauer in der Form - 
Pte 0 + 0? + +: = 20 QO" (19) 


vorgegeben. Setzt man in diese Reihe fiir Q =  (h) die Reihe (17) ein (in der also U durch h zu 
ersetzen ist), so erhalt man wegen (16) die Beziehung 


. ae a 
= 1 = ) ; AEs & »/2 

Wy ee ‘ee 

a 2) 2x0 are = Doth e)) 


Der weitere Gang der Rechnung ist jetzt folgender: Wir rechnen auf der linken Seite dieser 
Gleichung die n-ten Potenzen aus und ordnen die Reihenglieder nach Potenzen von yh . Dann 
lassen sich durch Vergleich gleicher Potenzen von yr in (20) Rekursionsformeln gewinnen, die 
es erméglichen, die Koeffizienten t, aus den gegebenen Koeffizienten Tp, T,, ... nacheinander zu 
berechnen. Mit den so gewonnenen Werten der t, gehen wir in die Gleichungen (17) und (18) 
ein und erhalten so den Zusammenhang zwischen der Riickstellkraft R und dem Ausschlag q 
in Form einer Parameterdarstellung mit U als Parameter. 

Zur Abkiirzung setzen wir voriibergehend 


v 
—)! 
Ve (>) 5 
st = es —s == v 
ie) Tn On ° ( se *): t, Sy, yh uw. ( 1) 
2 . 


Dann geht (20) iitber in 


oder, wenn wir ausmultiplizieren und dabei die Glieder bis > beriicksichtigen, 
Oy + (Sp + 81% + 8%? + 8x8 + 8, x4 + eS) #04 
+ [s§ + 2 59 512 + (2 Sq 8. + sf) x + 2 (59 83 + $1 $2) #8 + °°] 2? Oy + 
+ [8§ + 3 86 81% + 3 59 (Sq S2 + St) 4° + ++°] 29 03 + 
+ (s§6+ 4 s8s,%-+ --:) x404+ 


! 


=) 


v 
Sy 


+ (s8 + ++) x80, “35 
ae 


Der Vergleich der Koeffizienten gleicher Potenzen von x ergibt die Rekursionsformeln 


= 3 
$3 = $204 + 2 Sp $1 02 + $993 » 


to| oe 


e. ! (22) 
pyr ta $3911 (2 Sy 5. + s%) dy + 3 88s, 03 + 5904, 

I « 

= $5 = 840, + 2 (59 83 + $1 82) Op + 3 89 (8 8. + 87) 03 + 4 $6 $104 + 8995+ 


e) 


220 Kauderer: Bestimmung der Kennlinie eines nichtlinearen freien Schwingers. Ingenieur-Archiv 


Setzen wir jeweils in die Gleichung fiir s; (i = 1, 2, ..., 5) die Werte von s;_,,..., S9 aus den 
vorangegangenen Gleichungen ein, so erhalten wir schlieBlich, wenn wir nach (21) wieder zu 
den urspriinglichen Bezeichnungen t, und t, zuriickkehren, 


ty = T> 
: TT 
n¥2a °°’ 
ae 2 
Pees 1 TT TT) 
i oa . 4 ae 
“= aapge Olt Peles xz) ttt t tat , 
1 1 fe ae ; G4 Ue lore (23) 
UST ar eg o|3 7 5} st Ce TTT a 101 T3 + 35 TO Te TF oe TOT > 
-_ 1 1 4. 1 16 LGN ne, 
Ae tots 7 H(3 + | ty t8 ey + (3-4 —-\e8 22 t+ 


DEE Caer 16 oe Ts 
be a) 88 ae a) taut 8+ Tp Tz T3 + 


Mit diesen Werten von t, bilden wir nun gema48 (17) und (18) die Reihen fiir q(U) und 1/R(U). 
Ihre ersten Glieder werden, wenn wir die Reihen in zwei nach ungeraden bzw. geraden Potenzen 


von jU fortschreitende Teile aufspalten, 


q(U) = | YU [ot gal + yg 4l? + | + : u(t } +10 4 putt, (24) 


ae een 4 2s Meee U 3 Laas peda 

Zur praktischen Anwendung der Formeln haben wir somit lediglich aus den gegebenen 
Koeffizienten 7; mittels der Gleichungen (23) die G1éB8en t; zu berechnen und in die Gleichungen 
(24) und (25) einzusetzen. Diese liefern dann die Parameterdarstellung fiir den Zusammen- 
hang zwischen der Riickstellkraft R und dem Ausschlag qg mit dem Potential U als Parameter. 


5. Anwendungsbeispiele. Erstes Beispiel: Die Schwingungsdauer T sei als Funktion 
der Amplitude Q mit den positiven Konstanten T, und Q, in der Form 


T (Q) = T,Goj ¢ 


vorgegeben. Da die @o]-Funktion eine Reihendarstellung nach Potenzen von Q? zulaBt, emp- 
fiehlt sich die Anwendung des Verfahrens von Ziff. 3. Aus der Reihenentwicklung 


TO L a x (S}+ n(Z)'+ 2 | 


ergeben sich die ersten fiinf Koeffizienten von (11) zu 
eG ae gee ae 1 


—_ —? Tt =—- = =? ae 
2 Q3 i 2405 8 72005 


%=T), a 


Gehen wir hiermit in die Gleichungen (15) ein, so erhalten wir 


je ee R,—4A4ve1 _ 1747’ a 1 
ees: = 3S Tee Ot oe Ose 
248 47? a I 1382 4a7%a 1 


2835 T2 Qs’ 
und somit Jautet mit (14) die Potenzreihendarstellung fir die Riickstellkraft 
_ An aQy[ q 4Afaq\?, 17/q\5 248 Ca eels 82 9 
(9) T? e, 3 (3,) + 75 (3) 315 (5) 7 2835 (2,) ue |: 
Die sich hieraus ergebende Kennlinie des Schwingers ist in dimensionsloser Darstellung mit 
& = q/Q), 1 = RTG/4 2? aQ, in Abb. 1 wiedergegeben. Hierbei sind, um den Grad der An- 


naherung durch die ersten Reihenglieder beurteilen zu kénnen, zwei Kurven gezeichnet, Die 
gestrichelte stellt die Naherung dar, welche die mit q? abgebrochene Reihe liefert, wahrend in 
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a ausgezogenen Kurve auch noch das Glied mit q? beriicksichtigt ist. Man erkennt, daBin dem 
ereich, in dem die beiden Kurven praktisch tibereinstimmen, die Kennlinie unterlinear wird, 
was ja bei einer mit der Amplitude monoton wachsenden Schwingungsdauer auch zu erwarten 
war. 


Zweites Beispiel: Die Schwingungsdauer T sei als Funktion der Amplitude Q, wiederum 
mit Konstanten T, und Q,, in der Form 


TO) ee 
(Q) pg 


(26) 


vorgegeben, wobei zunachst Q) positiv sei soll. Da sich diese Funktion durch eine geometrische 


Reihe nach Potenzen von Q entwickcln 148t, kommt hier nur das Verfahren von Ziff. 4 in Frage. 
Aus der Reihendarstellung 


T(Q) = Tot Fat la Bll. | 


Qo 
ergeben sich die Koeffizienten von (19) zu 7 
Cyl ge. by tl Oy 3. os = [QO ++ 
7] 
45 §70;—_ Mo <a al eal 
G | = | 
oN 
xX Sj 
‘| Wo 
Q3 
005 T 
42 + +- 
07 sibs | 
F 
g é 
0 05 40 0 005 070 
Abb. 1. Kennlinie zu T = T, Coj Q . Abb. 2. Kennlinie zu T = tet tee. $ 
5 Qo 1 — (Q/Qo) 


Aus (23) berechnen wir jetzt zuerst die HilfsgiéBen t;. Mit der aus den vorgegebenen Kon- 
stanten a, T, und Q, gebildeten Abkiirzung 


wee 
m V2aQ 
werden diese 
fo = Ty, ty = hTy, ty => (4+) BTy, = 5 (5 +2) Ty, y= Ge (112-4 67x + Int) Ty, 
t= zyp (8056 + 2306 a + 447 n%) IT, 


Gehen wir hiermit in (24) ein, so erhalten wir 


au) = [YO + Atm) MU gy (112 + OT + 98) MUP +) + 
+4 Ul + EO +2 BU + 5p 8056 + 28060 + 447 28) BU + 


oder, wenn wir eine dimensionslose Darstellung mit A = k VU als Parameter und mit & = q/Q, 


beniitzen, 


eae 5 469 g 112 +672+97? 
EQ) =AQL+ DA+ ee Es {pian 50 a+ | 


AY 
3056 + 2306 a + 447 m? 4 
| siete eS 
Sa 4608 wv 


(27) 
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Auf ahnliche Weise folgt dann aus (25) die dimensionslos gemachte Riickstellkraft 
n=—~¥Q,R= Te R/Ar2aQ, Zu 


n (A) = 

is 7 ‘ 6 Be 2 3056 +2306 + 4472? +e 

Pa 4 + 2 5 +2 1124+ 672 si IG ge" se 
eerie ay ge 18 ae 768 A 


Den Zusammenhang zwischen & (A) und 7 (A) und somit ein Bild der dimensionslos gemachten, 
wiederum unterlinearen Kennlinie des Schwingers gibt Abb. 2 wieder. 

Ist Q, in (26) negativ, so ist wegen der Abnahme von T mit wachsendem Q eine iberlinare 
Kennlinie zu erwarten. Da dann k und (bei positivem q und R) auch & und 7 negativ werden, 


erhalten wir die Kennlinie fiir diesen Fall wiederum im ersten Quadranten, wenn wir 7 = k {U 
negative Werte annehmen lassen und in den Gleichungen (27) und (28) bei € und », die jetat 
durch & = q/|Qo| und 1 = 72 R/4 x? alQ,| zu ersetzen sind, die Vorzeichen umkehren. Auch diese 
Kennlinie ist in Abb. 2 eingetragen. 

Drittes Beispiel: Die MeBergebnisse fiir die Schwingungsdauer T der Unruhe einer 
Taschenuhr in Abhangigkeit von der im Bogenma8 gemessenen Amplitude Q lassen sich in 
guter Naherung durch das Polynom 


TF (Q) = % + TQ? + 1% O* + T, O° (29) 
wiedergeben, wobei (in Sekunden) 
T, = 0,4 — 487,8- 10-4, 7, = + 185,46-10-*, 1, = — 8,300 - 10°, 7, = + 0,1056-10°* (30) 


ist und alle 7; mit h6herem Index verschwinden. Auch hier ist wieder, da T(Q) nur gerade Po- 
tenzen von Q enthalt, das Verfahren von Ziff. 3 anzuwenden. Im Unterschied zum ersten Bei- 
spiel ist hier die Abweichung von dem ,,Sollwert* 0,4sek fiir die Schwingungsdauer im ganzen 
Bereich, der sich von Q = 0 bis etwa Q = 4,5 erstreckt, duBerst gering. (Sie spielt jedoch fir 
die Ganggenauigkeit der Uhr trotzdem schon eine erhebliche Rolle.) Wir werden also auch eine 
nur sehr geringe Abweichung der Kennlinie der Unruhefeder von der Geraden zu erwarten 
haben. 

Berechnen wir mit den Zahlwerten (30) die R; gemaB (15) und gehen hiermit in Gleichung 
(14) ein, so ergibt sich das auf das Tragheitsmoment a der Unruhe bezogene Riickstellmoment 
Rin Abhangigkeit von dem im BogenmaB gemessenen Ausschlag q zu 


= R (q) = (247,343 02 q— 0,306 19 g?-+ 0,016 714 g® — 2,6780 - 10-4 q? + 
+ 1,1180 . 10-%q9 — ++ ---) sek-2. 
Da in einem Diagramm fiir R (gq) die geringe Abweichung der Kennlinie von der Geraden nicht 


sichthar wiirde, ist in Abb. 3 nicht = R(q) selbst, sondern diese Ab- 


seks ARQ) 
+7 weichung 
B B I i! 

fee = AR (q) =|— R (q) — 247,34302 a) sek-2 

x wiedergegeben. Die gestrichelte Kurve entspricht wieder der mit 

% dem Glied mit q? abgebrochenen Reihe. 

5 6. Kontrolle der Genauigkeit der Lésung. Da es sich nicht 

ae immer leicht entscheiden 14Bt, welchen Grad der Genauigkeit wir 

-5 durch Abbrechen der Reihenentwicklungen (14) bzw. (24) und (25) 

o nach einem bestimmten Glied schon erreicht haben, diirfte eine 
\ Kontrolle des Ergebnisses erwiinscht sein. Eine solche ist méglich, 

4 eer wenn wir mit Hilfe des berechneten Diagramms der Kennlinie die 

J 


Gleichung (3) auswerten und uns iiberzeugen, ob der Wert von T(Q), 

Abb. 3. Abweichung der Kennlinie dey sich hieraus ergibt, mit dem vorgeschriebenen Wert noch be- 
vom linearen Verlauf bei der neat se ‘ Z P f r 

Unruhe einer Uhr. friedigend iibereinstimmt. Einen Weg, auf dem sich das uneigent- 

liche Integral (3) mit verhaltnismaBig geringem Aufwand berechnen 

1aBt, hat K. Klotter gezeigt, und es sei hier auf dessen Darstellung verwiesen!. Die 


1 K. Klotter, Technische Schwingungslehre Bd. I, S. 152, 2. Aufl. Berlin-Géttingen-Heidelberg 1951. 
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Gleichung von Klotter, der eine Auswertung des geeignet umgeformten Integrals (3) nach der 
Simpsonschen Rege] zugrundeliegt, 1a8t sich, wenn T als Potenzreihe in Q? vorgegeben ist, auch 
noch soweit umgestalten, daB wir eine Forme] fiir T(Q) erhalten, in welche wir die cewae (15) 
berechneten Koeffizienten R,, R;,... der Potenzreihenentwicklung von R(q) direkt einsetzen 
kénnen. Da eine solche Darstellung in mancher Hinsicht bequem erscheint, sei sie hier her- 
geleitet, und zwar gleichzeitig mit einer Verschirfung des Ergebnisses von Klotter, die beson- 
ders dann von Nutzen sein diirfte, wenn, wie im dritten Beispiel, die Kennlinie nae wenig von 
der Geraden abweicht. 

Wir gehen aus von der von Klotter aufgestellten Gleichung fiir T(Q), die in unserer Schreib- 
weise lautet . 


ul 1 1 3 


T(Q)=F12 00 | | | 
: W@R@ Yu@—u(Re) Vu@—uv(ze) ‘Teo —o(e0 
-|- j > (31) 


wobei wieder U das Potential der Riickstellkraft R (q) gemaB (2) bedeutet. Ersetzen wir in (31) 
die unendliche Reihe (14) fiir R(q) durch das tatsachlich berechnete Polynom 

R(q) = Ryq— Rg + R,@—+ °° + (—D* Reng g?**!, (32) 
und entsprechend das Potential U durch das Polynom 


q 
hoe peat 1 1 =1)s 
U(q)= [R@ Bie tee gg, gL ob OE Rang POD, (33) 


so erhalten wir an Stelle von T (Q) den Wert T’(Q), dessen Abweichung von dem vorge- 
schriebenen Wert T(Q) uns ein MaB fiir die Genauigkeit des Ansatzes R(q) fiir die Riick- 
_ stelikraft liefert. Es exgibt sich, wenn wir in (33) fiir die Verwendung in den einze]nen Radi- 


5 
kanden von (31) q= Q, ~ Q, - Q und = Q setzen und die Glieder ordnen, 


T ot. xa ee Pag! I 3 I , 
Cee im VW, "TW, Tw, | Al SS 


wobei wir fiir ein Pol ynom (32), das héchstens vom neunten Grad in q ist, 


W,= R,— BR, + R,A— R,O+ RO, 
W, = 0,060 5469 R, — 0,056 8809 R, Q? + 0,053 5110 R, Q* — 0,050 4101 R, Q% + 
+ 0,047 5540 R, Q°, 
W,, = 0,218 7500 R, — 0,170 8984 R, Q? + 0,137 0036 R, Qt — 0,112 4859 R, Q% + 
+ 0,094 3686 R, Q?, ¢ (35) 
W, = 0,404 2969 R, — 0,240 8409 R, Q? + 0,165 4979 R, Q* — 0,127 8322 R, Q% + 
+ 0,099 9743 R, Q°, 
1 1 1 il 1 
W, —, i — 7 #, Or Bs OC Q° + 79 Be O? 
zu setzen haben. Fiir Pol ynome niedereren Grades werden die entsprechenden Koeffizienten 
Ren+1 einfach gleich Null. 

Wenden wir diese Forme] auf das dritte Beispiel der vorigen Ziffer an, so erhalten wir etwa 
fiir Q? = 10 den Wert T’ (//10) = 0,400 817 sek, wahrend nach (29) T'(/10) = 0,400 642 sek vor- 
geschrieben ist. Dieser in Anbetracht der geringen Abweichung der Kennlinie von der Geraden 
verhalinismaBig groBe Unterschied der beiden Werte riihrt nun aber tatsachlich nicht davon 
her, daB das Naherungspol ynom R (q) fiir q = 10 schon einen schlechten Wert liefern wiirde, 
sondern von der Ungenauigkeit der Forme] (31), die ja nur eine Naherungsformel ist, da sie 
unter Verwendung der Simpsonschen Regel aufgestellt wurde. Es diirfte sich daher gerade fiir 
die Anwendung bei nur wenig von der Geraden abweichenden Kennlinien empfehlen, diese 
Forme] durch ein Korrekturglied zu erganzen. Dieses soll den Fehler, der dieser Formel schon 
bei der Anwendung auf eine gerade Kennlinie anhaftet und der hier den Hauptanteil des Ge- 
samtfehlers ausmachen diirfte, beriicksichtigen. Zu diesem Zweck gehen wir aus von der ge- 
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raden Kennlinie ya 
R(q)= Riq 
und setzen demgem48 in (35) alle Ry,41 auBer R, gleich Null. Dann liefert (34) den Wert 


7" (Q) = 6,285 920 lao const , 
al 


wahrend ja fiir den linearen Schwinger bekanntlich 


@ 
herauskommen miBte. Der Fehler betragt demnach 
6T= T — T= 0,002 7354 ; 
Ry 


Der Wert T’ (Q) nach (34) wird also fiir wenig von der Geraden abweichende Kennlinien ver- 
bessert, wenn wir ihn durch 
T (Q) = T’ (Q) — 6T = T’ (Q) — 0,002 735 /z 
ersetzen. Im Fall unseres Zahlenbeispiels wird so die Korrektur 
6T = 0,000 174 sek 
und daher der verbesserte Wert 
T (10) = (0,400 817 — 0,000 174) sek = 0,400 643 sek; 
dieser stimmt bis auf eine Einheit der sechsten Dezimale mit dem Wert von T (10) iiberein. 
Fir Q? = 20 erhalt man auf demselben Weg T (20) = 0,400 726 sek an Stelle von T (20) 


= 0,400 746 sek. Das Naherungspolynom fiir R (q) liefert also hier ein zu kleines T; d.h. die 
Unruhefeder ist nicht ganz so steif, wie dies die Naherungsformel] angibt, und tatsdchlich wiirde 
auch das nicht mehr beriicksichtigte Glied mit q!! negativ werden. 


7. Eine Lésung in geschlossener Form. Zum Schlu8’B wollen wir noch fiir sich den Fali 
untersuchen, daB die Schwingungsdauer T eine lineare Funktion des Ausschlags q ist: 
T (Q) = T+7Q, 


wobei Tt und zunachst auch T, positiv sein sollen. Gleichung (20) lautet dann wegen 1; = 0 
fiir: = 2; 


0 


Vy 
(eens (2): he Sai 
, heres = Dith ; 
= 
5 ! 


Der Vergleich der Koeffizienten gleicher Potenzen von Vh ergibt t) = T) und fiir »y >1 die 
Rekursionsformel fiir die t, 


ieee 
1 eh Ne 
pS 


Ee st 
"Oo xa ni 
=e 


Ura t (V-=d2ee ye 


Wenden wir diese Forme] y-mal hintereinander an, so ]a4Bt sich t, auf ty = T, zuriickfiihren. 
Es heben sich dann alle Fakultaten bis auf (z)! im Nenner und 0! = 1 im Zahler fort, und 


wir erhalten 


1 T, 
L= U4 Ws 


~ (Baars 
=) 


und diese Gleichung gilt auch fiir y= 0. Gehen wir hiermit in (17) ein, so folgt 
y+1 
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oder, falls wir die dimensionslose GréBe 


Ty 
ge U 
82a 
einfiihren, 
= . y+1 
7 > a 
Ty e=10 ieee) ! 


Die Funktion @ (2) spalten wir jetzt auf in eine Reihe mit ungeraden und eine mit geraden y: 


P (2) = 71 (2) + % (2). 
Hierbei wird 


ess | Ty 
und 
eX 1 ie 
ed ies pees 
A= — 
‘(A+ z)! 
Da 
1 Sb g's 2A+1 
ic 
(A+)! ae: 2 


ist, ]4Bt sich gy, (z) auch in der Form 
W23, ef wl 
92 (2) = eat e+e 4 
; 22 


schreiben. Der Ausdruck in der Klammer rechts ist! nun aber die Reihenentwicklung der 


on] vw 

No] a] ee 
XQ 
No 
+ 

SS ee 


konfluenten hypergeometrischen Funktion M(l, = A ) , so daB 
gy (z) = ets = Ve M(1352)—1| 
Ty x 2 
wird. Um zur Berechnung der Riickstellkraft gema®B (6) die Ableitung 


dp ti dp 


dU 8xadz 
zu bilden, beriicksichtigen wir, daB 


ime) sla) od) 


dpa t: i 4 
dU eae ye eM SA: 
Mittels der Hilfsfunktion 


ist und erhalten 


vse +7 f2M(1, 5.) (8) 


ergibt sich so schlieBlich der Zusammenhang zwischen R= 1/(dy/dU) und q= (2) mit z 
als Parameter in der Form 


iG ; 
(=SyH™—2y, Re—222—1__. (37) 
: y (2) + pee 

Zur numerischen Auswertung sei noch bemerkt, daB fiir groBe z (d.h. praktisch schon 


fiir z = 3), wie aus der Theorie der konfluenten hypergeometrischen Funktionen hervorgeht, 
die asymptotische Darstellung 


eV M(L 34) me (z > 3) (38) 


1 In der Schreibweise von Jahnke-Emde, Tafeln héherer Funktionen, S. 271, 4. Aufl. Leipzig 1948. 
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gilt, so daB man also yp (z) ~ 2e* setzen kann. Da fiir groBe z dann auch im Nenner des Aus- 
drucks (37) fiir R(z) der zweite Summand gegeniiher y (z) vernachlassigbar wird, laBt sich 
y (z) aus den beiden Gleichungen eliminieren, und wir erhalten die Riickstellkraft R direkt 
als Funktion des Ausschlags in der Form 


Sua. ae. Il “ 
R(q) = ee ze Ta) (fiir groBe z). 


ie | a L 


0 5 10 } 20 
Abb. 4a. Kennlinie zn T = 7%, + 7,Q. Verlauf im GroBen. 


Der Verlauf der Riickstellkraft als Funktion von q ist in dimensionsloser Form fir die 
Veranderlichen § = 1, q/t) und 7 = 7, R/8aa in den Abbildungen 4a und 4b dargestellt; 
Abbildung 4b gibt das Diagramm fiir kleine &-Werte wieder. Die Steigung im Nullpunkt ist 
(dy/dé), = 402 

Seither wurde t, als positive Konstante vorausgesetzt. Fir negative t, (also fiir eine 
mit wachsender Amplitude abnehmende Schwingungsdauer) ]48t sich die Rechnung in ana- 
loger Weise durchfithren; wir haben 


9 dann lediglich in (36) und (37) fir Jz 
den negativen Zahlwert einzusetzen. 
8 Wollen wir jedoch auch fir diesen 
7 Fall das positive Vorzeichen von yz 
beibehalten, so sind die Gleichungen 
6 (36) und (37) entsprechend umzu- 
; formen. Sie nehmen dann mit einer 
Hilfsfunktion p (z) an Stelle von yw (z) 
y folgende Gestalt an: 
3 a3 ees 2} m(i 3 
z) =e a pg a 3 
(2) le M (15 
T = 
WS Faery Be) (t, <0) 
7 1 
é = 8ma 1 
0 05 10 15 a) Vole ieet et 
rT parecer tee) (2) 
Abb. 4b. Kennlinie zu T = 1, + 7,Q. Verlauf fir kleine £, Vu Vz 


Mit wachsendem z muB wegen (38) » (z) gegen Null gehen; dann geht aber q (z) gegen 7%/|t,| 
und R(z) gegen Unendlich. Wir erhalten also eine sehr rasch steiler werdende Kennlinie, 
welche die Gerade q = 1,/|t,| zur Asymptote hat. Jenseits dieser Asymptote liegende Punkte 


sind ohne Interesse, da fiir Q> 7/|t,| die Schwingungsdauer bei negativem 1, selbst 
negativ werden miiBte. 


(Eingegangen am 25. September 1953.) 
Anschrift des Verfassers: Dozent Dr. H. Kauderer, Stuttgart-N, Viergiebelweg 11. 
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